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Methods  for  obtaining  iouiratiou  rate  contours,  fractionation 
caters,  atom-concentration  is  toss  it;.-  ratios,  and  particle  site  distribution 
fron  the  Ml  Her  fallout  sod  el  are  given.  Included  are  complete  computer 
pro frass  for  estimating  sublimation  pressures,  ionisation  rate  contours, 
and  particle  site  pa  rase  t  e  rr. 

Graphical  integration  of  activity  iron  fallout  over  the  watersheds 
serving  Houston,  Texas  asd  iTev  'fork  City  is  described.  Maximum  levels  of 
selected  contaminants  at  the  water  intake  for  these  cities  ur.der  nost  adverse 
wind  conditions  were  calculated  and  reported  In  yc/ul  as  follows:  Mew  York 
City;  Sr~=9,  £-3  x  10'3t  Sr-90,  5-9  x  10*5;  Ku-106,  7-1  x  10‘4;  1-131, 

7-6  x  10'2;  Cs-137,  x  10"5;  Ea-li-O,  5-3  x  lO*2.  Houston;  Sr-89,  1.9  x  10‘2j 
Sr-90,  1.9  X  10'4;  Hu- 10c,  2.5  X  10*3;  1-131,  2-5  x  10*1;  Cs-137,  1.0  x  10'4; 
Ba-lk1,  1.9  X  10*1. 

It  Is  concluded  that  the  contribution  of  Induced  radioactivity 
to  the  contf— l.r.stlon  of  water  supplies  would  not  be  significant. 

Many  commercially  available  Instruments  of  the  survey  type  can 
be  used  to  detect  adequately  the  presence  of  substantial  concentrations 
of  radioisotopes  In  water  under  emergency  conditions.  For  more  accurate 
measurement  of  such  concentrations  down  to  safe  levels  more  sensitive 
equipment  or  methods  of  preconcentration  of  the  samples  are  needed.  Such 
methods  or  instrument*  should  extend  the  range  of  sensitivity  by  a  factor 
of  the  order  of  100  and  should  be  reaiily  adaptable  for  field  use.  Outlines 
of  procedures  for  the  analysis  of  radioactivity  in  water  are  given.  Efficiency 
of  water  decontamination  by  various  methods  Is  discussed. 

B oiy  burdens  of  individual  radioisotopes  nay  be  determined  from 
a  consideration  of  the  rate  of  ingestion,  the  effective  decay  rate,  and  the 
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affinity  for  the  Isotope  of  the  critical  organ.  It  Is  heller**}  that 
uptake  of  radioactive  Isotopes  nay  he  reduced  or  prevented  by  selective 
blocking  of  critical  ore&ns  vlth  stable  Isotopes. 
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The  purpose  of  this  study  was  to  evaluate  and  s-.immarize  available 
Information  on  the  problem  of  water  contamination  bv  radioactive  fallout 
in  the  event  of  nuclear  war.  Consideration  has  been  given  to  the  current 
theories  of  the  formation  and  distribution  of  fallout  and  the  level  of 
fallout  that  might  result  from  a  possible  nuclear  attack. 

(l  ? ) 

Of  the  various  attack  models  that  have  beer,  recently  described'  ' 

(d' 

the  or.e  prepared  by  Technical  Operations,  Inc.  '  hn3  been  selected  aa 


most  appropriate  for  the  present  study.  This  model  appears  to  have  seen 


well  thought  out  and  carefully  developed,  giving  due  consideration  to  the 


relative  Importance  of  various  military,  industrial,  governmental,  and 
power  resource  targets.  This  model  assumes  all  detonations  to  be  surface 
bursts,  although  perhaps  a  more  realistic  situation  would  include  many 
air  bursto,  which  would  be  more  effective  against  an  unhardened  military 
or  industrial  target.  As  the  type  of  burst  has  a  considerable  Influence 
on  the  fallout  produced.  It  must  be  carefully  evaluated  In  a  study  at 
water  contamination. 

(M 

A  thorough  analycls  of  the  fallout  model  of  Miller  has  been 
made.  Several  Important  functions  derived  from  this  model  have  been 
utilized  In  this  study  and  are  discuGscd  in  detail  in  Section  II,  "Summary 
of  Analysis  of  the  Fallout  Model". 

Information  from  the  Miller  model  has  ulco  been  used  in  the 
study  of  possible  levels  of  contamination  by  selected  icotopco  In  the  water 
cupply  systems  of  Houston,  Texas  and  Hew  ’fork  City.  With  appropriate  con¬ 
sideration  to  the  assumptions  made,  these  results  may  he  applied  to  the 
watersheds  of  municipalities  other  than  those  shown  here. 
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This  3tudy  Is  corvee  reed  with  the  effect  of  Individual  nuclear 
blasts,  and  does  not  consider  the  effects  of  overlapping  fallout  patterns. 

The  complex  problem  of  overlap  has  been  treated  by  Technical  Operations, 
Inc.^  and  others,  but  additional  studies  on  this  subject  appear  desirable. 

Transport  of  fallout  by  surface  eater  has  been  studied.  Evslua- 
tion  of  the  effects  of  watershed  characteristics  on  the  concentration  of 
radioactivity  in  surface  water  would  be  valuable. 

Decontamination  procedures,  both  standard  and  emergency,  have 
beer,  described  previously .  ^  See  Appendix  A)  These  procedures  are  evaluated 
according  to  their  efficacies  for  the  removal  of  various  radioisotopes  and 
nlfo  according  to  the  feasibility  of  their  use  after  a  nuclear  attack. 

For  purposes  of  this  study,  it  has  been  assumed  that  waterworks  facilities 
have  suffered  only  minimal  damage,  and  that  personnel  will  be  available 
for  operation. 

A  study  of  published  methods  of  radiochemical  analysis  and 
existing  Instrumentation  has  been  made.  Due  to  the  question  of  availability 
or  personnel  and  facilities  it  appear*  that  this  information  will  be  of  more 
value  in  recognizing  long-term  hazards  than  in  determining  potability  la 
the  immediate  post-attack  period. 

‘Ric  biological  uptake  and  resultant  body  burdens  in  man  have  been 
evaluated  in  relation  to  tho  levels  of  water  contamination  considered  likely. 
Existing  deta  haa  been  put  into  a  workable  form  so  *hat  the  body  burden 
of  any  of  the  selected  biologically  important  radioisotopes  may  be  determined 
at  any  specific  time  after  the  initiation  of  intake. 

"The  internal  radiation  hazard  resulting  from  nuclear  attack  is 
of  importance  only  to  those  persons  who  survive  the  catastrophe  for  a 
period  of  several  years.  As  a  primary  onnihilator,  it  is  of  no  consequence. 
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as  the  external  dose  concomitant  with  a  lethal  Internal  dose  la  more  deadly 
by  several  orders  of  magnitude.  In  general,  radiation  from  sub-lethal 
concentrations  of  internally  deposited  radionuclides  act3  slowly,  and 
its  results  do  not  be core  apparent  for  years.  Therefore,  the  internal 
radiation  hazard  is  not  among  the  forces  debilitating  the  population  and 
reducing  its  Immediate  ability  to  cope  with  a  radically  altered  environ¬ 
ment.  Rather  it  has  long-term  effects  which  culminate  in  producing  diseases 
such  as  cancer,  shortening  of  the  life  spaa,  and  possibly  affecting  the 
gem  plasm  of  mankind  to  a  small  degree. 
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II.  SUMMARY  OF  ANALYSIS  OP  TOE  O.C.D.  FALLOUT  MODEL 

A.  Ba3lc  Assumptions  of  the  O.C.D.  Model 

During  the  initial  phase  of  this  project,  the  Postattack  Research 
Division,  Office  of  Civil  Defense,  requested  that  all  studies  be  based  on  the 
fallout  model  as  developid  in  "Fallout  and  Radiological  Countermeasures"^. 

A  thorough  and  quantitative  study  of  the  model  vas  conducted  because  the 
evaluations  of  vater  contamination  depend  directly  on  previously  established 
information  concerning  the  distribution  and  properties  of  local  fallout.  The 
results  of  the  detailed  study  are  now  presented  In  summary  form. 

The  model  was  based  on  both  established  physico-chemical  theories 
and  actual  test  data,  brought  together  by  the  scaling  method-  The  rationality 
and  simplicity  of  tho  basic  assumptions  provide  this  model  with  many  advan¬ 
tages.  The  basic  assumptions  include: 

(1)  The  model  applies  generally  to  a  "ground  surface"  burst,  but 
it  con  be  extended  to  other  types  of  detonations  by  suitable  transformations. 

(2)  The  required  Input  data  consist  of  the  weapon  yield  and  the 
wind  velocity. 

(3)  The  radioactive  decay  schemes,  chemical  properties  and  rela¬ 
tive  abundances  of  the  fission  products  are  based  on  the  data  developed  by 
Bolies  &  Ballou^,  Kntcoff^,  and  Miller  and  Loeb^. 

(4)  The  Increase  of  the  fireball  volume  with  time  Is  beGed  on 
the  adiabatic  expansion  of  an  Ideal  gas  with  the  thermodynamic  equations 
modified  to  Include  a  term  for  the  change  in  free  energy  with  altitude  and 
with  the  external  preooure  proportional  to  exp  (-mgZ/kt)  . 


*  For  explanation  of  symbols,  cce  Glossary. 


(j)  The  radioective  cloud  is  assumed  to  have  the  shape  of  an 
oblate  spheroid;  hence  the  Ionization  rate  contours  or  fallout  patterns 
under  a  constant  wind  velocity  are  generally  cigar -shaped.. 

(6)  All  values  for  the  dimensions  and  properties  of  the  fireball 
and  fallout  patterns  are  scaled  from  test  data  by  the  scaling  method  to  he 
a  function  of  weapon  yield  and  wind  velocity. 

(7)  The  condensation  process  of  fallout  is  divided  into  two 
time  periods:  (a)  the  first  period  is  characterized  by  the  presence  of 
gas  and  liquid  phases  and  cuds  when  the  bulk  carrier  or  substrate  material 
of  the  particles  solidifies,  nnd  (b)  the  second  period  i3  characterized 

by  the  existence  of  ga3  aad  solid  phases.  'Hie  temperature  of  l673°K  Is 
used  to  divide  these  two  periods  of  condensation.  The  radioactive  material 
vhlch  condenses  In  the  first  period  will  be  fused  inelde  the  fallout  parti¬ 
cles  and  generally  becomes  Insoluble  in  water.  The  material  which  condenses 
during  the  second  period  will  be  adsorbed  on  the  outside  of  the  fallout 
particle  and  becomes  readily  soluble  in  water.  With  these  basic  assump¬ 
tions,  the  model  will  yield  the  intensity,  the  activity  and  atom-concentra¬ 
tion  for  a  specific  dowr.wiud  area  and  other,  related  Information. 

B.  Formation  orvl  Ceoactry  of  Stem  and  Cloud 

After  a  surface  burst,  the  fireball  forms  the  shape  of  a  sphere, 
and  ita  radii,  both  horizontal  and  vertical,  expend  exponentially  with 
altitude  as  it  rises.  Since  the  expansion  coefficients  are  different  for 
the  two  radii,  the  rircball  grows  to  t.he  form  of  an  oblate  spheroid,  vith 
a  circular  tor-view  and  an  clllpticul  side-view.  When  the  fireball  reaches 
its  flnul  stabilized  height,  It  is  commonly  known  as  the  "cloud”.  A  stem 
io  shaped  by  the  trajectory  of  a  continuous!;/  expanding  fireball  and 
assumes  the  fora  of  an  inverted  exponential  horn. 
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The  Initial  spherical  radius  of  the  fireball,  R  ,  the  final 

s 

horizontal  semi-axis,  a,  the  final  vertical  semi-axis,  b,  and  the  final 
height  of  the  center,  h,  of  the  cloud  ■  ere  related  to  the  weapon  yield, 
V,  through  empirical  data  as  shown  In  the  following  sealing  functions: 


2.09  x  102  W0,333  ft. 

V 

•  1  to  105  KT 

(1) 

2.45  X  103  V0-1*31  ft. 

V 

»  1  to  105  KT 

(2) 

1.40  x  103  w0*300  ft. 

V 

-  1  to  105  KT 

(3) 

0.66  x  iou  v0"1*1*5  ft. 

V 

-  1  to  28  KT 

<«0 

1.68  x  10U  V0,164  ft. 

V 

-  28  to  105  KT 

(4a) 

At  a  given  altitude,  Z,  the  horizontal  and  vertical  serai-axle 
of  the  fireball  can  be  obtained  by 


*z  -  V 


\z 


\  -  V 


V 


(5) 

(6) 


where  a  ,  to  ,  x  ,  X.  can  be  solved  from  the  known  information  of  R  ,  a,  b,  h. 

O  O  ft  0  0 

The  geometry  of  both  stem  end  cloud  is  shown  In  Figure  1. 

C.  Particle  Size  Paraj-eter  and  Ionization  Rate  Contours 

As  mentioned  previously  under  the  basic  assumptions,  the  fission 
products  arc  either  fused  Inside  or  adsorbed  on  the  outside  of  the  dust 
particles  and  are  Inhaled  into  the  fireball.  The  ultimate  ionization  rate  • 
contours,  or  fallout  patterns,  on  the  ground  are  related  to  the  size  of  the 
falling  particles.  The  particle  size  parameter,  a,  is  defined  as  the  ratio 
of  wind  velocity,  V  ,  to  the  average  falling  velocity  of  the  particle, 
since  It  la  more  concerned  with  the  falling  rate  than  the  actual  physical 


size  of  a  particle 


9 


For  the  9tea,  It  is  assumed  that  the  particles  at  a  given  altitude  in 

the  rising  fireball  are  concentrated  on  a  horizontal  plane  through  the  center 

of  the  fireball.  Therefore,  all  particles  which  have  equal  values  of  a  fall 

out  of  the  stem  starting  at  the  same  altitude.  Any  particle  starts  to  fail 

when  Its  instantaneous  falling  rate,  V^,  equals  the  rate  of  rise  of  the  our* 
d2  _ 

rounding  air  m03S,  tt*  There  i3  an  empirical  relation  bctueen  V  ,  V  ,  and 

UU  Ct  l 

the  flrebalL  altitude  Z: 

V2/Vf  -  p  ♦  qZ  (7) 

vhere  p  •  0-95.  q  ■  1.02  x  10*^,  where  Z  ■  5»C00  to  50,000-  ft 

with  particle  diameter  d  ■  200  to  1,200  microns 
p  -  0.58,  q  -  1-71*  x  10*5,  when  Z  -  50,000  to  110,000  It 
with  particle  diameter  d  -  300  to  1,000  microns 


There  io  also  an  empirical  approximation  for  Z: 


•M 

Z0  (1  -  e  ^ 


for  t  »  20  to  •‘BO  seconda  (B) 


where  ZQ  lo  a  yield  dependent  multiplier 


-1 


>2  io  0.011  sec  . 


Itous, 


Z  '  g  -  >Z  <*o  -  2> 


(9) 


For  values  of  Z 


see  Figure  5- 


■IQ 


Therefore,  the  particle  site  parameter,  a,  for  the  particle  group 
falling  free  the  altitude  Z  la 


\  Vjp  ♦  q2) 

Vf  '  V2o  -_Zl 


'  (10) 


Moreover,  it  can  be  determined  from  geonetry  of  the  stem  that 


-  X  *  a^  (  ♦  for  a 


min* 


for  a  ) 
max 


(11) 


Therefore,  Q _ and  a  cay  be  obtained  for  the  particles  falling 

QLELX  mill 

at  a  distance  X  from  ground  zero.  In  addition,  the  rising  and  falling  time 
for  particles  of  a  given  value  of  a  can.  bo  determined  from  combining  Equa¬ 
tions  (8)  and  (10),  as  shovn  below 


*r  (sec) 


1_ 

*2 


In  ( 


■Vo  *  V*o 


j* 


tf  (oec) 


°(°V0  *  V} 

Vv(«^  ♦  Vyq7 


(12) 

(13) 


The  outer  dimensions  of  the  cloud  are  defined  by 


(14) 


Consequent ly,  and  a  n  for  the  downwind  distance  X  con  be 

derived,  as  shovn  by  the  equation 


Q 

m 


hX 


t  y? 


LL 


Va2) 


«  2  **  \  t  . 

(a  -  V  )  [h  ~  b  (1  - 


u 


h2  -  b2(l  -  y 


/a2) 


(15) 


This  calculation  usually  includes  the  l8o  sec.  overage  delay  time  due 
to  fireball  toroidal  circulation,  which  should  be  subtracted  in  order 
to  obtain  the  actual  rising  time. 


-11- 


(where  the  sign  Is  ♦  for  aB&x>  611(1  *  for  aniin) 

Since  the  cloud  is  assumed  to  he  stabilized,  all  particles  with 
the  same  values  of  a  fall  along  the  same  slope.  Hence,  the  tines  of  arrival 
and  cessation  for  particles  of  a  given  a  can  be  determined  from  the  geometry 
of  the  cloud 


ta  •  «  (h  ♦  za)  /  Vy 


(16) 


tc  -  «  (h  ♦  zc)  /  \ 


(17) 


vhere  z  and  z  are  the  smallest  and  largest  intercepts  of  the  slope  to 

a  c 

the  outer  dimension  of  the  cloud,  respectively. 


QtX  -  Oh)  b2  ♦  ab  J  (a2  *  a2!?2)  (1  -  y2/b2)  -  fx  -  oil)7 

«2  .  A2 


Z  m 


(19) 
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(where  the  sign  is  +  for  i  ,  and  -  for  z  ). 

c  a 

Assuming  that  the  ionization  rate,  or  fallout  intensity,  1^, 
vnrLeii  exponentially  with  downwind  distcnee ,  end  its  dlstrltortion  is  closely 
related  to  a  and  W,  also  confirmed  by  cr.pl r leal  test  data,  a  complete  set 
of  scaling  functions  which  defines  the  voriatioc  of  intensity  versus  down¬ 
wind  distance  and  the  bi-elliptical  outline  of  the  ionization  rate  contours 
is  tabulated  In  Chapter  3  of  "Fallout  and  Radiological  Countermeasures"^^ . 
Those  scaling  functions  arc  not  repeated  in  this  report,  but  a  corresponding 
computer  program  was  established,  as  presented  in  Figure  5.  A  typical 
ionization  rate  contour  for  a  10 MT  weapon  yield  and  a  15  mph  wind  velocity 
is  shown  in  Figure  6.  Since  the  only  inputs  are  weapon  yield  and  wind 
velocity,  this  computer  program  is  relatively  simple  to  understand  and 
to  use  by  others  to  obtain  the  information  on  ionization  rate  contours. 

It  should  be  pointed  out  that  these  contours  have  been  corrected 
to  a  standard  reference  time,  H  ♦  1  hour,  by  a  decay  correction  factor 


I  (t:  actual  time) 
d(t,  l)  »  jx(1:  «  ♦  1  hour) 


(19) 


This  decay  correction  factor  can  be  obtained  directly  from  a 
typical  close-in  fallout  decay  curve  ns  shown  in  Figure  2  , 

Knowing  the  ionization  rate,  the  exposure  dose,  D^,  in  the  cloud 
fallout  area,  which  is  of  greater  Interest  than  the  stem  fallout  a*-ea,  can 
te  estimated  by 


DK  */t  ¥t>dt 


(20) 


(u  os/noissu  a3d  bH  a  5 ( _oi >  3iva  NonvziNOt 


D.  Condensation  of  Fallout  or.d  Practlonatlon  tfumbe ra 


During  a  nuclear  detonation,  a  tremendous  amount  of  heat  Is  created 
vy  nuclear  fission  and  fusion.  This  heat  melts  and  also  vaporises  every¬ 
thing  inhaled  Into  the  air  mass  fron  the  ground  and  forma  the  familiar 
fireball.  When  the  fireball  rises,  Its  volume  begins  to  expand,  and  lta 
temperature  io  assumed  to  drop  according  to  the  Ideal  gos  lav.  The  material 
swept  Inside  the  fireball  starts  to  liquldlfy  and  finally  solidify  to 
particles,  in  or  on  which  the  fission  products  begin  to  condense.  The 
condensation  process  Is  generally  divided  Into  two  time  periods  as  mentioned 
In  the  basic  assumptions.  The  dividing  temperature  depends  largely  on  the 
composition  of  the  fireball.  For  a  model  surface  burst  on  a  soil  consisting 
of  silicate  minerals  this  temperature  Is  assumed  to  be  l6?3°K,  as  cost  ele¬ 
ments  will  have  started  to  solidify  e.t  this  temperature.  Since  the  two 
periods  are  charactarlzed  by  the  coexistence  of  either  gas  and  liquid,  or 
gas  and  solid,  the  vaporization  or  sublimation  pressure  plays  an  important 
role  In  the  condensation  process.  The  following  formula  for  the  evalua¬ 
tion  of  the  pressure  la  presented  In  "Fallout  and  Radiological  Counter- 
measures  '  : 


log  Pj  -  j  ♦  B  ♦  G  (21) 

where  Pj  Is  the  vaporization  or  sublimation  pressure  for  element  J, 
gm/ sq .  cm . , 

T  lo  temperature  of  the  fireball,  °K, 

A,  B,  C  are  the  empirical  constants  for  vaporization  or  sublimation 

reactions  of  elements  contained  in  the  fireball  which  were 

(4) 

tabulated  In  Vdller'a  manuscript  . 

.  The  established  computer  program  for  estimating  Pj  Is  shown  In 
Figure  3» 


Figure  3.  Cooputer  Program  for  Estimating  Sublimation  Pressure 
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Figure  3  (Continued).  Computer  Program  for  Estimating  Sublimation  Pressure 
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Figure  3  (Continued).  Computer  Program  for  Estimating  Sublimation  Pressure 
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Wgure  3  (Continued).  Omputer  Program  for  Katl  mating  Subllmticn  Pres*ur» 
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figure  3  (Osutlnued).  Cbaputer  Program  for  Estimating  Sublimation  Pro*sure 
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*.»4# 

.040404 

.11072?*  *• 

2 

JO** 

9 

*444 

.«*»!, 

>.»** 

•  8404 J9 

.3*7794,  if 

4 

"0 

10*4 

• 

*444 

•  11*4, 

.41*7)1,  «* 

t 

"0 

}<t** 

4 

*0«* 

10)1. 

•2.041 

.000004 

, )7004», -91 

»0 

1040 

*.»*•* 

•29*1.0. 

0.0*0 

.000004 

.3391 14,-9* 

"0 

to** 

4 

7(*(H 

•72**0. 

0.479 

..}17*»1 

,1>(«I7»-0I 

« 

HO 

J*** 

*444 

•40*, 

..  1.0*7 

-.071974 

.,10704*,  ll 

a 

"0 

1**4 

*444 

11*4*. 

..0*1 

.1.07*4** 

.l"77*t,  *9 

4 

HO 

10*4 

*0*4 

»mo. 

.4.1»* 

-1.1*74*4 

.l>t«l>,  J* 

5 

HO 

10*4 

*044 

*44*0 , 

-7.7** 

-1.7774?* 

.  .  .17*779,  »■  ... 

TC 

1*00 

*444 

.41*04. 

12.0*7 

.1404*1 

,*47*74.-1*  . 

TC 

1*44 

*444 

.*11*0. 

9.004 

.0«0404 

,377741,-9*  . 

TC 

44.1 

Ml 

.101*. 

1,2?* 

, *0*4*4  _  . 

./'HI**  9*  - 

cu 

*04 

131* 

.1741*. 

0.414 

.00*004 

.17*121,-9* 

cu 

1-15* 

*4*4 

.1*1*4. 

».*•! 

.0000*4 

.9437#»,-a4 

•  0 

*44 

1>J* 

.14*)*, 

0.1*4 

,00*4*4 

. ,3*210*, -If  .  . 

1?J« 

jam 

.11*0*. 

».0*1 

.0004*4 

.27*0*1, -Of 

1440 

14*1 

•*9011. 

0.10* 

.0094*4 

.1301*0,-1* 

fo 

1*  JJ 

*144 

•  1910*. 

4.0*9 

.0000*4 

.111077, -1» 

1*M 

1*24 

.4*10. 

4.240 

.0004*4 

.14)1*4,  9* 

**« 

41* 

*044 

•1**0, 

4.0»0 

• *0*0*0 

.274**3,  If 
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The  elements  In  the  firebaiL  are  radioactive  ond  hence  the  number 
of  moles  of  each  species  is  constantly  changing  with  time.  If  y,,(t) 
denotes  the  amount  of  radionuclides  of  element  J  ar.d  mass  number  A  present 
at  time  t  after  fission,  then  the  total  amount  of  element  J  present  at  tine  t 


'fj(t)  «  atoris  or  ®oles 


Correspondingly,  the  total  amount  for  the  entire  chain  yield  of 
mass  number  A  Is 


YJA(t)  atoms  or  moles 


in  which  16  constant,  except  for  masu  chains  containing  neutron  emitters. 

For  the  first  period  of  condencation,  the  material  balance  for 
element  J  between  gas  and  liquid  phase  la 

ty°  ■  Vj*(,)  •  V;a°<1>  *  Vja<‘>  (2‘ 

where  njA°(t)  la  the  number  of  moles  of  element  j  with  mass  number  A, 
which  is  mixed  with  moles  from  other  mass  chains  to 
fora  n  moles  of  vapor,  and  iu  abbreviated  to  n®,  and 
r.  .(t)  is  the  number  of  moico  or  eicment  J  with  msec  number  A 

JA 

dissolved  in  the  n(/)  moles  of  liquid  carrier,  which 
Is  the  particle  in  the  liquid  phase,  prior  to  solidi¬ 
fication,  and  is  abbreviated  to  n.. 

J 

For  a  p.-rftct  gas,  the  ratio  between  the  mole  fraction  of  clement  J 
in  the  vupor  pnasc  and  that  in  the  ,’iquld  phase  is 


-2L- 


N?  n®/n  n®  KT/p  V  k 

Nj  "  nj/n(/J  *  *  p^  (25) 


where 


n 


n(i) 


R 


T 


PJ 

V 


is  the  cole  fraction  of  element  J  in  the  vapor  phase 

is  the  cole  fraction  of  element  J  in  the  liquid  phase 

is  the  number  of  moles  of  vapor 

is  the  total  moles  of  liquid  carrier 

Is  the  nolar  gas  constant 

is  the  absolute  temperature 

is  the  total  vaporization  pressure 

is  tiie  molar  volume,  and 

is  the  Henry's  Law  constant 


From  Equation  (25),  the  following  relation  between  n®  and  n  can  oe 

J  J 


obtained: 


n.k. 


NO/ 


VI  RT 


"jkj 


(26) 


where  k°  »  k^/[n(/)/vj  RT 


Substituting  Equation  (26)  into  Equation  (24)  yields 

■  <‘  *  9  Vj *<'> 


(27) 


The  fractlor.atlon  number  of  the  first  period  of  condensation, 

r  (A,t),  is  defined  to  be  the  fraction  of  the  moos  chain  condensed,  or  in 
o 

mathematical  fora 


rQ(A,  t)  -  ^  Ej  nJA(0 


y.1A(t) 
l  ♦  k? 


(28) 


vh*  re 


where  t,  in  seconds,  is  related  to  the  fireball  temperature,  T,  by 

T  ■  U .66  x  103  w’°*°l°  exp( >0.546  «*°*373t)°K  (33) 


and  the  data  of  yjA(t)  verc  calculated  and  tabulated  as  Hj(A,  t)  by  holies 
and  Ballou^  .  ;ter:cc,  the  fractionation  number,  ro(A,  t),  can  be  evaluate  i,r 
the  elements  contained  in  the  fireball. 

During  the  second  period  of  condensation,  if  an  exceed  of  solid 
surface  area  iu  present,  the  reiml  condensed  at  r.ny  time  after  colldlflea- 
tior.  of  the  carrier,  io  given  by 


-23' 


ai  “  ai  ’  °j 


(5«0 


where  o.'  Is  the  amount  of  element  J  condensed  on  the  surface  of 

the  solid  particles 

n?  is  the  amount  of  element  J  that  has  not  condensed  in  the 

V 

liquid  carrier,  and 

n^j  is  the  amount  of  element  J  In  the  vapor  phase. 


For  a  perfect  gas 


ni  ‘  iar 


(35) 


where 


1 3  the  sublimation  pressure  and  may  be  evaluated  according 
to  Equation  (2l). 


The  material  balance  for  element  J  is 


7j(t)  •  Dj(t)  ♦  nj(t)  ♦  n”(t) 


(36) 


Therefore,  employing  Equations  (27)  and  (35) 

.  £ 

RT 


(37) 


or 


VjA 


1  ♦  K 


ZAyJA  ’  K?  EAP«A 


(38) 


In  thla  equation  the  partial  pressure  of  each  nuclide  is  propor¬ 
tional  to  Its  abundance  at  uny  given  time,  ar.d  is  given  by 


■  <x 


(39) 
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‘  (I  ’  ro)yJA* 

vO  -  o 
■J  *  Zkyik 


The  f  rac  t  i  on*.  1 1  on  number  of  the  second  period  of  condensation 
r^(A,  t)  Is  defined  to  be  the  fraction  of  the  mass  chain  condensed  up  to 
a  given  time  during  the  s-.cond  period,  or  expressed  mathematically 


r;  (A,  t) 


h 


V  r  a  i 
0.23  SVRT  JPJAJ 


(40) 


(A  t\  -  \  „  ^  /  A  f  \  m  - - - 

‘o  lA’  }  1  O’’  z;  O^jtWK? 


£/3A/y>°A 

EJ  yjA 


(41) 


vhere 


is  the  fireball  volume  and  is  given  by  scaling  the  function 
V  -  7.62  x  1011  W  exp  (0.510  W*°'373t)  cmJ 


B  le  the  ratio  of  fl66lon  to  total  yield,  and 

0.23  Is  a  conversion  constant  changing  weapon  yield,  W,  to 

r.urter  of  fission-moles  per  atom. 

Physically  In  Equation  (-1),  p1?  V/RT  glveo  the  number  of  molts  of  clement  J 

J 

In  the  vapor  phase  sr.d  0.23  Wl?  gives  the  total  n'lmter  of  moles  of  the 

V 

element  produced  nr.d  In  existence  at  a  given  time  In  the  second  period. 

The  combined  tern  in  Elation  (4l)  gives  the  fraction  still  remaining  In 
the  vapor  -tats  during  the  second  period  of  condensation.  The  1  -  rQ(A,  t) 
torn  glv<  j  the  fraction  which  has  rot.  condensed  in  the  liquid  carriers  during 


their  existence  In  the 


rot  period  of  condensation. 
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The  fractionation  numbers  yield  the  information  about  the  fractions 
of  radioactivity  condensed  in,  or  on  the  particles,  at  any  tine  after  a  burst. 
Therefore,  the  fractionation  numbers  not  only  make  the  conversion  between 
Intensity  and  activity  possible,  but  also  lead  the  vay  to  the  evaluation  of 
soluble  and  Insoluble  activity  vhich  is  important  to  biological  availability 
and  radiological  counte  measures  research. 

S.  Conversion  between  Intensity  and  Activity 

The  conversion  between  Intensity  (measured  as  ionization  rate)  and 
activity  density  in  fissions  per  unit  area  is  presented  in  "Fallout  and 
Radiological  Countermeasures",^**)  based  on  the  simple  assumption  that 


Ix(t) 


Ix(t)(ideal) 

ideal) 


yt) 


o*> 


where  the  conversior  factor,  K^t),  *s  further  broken  down  into  the  fol¬ 
lowing  components 

K^t)  -  Dqx  (rx(t)  ifp(t)  ♦  l1(t)l 


(>3) 


in  which 


rx(0 


it(t) 


is  an  instrument  response  factor,  usually  assigned  a 
value  of  0.75 

is  the  terrain  shielding  factor,  usually  assigned  a 
value  of  0.75 

l 8  the  gross  fission  product  fraction  number,  sometiaea 
aloo  designated  as  rQ(t) 

is  the  air  ionization  rate  per  fission  at  3  ft  above 
an  Infinite,  ideal  plane  for  a  uniform  distribution  or 
the  normal  fission  product  mixture,  and 
Is  the  cone  unit  for  neutron  Induced  activities. 
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Details  about  the  neutron  Induced  ionization  rate,  i^(t),  are  almost  entirely 
classified,  but  generally  assigned^  as  1,9$  of  i^p(t).  The  other  components 
In  equation  (^3)  can  be  readily  derived. 

•  Detailed  information  about  the  activities  of  the  U-235  fission 

(7) 

products  has  been  provided  by  Bolles  and  Ballou  .  These  data  can  be 

(8) 

correlated  to  the  data  supplied  by  Katcoff'  and  applied  to  other  types 
of  fission  elements,  such  as  U-238,  Pu-239  by  taking  proportions  of  the 
corresponding  mass  number.  This  conversion  has  been  systematized  and 
tabulated  in  "Fallout  and  Radiological  Countermeasures".  Moreover,  Hiller^14  ^ 
translated  the  activity  data  from  Bolles  and  Ballou  into  Ionization  rate  by 
a  disintegration  multiplier,  m.  Both  tables  of  correlation  and  of  disinte¬ 
gration  multipliers  are  reproduced  from  Miller's  manuscript  as  Tables  I 
and  II.  Therefore,  the  ideal  or  normal  ionization  rate  of  fission  product 
mixture  may  be  computed  according  to  the  following  formula: 

ifp(t)  -  Bn At  (kk) 

4 

where  Afc  is  the  activity  of  the  nuclide  at  time,  t,  per  10  fission. 

In  accordance  with  the  earlier  discussion,  the  fission  products 
do  condense  and  fractionate  with  time  instead  of  condensing  completely  at 
the  very  beginning  as  in  the  ideal  situation.  Therefore,  the  real  ioniza¬ 
tion  rate,  if*(t),  should  be  multiplied  individually  for  each  radionuclide 
by  its  fractionation  number  at  any  specified  time.  Mathematically,  this 
may  be  expressed  by 

if*(t)  -  Bn  At(rQ(A,  t)  ♦  r^(A,  t)J  (Uj) 

The  (A,  t)  term  vanishes  in  the  first  period  of  condensation. 

As  time  Increases,  the  sum  of  the  two  fractionation  numbers  approaches  unity. 


i 
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TABUS  I 

(4) 

Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  la  per  cent  of  fissions} 


Mass 

U2^ 

U2^ 

Pu2^ 

number 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Pi  33  Ion 
Neutrons 

72 

1.6xlO*5 

4.6xl0’U 

5.0xl0"6 

m 

•4* 

1.2x10 

• 

73 

-4 

l.LxlO 

0.0012 

3.7xl0”5 

« 

2 .2xl0*U 

- 

74 

(3-2xlO*V 

0.0034 

1.1x10 

0.001 

4.1xlO~4 

0.0011 

75 

(8.8x10^) 

0.0062 

8.3x10*^ 

0.0040 

J4- 

7.xl0 

0.0023 

76 

(0.0029) 

0.012 

0.0012 

0.0078 

0.0014 

0.0051 

77 

O.OO83 

0.023 

0.0038* 

0.014 

0.0026 

0.011 

78 

C.021 

0.048 

0.0095 

0.026 

0.0049 

0.025 

79 

(0.041) 

0.096 

0.019 

0.053 

0.0090 

0.043 

80 

(0.077) 

0.19 

0,045 

0.096 

0.016 

0.075 

8l 

0.14 

0.21 

0.088 

0.18 

0.030 

0.14 

82 

(0.29) 

0.50 

0.20 

0.35 

0.056 

0.23 

83 

0.544 

0.80 

0.40* 

0.66 

0.10 

0.37 

84 

1.00 

1.3 

0.85 

1.02 

0.17 

0.60 

85 

1.30 

1.85 

0.80 

1.45 

0.28 

0.92 

86 

2.02 

2.5 

b+ 

WJ 

CD 

+ 

1.9 

0.45 

1.15 

87 

(2.94) 

3-3 

l.S» 

ro 

0.73 

1-5 

88 

(3.92) 

4.2 

2.45 

2.7 

1.2 

1-9 

89 

4.79 

5.1 

2.9* 

3.17 

1.9* 

£  .  ^ 

Continued 


TA3LE  I  (continued) 


C-JMulatlve  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  per  cent  of  fissions) 


y235  -j23S 


Thermal 

lieutrons* 

Flsslor. 

2,’eutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Nc-utrcas 

Fission 
Neut rons 

20 

5.77 

5.8 

3.2* 

3.7 

2.4 

3-0 

21 

5-e\ 

3.85 

3.6 

4.3 

3-0 

3-7 

9£ 

6.03 

6.0 

4.1 

4.8 

3.7 

4.4 

93 

6.45 

6.4 

4  >85 

5.2 

4.6 

5.0 

94 

6.40 

6.4 

5.3 

5-45 

5.5 

5A 

93 

6.27 

6.3 

5-7* 

5-6 

5.9* 

5.6 

9^ 

6.33 

6.3 

5.8 

5-7 

5-7 

5.3 

97 

6.0? 

6.1 

5.7 

5.64 

5-6* 

5.2* 

9* 

5.78 

5.8 

5-7 

5.6 

5.4 

5.4 

92 

6.06 

6.1*-* 

6.3* 

6.2*» 

5*9* 

5-9* 

100 

6.30 

6.7 

6.1 

6.4 

6.0 

6.4 

101 

5.0 

5-3 

5.5 

6.5 

6.0 

5-9 

102 

4.1 

2.9 

5.6 

5-9 

5.9 

5.3 

103 

3.0 

1.7 

6.6 

5.0 

5.8* 

4.6 

104 

1.8 

0.95 

5.4 

3.2 

5.0 

3.5 

102 

0.90 

0.54 

3.9 

2.2 

3.9* 

3.2 

106 

0.  j3 

O.jO 

2,7* 

1.5 

5.0* 

3.6 

107 

0.19 

0.17 

i.35  ' 

1.0 

4.0 

3.1 

103 

(0.035) 

O.095 

0.67 

0.70 

3.0 

2.6 

109 

(0.0;9) 

0.0>3»** 

0.32* 

0.43 

1.5* 

1.9* 

110 

(0.00) 

0.030 

0.15 

0.33 

0.65 

0.81 

111 

(0.015) 

0.0..2*** 

0.0/3* 

0.23'** 

0.27* 

O.34 

li: 

(0.013) 

0,0_0*»» 

0.046* 

0.19 

0.10* 

0.14* 

113 

(0.0L?) 

0.018 

0.043 

0.17 

0.055 

O.OpO 

.29' 


TABLE  i  (continued) 

Cumulative  Mass-Chain  Yields  of  Fission  Products 


(Values  are 

in  per  cent  of  fissions) 

Mass 

Number 

U2*5 

U2^ 

Pu2*> 

Thermal^ 

Neutrons 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thi-rrrifll 

Neutrons 

Ficslon 

Ne-tron; 

114 

(O.Oll) 

0.017 

0.041 

0.16 

0.046 

0.075 

U5 

0.0104 

0.017"** 

0.040* 

0.15*»* 

0.041* 

0.069* 

116 

(o.oio)b 

0.017b 

0.039 

0.14 

0.039 

O.O65 

117 

(0.010) 

0.017 

0.039 

Q.l4b 

O.O38 

O.C65 

118 

(0.010) 

0.017 

0.4pb 

0.14 

O.O^ 

0.064b 

119 

(o.ou) 

0.017 

0.041 

0.14 

0.039 

0.064 

120 

(o.ou) 

0.018 

0.042 

0.15 

0.041 

0.065 

121 

(0.012) 

0.020 

0.044 

0.16 

0.044* 

0.066 

122 

(0.013) 

0.022 

0.046 

0.17 

0.047 

0.069 

123 

(0.015) 

0.030 

0.050 

0.19 

0.052 

0.076 

124 

(0.017) 

0.053 

0.055 

O.23 

O.O58 

0.082 

125 

0.021 

0.095 

0.072 

0.33 

0.072* 

0.14 

126 

(0.053) 

0.17 

0.175 

0.48 

0.175 

0.35 

127 

(0.H5) 

0.30 

0.39 

0.70 

0.39* 

0.80 

128 

0.37 

0.54 

0.77 

1.0 

0.77 

1.9 

129 

0,90 

0.95 

1.45 

1.5 

1.45 

2.5 

130 

2.0 

1.7 

2.5 

2.2 

2.5 

3-2 

131 

(2.63) 

2.9 

3.2* 

3.2 

3.8* 

3-8 

132 

(4.31) 

4.3 

4.7* 

4.4 

5-0 

4.6 

133 

(6.43) 

6.1 

5.5* 

5.4 

5-27* 

4.9 

134 

(7.80) 

7.3 

6.6* 

6.5 

5-69" 

5.2 

135 

(6.40) 

6.3 

6.0* 

5.9 

5.53* 

5.1 

136 

(6,36) 

6.4 

5-9* 

5.3 

5.06" 

5-3 

137 

(6.05) 

6.0 

6.2 

5.85 

5-24' 

6.4* 

Continued 
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IA3LE  I  {continued) 

Cumulative  Mass-Chain  Yields  of  Fission  Products 


(Values  are 

In  per  ednt 

of  flpelons) 

Mass 

1' usher 

,235 

1 

238 

Pu239 

'Ihemal 

Neutrons 

Fission 

t’eutro.-.L= 

Fission 

Neutrons 

8-Mev 

Neutrons 

Chennai 

Neutrons 

Fission 

Neutrons 

133 

5-7** 

5-7 

6.4 

5.9 

5-5 

5.4 

139 

(6.34) 

6.4 

6.5 

6.0 

5-7* 

5-2 

l4o 

6.44 

6.4 

5-7* 

5-6 

5.68* 

5.o« 

14  i 

(6.30) 

6.3 

5-7 

5-5 

5*2* 

4-7 

142 

(5  85) 

5-9 

5*7 

5.4 

6.69* 

4.9 

1^*3 

(5-87) 

5.8 

5-5 

4.97 

;.4* 

5.0 

ill* 

5.67 

5.1** 

4.9* 

4.3“ 

5-29* 

4.6 

145 

3-95 

4  .2 

3.7 

3-7 

4.24* 

4.4 

1!*6 

3-07 

3.3 

3.1 

3.17 

3-53* 

3-7 

l4? 

2.38 

2.5*# 

2.6“ 

2.7** 

2.92* 

3.0 

l'i8 

I.70 

1.85 

2.0 

2.27 

2.28* 

2.36 

H9 

1.13 

1.3** 

1.45 

1.9“ 

1*75 

■  1.86 

150 

0.67 

0.80 

1.05 

1.45 

1.38* 

1.48 

151 

0,45 

0.50 

0.74 

1.02 

1.08 

1.16 

132 

0.265 

0.31 

0.50 

0.66 

0.83» 

O.92 

153 

0,15 

0.19** 

0.32 

0.41“ 

0.5? 

0.60 

151 

0.077 

0.0-56 

0.19 

0.25 

O.32* 

0.37 

155 

0.033 

0.048 

o.ll 

0.15 

0.20 

O.23 

15' 

0.014 

0.023“ 

0.0  66* 

0.092-*' 

0.12* 

0. 14 

l?  ( 

0.0078 

0.012 

O.034 

O.C57 

0.064 

0.075 

Ir»- 

0.002 

C.9062 

0.016 

0.032 

0.034 

0.043 

131 

0.00107 

0.00}-" 

O.OC90** 

0.017“ 

0.020““  0.025 

Con  timed 
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TABLE  Z  (continued) 

Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  In  per  cent  of  fissions) 


Mass 

u235 

^38 

PuJ 

39 

Number 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission 

Neutrons 

8-Mev 

Neutrons 

Thermal 

Neutrons 

Fission 

Neutrons 

160 

3-5x10-4 

0.0012 

0.0036 

0.0085 

0.0092 

0.011 

l6l 

7.6x10-5 

4.6x10-4** 

9-4x10-4 

0.0044** 

0.0 038*^»* 

O.OOJi 

■•Seymour  Katcoff,  Fission-Product  Yields  From  U,  Th  and  Pu,  Nucleonics  , 

Vol.  16,  No.  4,  p  78-85  U95B) 

r.  Bunney,  E.  M.  Scadden,  J.  0.  Abriam  and  N.  E.  Ballou,  Radiochemical 
Studies  of  the  Fast  Neutron  Fission  of  L'-35  and  U^38t  Second  Inter¬ 
national  Conference  on  the  Peaceful  Uses  of  Atonic  Energy,  A/Conf.  15/F/643* 
USA,  June  1958 

«**G.  P.  Ford,  J.  S.  Gilmore,  et  si,  Plsclon  Yields,  LADC-3083,  1958 

*»«L.  R.  Bunney,  E.  M.  Scadden,  J.  0.  Abrlam,  N.  E.  Ballou,  PI  scion  Yields 
In  Neutron  Fission  of  Pu^^°,  USNTOL-TR-268,  1958,  Unci. 

a.  Parentheses  indicate  estimated  values  or  where  Katcoff 8  value  was  altered 
in  order  to  adjust  the  yields  to  a  gross  sum  of  100  in  each  peak. 

b.  Line  Indicates  division  of  two  peaks  that  woo  used  for  Individual  peak 
sumo. 


Nuclldr 


Zn-72 

2n-7 4 

Ca-72 

Ca-73 

Ca-74 

Ce-75 

0-77 

Ce-78 

Aa-77 
As -78 
Aa-79 

Se.-8ljB 

Sa*-8i 

Se-83 

Br-83 

Br-84 

Kr-833! 

Kr.-85« 

Kr^-85 

Kr-87 

Kr-88 

Rb-88 

Rb-89 

Rb-91 

8r-89 

8r-90 

Sr-9l 

8r-92 

3r-93 

V-90 

Y,-91n 

Yp-91 

Y-92 

Y-93 

Y-9i* 

Zr-95 

2r-97 


•S. 

TABLE  11 


Air  Ionization  Disintegration  ftaltlvlien  tor 
-na_  Fisa  ion  Product  and  oTh'er  RaJlon^n?^^) 


n 

10*9  7/hr-ft2 

v^To/SecT”’ 


ITuclide 


€.75 

4,22 

14.3 

2.11 

12.4 

0.216 

13.1 

2.74 

0.0735 

2.51 

0.0 

0.LL7 

0.0 

12.8 

0.0174 

8.12 

0.0959 

0.973 

0.0316 

7.03 

7.40 

4.30 
12.0 
6.6  7 

0.0 

0.0 

5.03 

6.93 

1.20 

0.008 

3-24 

0.0135 

6.56 

0.817 

6.67 

fc-32 

0.346 


JTb. -95« 

Nbi*95 
Kb, -97b 
Nbp-97 
Nb-98 

Mo-99 

Mo-101 

Mo-102 

Tc-99b 
Tc-101 
Tc -102 

Bu-103 

Pu-105 

Ru-106 

Rh.l03« 

Rh.-105m 

Rhi-105 

Bh-106 

Rh-107 

W-109 

Pd-Ul 

Pd-112 

Ag-109* 

Ag-lii 

Ag-112 

Afi-113 

Ag-115 

Cd  -115 

Cdf-115 

Cd>U7m 

Cd-118 

Cd-120 

In-115 

In-117 

In-ll8 

In-119 

Sn-121 
Sn-123 
Sn-125 
Gn-126 
Sr.- 127 


a 

10'9  r/hr-ft2 
(d is/ sec) 


1.41 
4.56 

4.42 
4.02 

11.7 

0.772 

6.85 

0.0 


0.782 

2.03 

2.45 


2.99 

4.36 

0.0 


0.150 

0.283 

0.559 

1.8l 

5-08 


0.0 

5-00 

0.0225 

0.172 

0.175 

*♦.59 

O.158 

1.28 

O.I87 
1.11 
8.H 
0.0 
3. 08 


1.12 

1.36 

13.5 

0.0371 


0.0 

0.0 

6.19 

0.0 

3.08 
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Table  11  (continued) 


Air  Ior.l2atlon  Disintegration  Multipliers  for 


Ruellde 

the  Fission  Product  and 

Other  R'idlonuclldea 

a 

10*9  r/hr-ft2 

a 

20’9  r/hr-ft2 

rtuellde 

(dlt/see) 

(dls/aec) 

Sb-125 

2.78 

Ce-l4l 

0.470 

Sb-126 

16.1 

Ce.l43 

2-34 

Sb-127 

2.70 

Ce*l44 

0.188 

Sb-i28 

13.0 

Ce-145 

4.22 

Sb-129 

6.10 

Ce-146 

1.60 

Sb-131 

3.66 

Pr-l43 

0.0 

Te.  -125** 

0.0525 

Pr-144 

0.163 

Te,-127 

0.181 

Pr-145 

0.0 

Tei-127 

0.0249 

Pr-146 

6.42 

Tcf -129a 

0.398 

Rd-l47 

O.968 

Te^-129 

1.52 

Te,-131 

8.U2 

Nd-149 

2.31 

Tei-131 

2.08 

IM-151 

6.44 

Te-132 

1.52 

Ra-l47 

Te. -133a 

0.790 

0.0 

Tei-133 

9.12 

P»-l49 

7.25 

Te-13l* 

7.89 

Pb-150 

4.99 

Pb-151 

2.21 

1-131 

2.3S 

Pm- 152 

3-79 

1-132 

13.0 

Pm-153 

4.76 

1-133 

3.59 

0.0544 

1-134 

9.63 

Sou  151 

1-135 

9.60 

Em-153 

0.455 

Sm-155 

1.80 

Xe  -131a 

0.207 

Sm-156 

1.20 

Xe*-133« 

0.1*13 

Sm-158 

3.08 

XeJ-133 

0.304 

0.314 

Xe^-135. 

2.69 

Eu-155 

XOp-135 

1.55 

Bu-156 

5-55 

Xe-138 

7.09 

Eu-157 

3-78 

Eu-158 

7.21 

Ca-137 

0.0 

0.430 

Co-130 

U  .2 

Od-159 

Co-139 

4.78 

Tb-l6l 

0.123 

Ba-137« 

3.65 

Ba-139 

0.888 

tl-237 

1.13 

Ba-l40 

1.10 

U-239 

0.338 

Ba-lUl 

4.22 

U-040 

0.0 

Ba-l42 

5.Y6 

:ip-239 

1.11 

U-lUO 

13.0 

Kp-240 

2.07 

Lo-l^l 

0.400 

Kn-56 

9.40 

La-l42 

11.6 

La-ll*3 

6.67 
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The  gross  Ionization  fractionation  neater,  r^(l),  is  defined  as 
the  ratio  between  the  ionization  rates  at  H  +  1  hour  of  condensed  and  normal 
fission  products.  It  nay  be  expressed  as 

rx(l)  - 


*f^> 


Thus,  all  the  components  for  the  conversion  factor,  K^t),  have  teen  cited 
and  derived.  The  intensity-activity  conversion,  shown  as  Equation  (h2), 
is  very  useful  in  fallout  contamination  studies  and  the  determination  of 
other  properties  of  fallout. 

F.  Biological  Availability  and  Atom-Concentration  Intensity  Ratio 

The  fraction  of  radionuclides  condensed  on  the  outside  of  fallout 
particles  Is  potentially  available  for  biological  uptake.  By  an  analogy 
to  the  Intensity -activity  conversion  presented  in  the  foregoing  section, 
the  ratio  of  the  number  of  atoms  condensed  on  the  exterior  of  the 
particle  to  the  standard  Intensity  at  H  ♦  1  hour  is  given  toy  the  following 
formula 


Na<A) 

rnr 


YAr;(A,  t) 


(**7) 


where 


*&A) 


Is  the  number  of  atoms  (of  the  end  member  of  the  mass 
chain  and  condensed  on  the  outside  of  the  particle) 
which  land  per  square  foot  of  ground. 

The  other  terms  in  Equation  (^7)  are  the  some  as  defined  previously, 
ihe  ratio  Is  colled  the  atom-concentration  intensity  ratio. 

Ac  shown  In  Equation  (12)  the  fireball  rising  time,  t,  is  a  function 
of  the  partlcle-slzo  part  meter,  a.  Because  of  tine  dependence  this  ratio  is 
also  dependent  on  Q.  However,  the  Or  can  be  approximated  and  simplified  by  an 


average  value  o^,  which  is 
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X 

h 


Since  ci^  does  not  vary  with  crocs-wlnd  dis*ftr,ce>  atom-concentration 


intensity  ratio  holds  constant  &s  a  first  approx inatlth*  f°r  ®H  valaco  or"  cross- 

Moreover,  ns  shall 


wind  distance  at  a  given  value  of  downwind  di6tnr.ce  X 
be  seen  in  the  following  section,  a  given  value  of  ci  yepresento  a  group  of 
particles  with  a  snail  range  in  diameters.  Therefore;  fhc  ratio  N^/l  ia  essen¬ 
tially  associated  with  a  given  particle  size. 


The  values  of  Na  (A)/l(t)  for  six  biologically  important  icotopeo  at 


different  downvind  distances  for  a  typical  10  MT  burst  have  been  calculated 
as  shown  in  Table  VI. 

This  ratio  can  be  applied  to  water  contaxnllthtl°n  studies  to  yield 
the  activity  of  soluble  radionuclides  in  water  cupplie8  >  W  tq(^>  O  1° 
replaced  by  rQ(A,  t),  the  ratio  WQfl  will  provide  the  Information  for  the  actlvltyj 
of  insoluble  radionuclides,  since  those  radionuclides  fased  into  the  silicate 
particles  are  assumed  to  be  not  readily  soluble  in  wat®r* 

0.  Ionization  Rate  Contour  Ratios  and  Particle  Size  £><  strlbutlon 

Several  ionization  rate  cootour  ratios  havs  l  ecu  defined  to  deter¬ 
mine  special  properties  of  fallout.  They  aro: 

(a)  The  mass  contour  ratio,  Mr(t),  which  Id  tho  ratio  of  the 
fallout  mass  per  unit  area  at  any  downwind  distance  10  the  ideal  ionization 
rate  at  that  location,  or  mathematically: 


Mr(t) 


where 


Is  the  mass  of  fallout  per  unit  area  ft!  tuV  downwind 
distance,  X 

I^t)  is  the  ideal  ionization-rate  at  3  ft  wLove  an  extended  open 
area  covered  with  fallout. 


(**9) 
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The  mess  of  fallout,  m,  empirically  determined  to  be  a  function 
of  both  a  and  V,  la 


a  -  f(t»)  W°‘9L7  Rg 

(50) 

where 

f(a)  -  5.33  X  10* 12  u’l‘25  mg/f ioslon. 

a 

■  0.1  to  0.9 

-  5.64  x  10* 12  mg/ f loslor.. 

a 

-  0.9  to  2.0 

*  7.14  x  10*13  H^/f lesion 

a 

>2.0 

According  to  Equation  {’1^)1  lb  *8  found  that 
yi)  *  Area  -  D(l)qx  Ua(l'  *  M01  W 


Therefore,  the  mass  cohU™**  ratio  associated  with  H  ♦  1  hour,  Mr(l), 
for  a  ground  surface  burst  will  to# 


Mr(l) 


f(a)  W^’083  _ _ mg/sq  ft 

D^l)  qj^Tr^Cl)  lfpu)  *  l^ljj  r/hr  at  1  hr 


(52) 


(b)  The  fisulon-produOl  contour  ratio,  FPr(t),  which  is  the  ratio 
of  the  number  of  atoms,  or  moles,  of  fission  products  per  unit  area  at  any 
do.wlnd  distance  to  the  Ideal  iohl^Rllon  rate  at  that  location,  or  matho •• 
natlcally; 

PPr(t)  *  *rjhlt)  (53) 

where  R  1<J  the  numbef  >'f  atoms,  or  moles,  of  fission  products 

fp 

per  unit  areft 1 

From  the  dlscuoelon  of  Ruction  6.,  this  contour  ratio,  corrected 


to  B  ♦  1  hour,  for  a  ground  surfsllu  burst  may  be  expressed  as 
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P?r(l) 


1-16  x  10~~*  ru(ro(A)  *  ro{A)1-A 
"  Diir^irQ(i)  i.-jrijvi^ijr 


taoles  fp/sq  ft  /ci  •. 

r/hr  at  I  hr '  (5^ 


(c)  CTe  fraction  of  device  contour  ratio,  FD^ft),  which  Is  the 
ratio  of  the  fraction  of  veapon  device  per  unit  orea  at  any  downwind  distance 
to  the  Ideal  ionization  rate  at  that  location.  This  ratio  when  corrected  to 
H  1  hour  may  be  expressed  as 

.o), 

_ ^1  x  10 

■*X  c*t:u-~'  *i 


™r{l)  “  PdTq,  cv"wi)T;-ri)%  '1(01  r/hr/3« ft  at  1  br- 


(53) 


These  conto-r  ration  predict  various  characteristics  of  local 
fallout,  as  they  yield  information  about  the  mass,  the  number  of  atoms  of 
fission  products  per  ’-tit  area,  the  fraction  of  weapon  yield  per  or.lt  Area 
and  other  properties  of  fallout  at  a  specific  location.  This  information 
is  valuable  for  the  design  of  fallout  shelters  and  other  radiological 
countermeasures.  More-over,  the  use  of  these  contour  ratios  ur.i  their 
scaling  functions  maker  possible  the  extrapolation  of  a  limited  amount 
of  experimental  data  to  cover  a  number  of  operational  cases  since  the  con¬ 
tour  ratios  are  not  constant  but  are  point  functions  whose  values  depend 
oa  many  variables.  Furthermore,  the  average  density  of  the  fission  products 
can  be  estimated  from  the  ratio  of  FPr(l)  to  Mr(l). 

Another  advantage  or  the  C.C.D.  fallout  model  ic  that  particle 
site  groups  for  any  giver,  downwind  location  can  be  estimated  by  a  graphical 
method,  provided  that  the  particle  fall-rate  is  known.  Since  the  fall-rate 
for  a  given  else  parameter  a  does  denote  a  group  of  particles  with  a  range 
of  diameters,  with  the  fell-rate  for  a  given  particle  size  parameter  known, 
the  range  of  diameter;  for  a  group  of  ps:*iclcs  can  be  determined,  depending 
on  the  thickness  of  the  cloud  and  the  altitude  from  which  the  group  falls. 


A  family  of  curves  which  shows  the  size  of  ideal,  spherical  particles  falling 
through  a  standard  atmosphere  from  various  altitudes  as  i  function  of  fall- 
rate  was  prepared  by  D.  E.  CLark^^  and  appears  as  a  sketch  in  Pigure  8* 

When  estimating  the  particle  size  range  at  a  downwind  distance,  it 
should  be  first  decided  whether  the  source  of  the  particle  group  originates 
from  the  stem  or  the  cloud.  The  origin  of  the  particle  group  can  be  generally 
Judged  by  comparing  the  downwind  distance  with  the  downwind  range  of  the 
stem  fallout.  In  case  the  distance  falls  within  this  range,  the  particle 
group  predominant. ly  originates  from  the  stem;  otherwise  it  is  essentially 
derived  from  the  cloud.  Within  the  stem  range,  the  maximum  and  minimum  a 
vnlucs  may  be  calcul-ied  by  Equation  (11),  which  is  r.  laborious  procedure. 

A  computer  program,  cs  shown  in  Figure  7,  has  been  established  to  solve  this 
equation.  The  height  of  fall  for  the  stem  fallout  is  derived  from  Equation  (M) 
which  is 


Z 


{Z0QXZ  *  V> 

°?8  *  V* 


(56) 


For  particle  groups  falling  inside  the  cloud  range  at  nny  down¬ 
wind  location  (X,  V),  the  crossvind  distance  Y  should  first  be  converted 

to  the  cloud  coordinate,  y,  according  to 

* 

y  -  a  (57) 


where  Y  is  the  raxlri'in  half-width  of  the  croscwir.d  distance  on  a 

m 

1  r/hr  contour,  and 
a  is  the  cloud  radius. 

Following  this  coordinate  transfer,  use  equations  (l$ )  and  (l3)  to 
Ktemine  Q  and  the  two  extreme  heights  of  fall,  z  and  ;  .  For  the  actual 
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distance  of  fall,  the  height  of  the  cloud,  h,  should  be  added.  The  fall-rate 
is  evaluated  from  Equation  (10),  which  i6  the  defining  equation  of  a  and 
which  yields 


Therefore,  with  weapon  yield  and  downwind  distance  known,  and  with 
the  help  of  the  fall-rate  curves,  the  approximate  particle  size  distribution 
at  (X,  Y)  can  be  estimated.  This  information  on  particle  size  distribution 
is  useful  for  studies  of  transport  phenomena  in  fallout  investigations. 

H.  Conclusion  -  The  Reality  and  Applications  of  the  O.C.D.  Model 

The  O.C.D.  fallout  model  mainly  uses  the  technique  of  the  scaling 
method,  hence  most  of  its  formulae  are  derived  from  realistic  experimental 
data.  Therefore  it  serves  as  a  very  good  bridge  between  the  ideal  theoretical 
hypothesis  and  practical  empirical  data.  Moreover,  this  model  has  been 
applied  to  several  real  test  shots,  and  was  found  to  yield  very  satisfactory 
results.  On  the  other  hand,  the  fundamental  thermodynamic  theory  on  which 
this  model  is  based  is  quite  simple  and  convincing,  and  all  the  computations 
Involved  are  relatively  easy.  Another  advantage  of  this  model  is  its  vide 
application,  since  it  gives  not  only  information  on  fallout  patterns  tut 
also  that  of  activity,  hence  it  is  very  valuable  to  radiological  research. 
Although  there  arc  still  a  few  defects,  such  as  constant  wind  velocity, 
the  sole  dependence  on  weapon  yield,  and  the  cigar-chape  contour  assump¬ 
tions,  etc,,  this  model  is  generally  considered  very  realistic  and  applicable 
to  the  problems  of  fallout  research.  Its  main  applications  will  be  radio¬ 
activity  prediction,  countermeasure  design  and  biological  protection. 


Ja¬ 

in.  COMPUTATIONAL  METHODS  OF  THE  FALLOW  MODEL 

A.  location  Rate  Contour* 

Although  the  physical  characteristics  and  the  mathematical  approxi¬ 
mation  for  the  rising  stem  and  the  cloud  have  been  discussed  In  detail, 
little  has  been  said  about  how  the  fallout  contours  are  determined.  Even 
though  the  available  data  from  actual  land  surface  shots  Is  small  ar.d  limited 
primarily  to  low  weapon  yields,  some  general  conclusions  have  been  drawn 
from  the  teat  data. 

Because  most  of  the  observed  properties  of  the  fallout  patterns 
are  In  terms  of  Intensity,  the  fallout  contours  are  usually  given  In  units 
of  roectgens/hour  (where  the  time  was  corrected  to  a  standard  time  of 
H  ♦  1  hour).  The  general  shape  of  the  fallout  pattern  can  be  approximated 
by  the  overlapping  of  ellipses,  for  both  the  stem  and  cloud  fallout. 

The  dimensions  necessary  for  constructing  the  ellipses  are  the 
distances  from  ground  zero  along  the  downwind  axis  where  the  lateral  dis¬ 
placement  la  zero,  (Y  «  0).  First,  to  construct  the  intensity  vs.  distance 
curve  (see  Figure  **)  the  downwind  distances  of  interest  are  the  highest 
intensity  and  a  preselected  low  intensity,  which  varies  with  wind  velocity, 
but  for  a  15  mph  wind  la  1  r/hr  at  1  hr. 

For  the  stem  these  locations  have  been  labelled  as  follows: 

X^  and  are  the  upwind  and  downwind  distances,  respectively,  from 

ground  zero  to  the  1  r/hr  at  1  hr  Intensity. 

Xg  and  are  the  distances  to  the  upwind  and  downwind  shoulders, 

respectively,  of  the  intensity  ridge. 

For  the  cloud  the  following  notation  is  used; 

Xj  and  Xj  are  the  arrival  location  for  the  group  of  particles  of  the 
the  selected  low  intensity  from  the  rear  and  front  of  the 
cloud,  respectively. 


DOWNMNO  DISTANCE  X  (MU.ES) 


X,  is  the  distance  from  ground  zero  to  an  intermediate  point 

o 

between  Xj  and  X^. 

X^  la  the  downwind.  distance  from  ground  zero  to  the  peak  of 

the  Intensity  contour. 

Xg  is  the  downwind  distance  from  ground  zero  where  the  contour 

Intervals  has  a  maximum  half-width. 

The  values  for  these  quantities  are  determined  by  equations 

appearing  In  the  Ionization  or  Intensity  rate  computer  program.  Eecause 

of  the  limited  characters  available  on  the  computer  printer,  the  notation 

used  In  the  analysis  of  the  model  had  to  be  changed.  The  corresponding 

notation  between  the  model  and  the  computer  program  is  given  in  Table  III. 

The  following  discussion  describes  the  method  for  calculating  the 

X * 3  and  how  the  scaling  functions  can  be  found  in  the  computer  program  in 

Figure  5.  First,  yield-dependent  scaling  functions  for  the  size  parameters 

a  ■  a,  y  o^,  a^,  ,  .  .  ,  have  been  determined  from  data  that  has  an 

effective,  or  average  wind  epeed  of  15  mph.  Thece  ocaling  functions  are 
VV 

then  multiplied  by  to  correct  for  the  wind  velocity  V^.  I2  ytho 

ionization  rate  between  Xg  and  X^, is  calculated  by  two  formulas.  For  yields 

W  *  100KT,  I2  3  -  Kg>3U)  VM°  Where  fche  values  for  K2,3(l)  *^2,3 

and  n  were  generated  by  fitting  curves  through  the  values  of  a  Table  in 

Chapter  Three  of  Carl  F.  Miller's  manuscript  "Fallout  and  Radiological 

Countermeasures" .  For  W  <  100  a  curve  was  fitted  through  the  first  three 

values  Of  I,  ,  given  in  another  table  from  Miller  and  adjusted  for  wind 
**  >  j  v 

velocity  by  multiplying  by  The  resulting  equations  appear  on 

lines  07^1  through  08C2  in  the  computer  program  appearing  in  Figure' 5. 


In  all  the  equations  that  follow  the  scaling  function  formulas 
and  constants  are  given  in  Table  IV.  For  any  scaling  function  dependent 
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figure  5  (Continued).  Computer  Prog  ran  for  Ioniration  Rate  Contour* 
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Figure  5  (Continued).,  Computer  Program  for  Zonitation  Rate  Contour. 
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Figure  5  (Continued).  Coraputer  Progma  for  Ionltatlon  Rate  Contour# 
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figure  5  (Continued).  Computer  Program  for  I o filiation  Itete  Contour* 
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Figure  5  (Continued).  Computer  Program  for  Ioniiation  Rate  Contour* 
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Figure  5  (Continued),  Computer  PrograHl  tor,  Xonlx&tlon  Rate  Contour# 
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ALP  (  ) 

V  “5*  °fe»  V  <ty  a9 

F-0  If  V  s  28  KT 
(used  to  differentiate 

Af4^3 

%3 

scaling  functions) 

ALP5 

°5 

G 

G«1  If  V  <  9  KP 

G=0  If  W  i  9  Kf 

AM^P 

"5 

a 

h 

A9 

«_ 

i  (  ) 

6 

■ 

b 

Il»  h>  x5»  r6»  I7'  T8*  *9 

ic  (  ) 

Intensities  of  chosen 

MfU.  DYS 

* 

(0,  I7),  (D,  I) 

contour's 

MY'',  by3 

* 

(0,  L  ),  (Y  •  1) 

u 

Intensity  of  chosen 

*=*0  • 

contour  I ?  , 

H12„  M2 

« 

(*1»  Il^>  *2,3) 

I0R3 

K2  (1) 

HJ't.  B34 

* 

^3*  *2,3^’ 

123 

J  *  CIJ 

*2  3 

H?(jt  B56 

* 

<X5,  Ic),  (V  I6) 

J 

Dummy  variable  used  as 

B67 

* 

(X6.  I6),  (*J,  Ij) 

"1"  in  I  ,  etc. 

H  B  789* 

(*J,  Ij),  (Xy  I9) 

XA 

Ka 

OA 

CA(I) 

KM  (  ) 

Kid)Ja,  Kl(l)Ha 

*c 

K 12 

V2 

01 

U7 

L(  ) 

Array  containing  constants 
for  scaling  functions 

B 

D 

LA 

In  a 

Df-P 

O’ 

LARS 

In  a/R 

Number  of  contours 

0 

studied 

LA3 

In  & 

0 

Table  III  (Continued) 
LIST  C.  SYMBOLS 


Rotation  of 
co^iter 
Tr-.o-raq 

LAO 

LI  {  ) 

“7 

IM 

M  {  J 


PHI 

<K 

RS 

SA 

8VW 

V 

vw 

VUB 

WSI 

W 

X(  ) 

XCl 

XC3 

xcu 


notation  due 
primarily  to 
_ Miller 

la  a 

o 

In  1° 

1 

In  V 

Array  containing 
constants  for 
scaling  function* 

Uccd  in  equations 
involving  X* 

V  *6'  f7 

°^,3  or  ai  *  5 

R. 

SA(I) 

8(^v) 

V15 


Rotation  of 
computer 
program 

XC6 

XJPRC 

XPRS 

XX 

xxc 

xxrr 

X4U 

X7 

X8 

X07 

X9 

X98 

yicc 


Rotation  duo 
primarily  to 
Miller 

7,  9  ~ 

X* 

X* 

X 

X 

lajfA 


In  (Vr9) 

X982  .  X872 
X?  in  feet 
Xq.4h  feet 

*8  '  *7 
Xj  in  feet 

X9  *  Xfl  ln 

Maximum  cloud  width  at 
chosen  contour 


a 

YIC3 

Maximum 

chosen  c 

T 

V 

YS 

Y 

V 

s 

YY 

Y 

Xl’  X2’  X3»  \>  x5* 

YOS 

Y° 

V  X7'  X8’  X9 

Y8 

l,  2  •* 

Y015 

V«) 

2. 

3,  «*  •* 

23 

5,  7  ** 

s 

20 

2 

0 
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Table  III  (continued) 


LIST  OP  SYMBOLS 
Continued 


Notation  of 
computer 
program 


Notation  due 
primarily  to 
Miller 


Notation  of 
computer 
program 


Z99WDlDVELOCrmiPH  vy  299YIEIDKX 


Notation  due 
primarily-  to 
Miller 


V 


*  The  slope  and  intercept,  respectively,  for  a  semi -logarithmic  line 
through  the  two  points  whose  coordinates  are  given. 

**  The  downwind  distance  (for  the  chosen  intensity)  to  the  fallout  pattern 
between  the  points  indicated. 
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table  iv 

Data  Table  for  Ionization  Rate  Contour  Program 


for 

rt 

Hui 

m 

(a/b)2 

i 

0.486 

0.262 

In  a/r 

a 

2 

1.070 

0.098 

In  A 

3 

3.389 

0.431 

R 

s 

4 

2.319 

0.333 

H 

5 

3.820 

0.445 

6 

4.226 

0.164 

°^3 

7 

.O.509 

0.076 

In  a 

0 

6 

2.880 

0.348 

X. 

9 

3. 308 

0.496 

i 

10 

3-564 

0.319 

°5 

11 

-0.054 

0.095 

X, 

12 

3.850 

0.481 

6 

13 

4.255 

0.200 

X. 

14 

3.862 

0.586 

7 

15 

4.268 

C.305 

x„ 

16 

4.005 

0.596 

8 

17 

4.410 

0.315 

x„ 

18 

5.190 

0.319 

9 

19 

5.202 

0.3H 

< 

20 

3.223 

0.400 

Tor 

V 

A* 

Slim  a  1 

a 

23 

1.371 

-0.124 

9 

24 

O.96O 

0.146 

*4 

25 

0.270 

0.089 

26 

-O.176 

0.022 

a6 

27 

0.030 

0.036 

°7 

28 

0.043 

0.141 

“8 

29 

O.185 

0.151 

kJT 

30 

-3-286 

-0.298 

ya 

31 

-2.889 

-0.572 

vlX 

32 

-3.185 

-0.406 

5*a 

33 

K,X 

34 

-1.134 

-0.074 

6  a 

35 

-1.225 

-0.022 

K.JI 

36 

-O.989 

-0.037 

ra 

37 

-1.079 

-0.020 

Vc 

38 

-2c 166 

-0.552 

i/H 

39 

-0.431 

-0.014 

40 

-0.837 

0.267 

41 

1 

fO 

vn 

O 

-0.404 

*12 

42 

-2,600 

-.337 

v  21  3-644  0.467 

5  22  4.049  0.186 
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variable,  for  example  H  (height  of  cloud),  the  program  determines  which 
formula  should  be  used,  The  equation  for  H  appears  on  line  06l3  of  the 
contour  program. 

The  following  chain  is  then  started 

3  -  r?  (5  +  F,  V)  line  0613 

FF(N,  W)  -  EXP(lfT(X,  V))  line  0454 


UFO?,  W)  -  2.3055851  L(h)  *  HI  *  M(tl)  line  0U67 


The  program  then  finds  for  W*  $  28  KT  that  F  *  0  since 
V  greater  28  is  false  and  F  Is  a  Boolean  variable.  Therefore  for  W  •  5, 
one  finds  in  the  Table  IV  that  L(5)  •  3*820  nnd  M( 5)  «  0.445.  The 
complete  equation  of  H  for  W  $  28  KT  la 

H  •  exp((2.302585l)(3*82O)  ♦  log^l#  •  0.445)  , 

log^H  -  (2.3025851)(3*^0)  ♦  log^w  •  O.W»5  , 

or 

log^H  •  3*820  ♦  0.445  logyjW. 


The  variation  of  and  X^  with  wind  speed  depends  on  the  wind  modified 
parameter  or  (lines  0713»  0717)*  The  equations  for  X.  end  X,  are  then 


*2  • 


2*303  logor2 


*1,2 


VWtf  ’  86;36) 

(1  ♦  9.273  X  10*\vy) 


(59) 
(59a) 

(59b) 

(60) 


(61) 
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or 


a  (a  2  -  1900) 
«•  *+  o 

~~a  *  O'.  02  0‘ 


(6la) 


for  1=  npH  *lnd  velocity . 

<lh«  equation-  for  the  distances  X^  through  X?  (lir.e  0995)  are 


given  by 


The 


For  i  -  6,  7,  8  and  9 

xt  -  6.C0  x  loV3’**1*5  at,  tf  -  1  to  28  KT 

X:  -  1.68  X  10Uy°*l6U  at,  «  -  28  to  105  KT 

atove  equations  also  hold  for  ^  when  0^  le  greater  than  a/h.  When 


q5  Is  lc-ee  thai  or  equhl  to  a/h,  then  X§  lo  determined  by 


(62) 

(62a) 


X5  -  6»C-d  X  loV*^  -1.1*0  x  loV'*0  <n.OC#-**  +  *?  »  ^  -  1  to  28  KT  (63) 


or 


lit/,  N  loV-^oi  ,  I. to  x  loV*300  /  3.06W0-262  +  ;  V-28  to  105  KT  (63a) 

?  5 

vw 

jiL*t  the  Integrities  1^  I4,  and  are  set  equal  to  33  mpS  • 

The  eqUA*.  i0»->  for  I;,  «d  1?  are  on  lines  0951  through  0961.  Iq  Is 
determined  fn**  a  sesi-lornrl  thole  line  between  (X?,  I?)  and  (X?,  1^ 
with  the  V*4'»e  o:  Xg  as  previously  calculated. 

Al’U  r  x^  through  X^  and  their  corresponding  intensities  have  been 
,,  cun  then  tc  drawn  of  the  Intensity  profiles  or  elevations. 

x  OUDG  ^  il  F  •  * 

,  „  rh,  1(0  u^ir.-.t  distance  X  for  a  10  CT  weapon  is  given  in 
A  plot  or  h-'t'io'  V  ~ 

Figure  U.  t'or  any  set  o:'  intensities  the  intersection  (downwind  distance  X) 
of  the  concur  cUlpse*  v-‘h  V  -  0  can  be  read  off  the  graph. 


•59' 


In  the  computer  program  semtlog  equations  were  found  for  the 
Intensity  vs.  Distance  curve  so  that  any  pre-selected  intensity  contour 
intervals  were  Included  la  the  computer  printout. 

Next,  to  find  the  lateral  locations  or  vidth  of  the  contour 
patterns  for  the  stem,  the  half  width  Yg  for  the  stem  fallout  is  needed. 

Yg  is  defined  as  the  lateral  distance  from  the  center-lir.e  of  the  sten 
pattern  to  the  lr/hr  at  1  hr  (for  15  mph  wind)  contour  (lines  1216-1220). 

A  seal logarithmic  line  from  (l2  y  Y  ■  0)  to  (1,  Yg)  is  then  calculated 
on  line  1315-  Using  this  equation  the  lateral  displacement  of  any  intensity 
contour  of  the  stem  can  he  found. 

The  maximum  pattern  half  width  for  the  cloud,  Yg,  Is  the  maximum 
distance  from  Y  -  0  to  the  lr/hr  at  1  hr  contour  when  the  wind  velocity 
is  15  mph.  The  calculation  of  Yg  with  its  variation  due  to  wind  speed 
occurs  on  lines  105*»  to  1092. 

To  find  the  lateral  displacements  of  the  contour  locations 
between  (X^,  0)  and  (Xg,  Yg)  a  serailogarlthic  line  from  (ij,  Y  •  0)  to 
(1,  D  -  7  Yg2  ♦  (Xg  -  X?)2)  is  calculated.  In  lines  1236  through  13**7 
the  X  ♦  Y  and  foci  are  then  found  for  the  ellipses.  Using  the  information 
lu  the  computer  printout  on  ionization  rate  contour  map  can  he  drawn. 

(figure  6). 

B.  Fractionation  Humbert 
1.  First  period  of  condensation 

For  a  ground  surface  burst,  the  ratio  of  the  total  moles  of 
liquid  carrier  to  the  molar  volume  n(x)/V  is  computed  by  Equations  (32)  and 
(33)  with  y  -  10  MT  or  10,000  XT,  and  T  -  l673°K.  They  yield 
t  ■  53  see 

(n(0/V]KT  -  3-25  x  10*2  atm. 


\0 


For  &  typical  element  Sr-89,  the  procedure  of  finding  its  fractionation 
number  of  flrat  period  of  condensation  rQ(89)  is  given  below: 

Tecay  chain  of  Sr-89: 

88.2 < Xr — Kb  — w  Sr 

Se-*-3r  - ' 


Here  a  neutron-emitter  Is  involved.  I  at  the  Police  and  Bnllou 
date  have  already  included  this  effect,  hence  no  alteration  needs  to  be  put 
on  the  number  of  atoms  N(A,  t)  in  their  work.  Moreover,  the  valuca  calculated 
by  them  according  to  Clecdenin  theory  are  generally  considered  to  be  more 
reliable  it.d  also  used  in  the  following  calculations. 


Se 

Br 

Kr 

Kt 

Sr 

V’o  k  101 1,106 

1.04 

0.603 

4.6l  x  10*13 

kj  ■  Uuij/vTKf  3,108  *  i0 

32.00 

18.62 

1.4i8  x  10*11 

N>,  t) 

0.21 

309 

94 

3 

C<'A,  t) 

0.21 

309.21 

403.21 

406.21 

MtJil 

1  +  *J 

6.363  x  10*  3 

0 

4.791 

3 

r  "LbJl 
u  1  *  V 

1  Ko 

6.363  x  10" 3 

(2)  6363 

4.797 

7-797 

r  (39) 

0 

0.0)03 

0 

0.0H9 

0.01919 

Vito:  data  un.i<;tcd,  usually  r.'ellgible 


flic  r  \ O'?)  for  ir,  Kr,  Pb  un  the  by-product.,, 
o 

2,  Stcoil  period  of  coadvnrntlou 

The  frvctior.stlon  number  of  second  period  of  condencatlon  '.A,  t) 
>p:nds  very  ra.tch  or.  the  time  and  the  e  It  rents  involved  in  the  decay  chain. 
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When  the  time  is  the  sublimation  pressures  of  moot 

P.V 


elements  are  ctill  quite  significant,  ar.d  the  values  of  are  usually 
greater  than  those  of  0-23  EW®.  In  the  second  period  of  condensation 

V 

4  portion  of  each  such  element  remains  In  the  vapor  phase.  Therefore, 
In  computing,  the  ratio  of  these  two  values  should  be  set  equal  to  unity. 
As  time  increases,  the  sublimation  pressures  of  all  elenenta  except  the 
rare  gases  and  As  and  Se,  decrease  sharply  and  eooo  become  negligible. 
Therefore,  the  second  period  fractionation  number  for  a  mass  chain  not 
Involving  the  rare  gases  and  Aa  and  De  will  be  greatly  simplified} 


r;  (A,  t) 


1  -  rQ(A»  t). 


«*) 


Four  typical  examples  are  given  belovi 

(a)  10  KT  ground  surface  burst,  60  sec  after  burst,  Kr  (83)' 

Sec ay  chain j 

.Se, 


Ce-*>As 


Br-*»Kr 


This  Involves  a  split  chain,  but  the  Bolles  and  Ballou  data  also  have 
taken  this  effect  into  account. 


Ce 

As 

*2 

Br 

Xr 

r0(83) 

— 

0 

0 

0 

0.0066 

0.0066 

l*ro 

— 

1 

1 

1 

0.993^ 

0.993k 

W(83-6o) 

— 

0,13 

21.80 

16.80 

12.30 

0 

— 

0.13 

21.93 

16.93 

51.03 

51.03 

K°.(l-ro)H 

— 

0.13 

21.80 

16.80 

12.20 

0 
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n°  - 

0.13 

21.93 

16.93 

50.93 

50.93 

Pj  (2)8216 

(9)3578 

34.53 

34.53 

0.2086 

* 10  7Al- 

7  AX 

7.41 

7.4i 

7.41 

7.41 

ZN"  _ 

m 

1 

1 

1 

0.9980 

0.9980 

r;(83)  - 

0 

0 

0 

0* 

0 

*  Rotes  r^(A,  t)  can  not  be  a  negative  value. 

henoe  the 

smallest  Is 

zero. 

(b)  10  MT  ground  surface  burst,  1.47  nln  after  burtt,Sr('i)(90). 

Decay  chains 


er - — Kr- 

Br 

Kr 

—sr - 

Rb 

Sr 

V 

ro(9°) 

— 

0 

0.036 

0.1378 

0.1378 

X*ro 

— 

1 

0.964 

0.8622 

0.8622 

H( 90,88) 

— 

67 

338 

106 

0 

a» 

— 

67 

405 

511 

511 

iT-d-r^l! 

— 

67 

325-8 

91-4 

0 

in" 

— 

67 

392-8 

392.8 

392.6 

(4)3050 

00 

(6)4338 

0 

0 

y  iA0 

0.23BWKf  X  10 

1.92 

1.92 

1.92 

1.92 

1.92 

ZN" 

231 

— 

1 

0.9698 

0.7686 

0.7686 

r;(90) 

— 

0 

0 

O.0936 

0.0936 

(e.)  10  MT  ground  surface  bu 


ro(T?) 

1-r 

o 

:<(  ?5,68) 
,N°-(l.ro)N 

*  108 

v  * 

J 

Vp  K° 

w 

ro<75) 


rst,  1.47  min  after  burst, Ge(75), 
t^cay  chair,. 

C  u  ■"  m  in  *•  Ga-  ~  ■■  Ce 


Cu 

Zn 

Ga 

Cc 

— 

l 

1 

0.6716 

0 

0 

0.3284 

— 

(3)12 

0.0459 

0.0536 

— 

(3)12 

0.0460 

0.0996 

— 

0 

0 

0.0176 

0 

0 

(11)2208 

i.9e 

1-* 

1-92 

1.92 

(3)4263 

— 

0 

0 

0.0018 

— 

0 

0 

0.3266 

Koto:  In  computing  Yy  the  individual  mass  number  dependent  conversion 
iuctor  QVA)  for  U235  to  u23®  should  be  multlpUed  to  each  radionuclide. 

“J  "  O*  a(A)  can  be  calculated  from  Table  I 


(d‘)  10  W  ground  6Jrfae*  6.77  »in  after  buret,  Sr(89). Decay  chain 

Wt-s  fhewn  in  Section  (a.) 


Br 

Kr 

Kb 

Sr 

0.0303 

0 

0.0119 

0.0192 

O.9697 

1 

0.9881 

0.9808 

■‘.^4.^1  fx.'f i,'  1.', 
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N(89,4o6) 

— 

=** 

84 

267 

36 

a 

— 

— 

84 

351 

ko  7 

U"«U- r0)X 

— 

- 

84 

O* 

84 

84 

6a 

os 

0 

0 

cr  ci 

1 

0.2393 

O.'OO* 

r  (89) 

— 

=. 

0 

0.71*86 

0.7 6-04 

3.  Cross  Fractionation  Number 

The  gross  fractional  loti  number,  r^,  or  rQ,  is  defined  by  Ejuaticn  C'?). 
The  tvo  Ionization  rates,  lf  (0  **nrt  lfp  ln  tht3  formula  are  defined  ly 
Equations  (45)  and  (46)  rocpectlVe^y*  The  computations  are  straightforward, 
though  tedious,  once  the  values  of  the  fractionation  numbers,  disintegration 
multipliers  and  activities  an?  fclitwn.  The  variations  in  the  ionization  rates 
of  normal  fission  products  with  tlne>  or  c Imply  the  decay  of  normal  fls- 
products,  for  U-235,  y.238  and  Fus239  hive  been  calculated  and  are  tab.  _  ■ 

in  Table  V,  From  this  tabic,  tit*.*  normal  ionization  rate  of  1J-238  at  H  ♦  1 
hour  can  be  Interpolated  and  ol/lniftad  as 

ifp{i)  •  6.973  x  10' 11  (r/hr  at  1  hr)  /  (fission /cq.  ft.)  (65) 

For  a  10  KT  surface  burst,  the  condensed  ionization  rates  for  four  different 
fireball  rising  times  have  been  *‘t»»'p-*tcd  and  are  listed  with  the  values  of 
gross  fractionation  numbers  comp  j*  *  d  i  rom  them: 

t  1  *{l)  r  (t' 

’M*  X'  ' 

2,'|08  x  10‘13  0.388 


53  sec 
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60  sec 

3.255  x  10‘13 

0.467 

1.47  nla 

3.873  x  10‘13 

0.554 

2.15  ms 

5.0P9  x  10'13 

0.721 

{since  the  fireball  rlrirr 

time,  t,  1 3  related  to 

the  particle 

meter.  or,  ty  Equation  *2x-r,  r^(t)  is  also  a  function  of  0(,  eooetlmea 
denoted  as  r„(t). 

•j'i  Act  1  v  ity  •.Concent  rvtf  an  tr.  tensity  Ratio 

Vlth  Equalling.  tiA  (48),  the  activity  concentration  Intensity 
j u.t lo  for  any  Isotope  at  v,y  downwind  location  can  be  readily  computed.  The 
values  oaf  the  ratio  for  ’Ax  biologically  Important  isotopes  in  a  10  MT 
surface  burst  are  c  tle-JUt.  &rA  shown  in  Table  VI, 
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TABLE  7 


Decay  of  floras!  Fission  Producte  from  U-235.,  U-233  and  P*->»239 
(r/hr  at  3  ft  above  an  Infinite  plnnc  Tor  10*  fissions  r*~r  sqT  ft.) 


Age _  _ U.235  U-2  38  _ Pu-239 


Year* 

Day* 

Hours 

Thermal 

Fission 

(8  Mev) 

Thermal  ] 

Fission 

0.763 

(8) 

(*) 
9977v  ' 

(8)9960 

(8) 

9329 

(8)8907 

(8)8750 

1.12 

(8) 

6648 

(8) 

6632 

(8 

6172 

(8)5966 

(8)5761 

1.64 

(8) 

4149 

(8) 

4153 

(8) 

3827 

0)3592 

(8)3537 

2.40 

(8 

2453 

(8 

2484 

(8 

2256 

(3)2071 

(8)2065 

3.52 

(8 

1410 

(8) 

1450 

(8) 

1303 

(3)1166 

(8)1196 

5.16 

(9)8079 

(9) 

6418 

(9) 

75 82 
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>1039  ( 
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(11) 

16774  ( 
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16899  ( 

11)7244  (11)6655 

i4.4 

(11) 

14490  ( 

:  11)4587  < 

;n 

>4462  ( 

11)4650  ( 

11)4226 

21.1 

(11 

12910  ( 

In; 

>2940  ( 

[nj 

>2845  (11)2953  ( 

11)2660 

30.9 

( 

;n; 

I1813  (11)1807  < 

:u: 

>1762  (11)1837  ( 

11)1645 

45-3 

1 

11 

(1061  (11)1039  ( 

I11 

>1034  (11)1092  ( 

12)9777 

66.4 

12 

16055  (12)5807  ( 

12 

>5910  (12)6360  (12)5728 

97.3 

:i2 

»3676  (L2' 

)3497  ( 

[12 

>3559  (12)3896  ( 

12  )  3543 

143 

;i2 

>2170  1 

12)2090  < 

12 

>2079 

12  2320 

12)2200 

208 

.12! 

>1194  1 

[12, 

)ll64  | 

k! 

>1133 

[12)1287  ( 

12)ll80 

301 

[13; 

>4874  l 

[13)4790  1 

[13[ 

>4733 

13)5707  ( 

13)5170 

1.2 

438 

[13, 

)1399  1 

13)1373  1 

13, 

1525 

13)2135  ( 

,13)1864 

•  78 

650 

14 

3884  l 

14)3758  1 

14 

5517 

14)9083  ( 

,14)7690 

2.60 

14 

2031  1 

14 

)1975  1 

14 

3160 

14)4964  f 

,14)4352 

3.00 

14 

1444  1 

(14)1432  1 

14 

2213 

14)2692 

,14)2594 

5-58 

(l4 

>1154  i 

(14 

)1158  1 

(14, 

)lC03 

[l4)l442  ( 

[l4)l6u 

8.18 

(14 

)1026 

(14 

)1021  1 

(14 

)1291 

(15)9971  (14)1225 

12.0 

15 

9432 

(15 

)9293  1 

(l4 

)l094 

15  3452  ( 

[14)1057 

17.6 

15 

)8310 

(15 

)8211  i 

(15 

)9l64 

(15)7377  ( 

,15)9160 

25.7 

(15 

)7l83 

(15 

)6987 

(15 

)7431  < 

(15)6219  1 

[15)7668 

Numbers  In  parentheses  indicate  the  number  of  zcroo  between  t..e  decimal  point 
and  the  first  significant  figure. 
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D*  Particle  Site  Distribution 

ttiere  la  almost  no  data  available  on  particle  size  distribution, 

nor  has  any  rigorous  theoretical  approach  been  made.  However,  using  the  scaling 

functions  of  the  previously  discussed  simplified  fallout  model,  estimates 

cm  be  made  of  the  particle  size  groups  that  fall  at  any  downwind  location. 

Given  the  fall  velocity  and  the  initial  altitude,  the  size  of 

spherical  particles  falling  through  a  standard  atmosphere  can  be  found 

in  curves  prepared  by  D.  E.  Clark. At  any  downwind  distance  X  the 

first  information  needed  is  the  maximum  and  mimimum  a  size  for  the  particle 

groups  arriving  at  that  location.  Because  the  formulae  for  a  and  a  . 

max  min 

depend  on  whether  the  a  group  descends  from  the  stem  or  the  cloud,  it  was 

decided  that  the  X -coordinate  intersection  of  the  logarithmic  line  connecting 

(X^,  Ij)  to  (XU.  1^)  (point  3  to  point  U  in  Figure  b)  with  the  logarithmic 

line  from  (X^,  1^)  to  (Xg,  Ig)  -  extended  if  necessary  -  would  be  the 

division  mark.  For  any  downwind  distance  lose  than  this  intersection  X 

coordinate,  the  source  of  the  particle  group  ie  assumed  to  be  the  stem. 

Then  and  for  the  stem  may  be  solved  by  graphical  methods  or 

by  the  following  computer  approach.  Essentially  a  Hevton-Raphsoa  Iterative 

procedure  is  performed  on  the  following  function  to  find  separately 

a  and  or  . 
max  mi  a 


f(a) 


».(a  •  aft 

2.303(0  *  ~) 
*Z 


a(zoa  -  V) 

log  (X  -  - jy-r - 1  *  log  a„ 


a  ♦ 


^T 


(66) 


where 


k&  is  the  exponential  constant  for  the  fireball  major  axi3  exj/tnslon, 
ZQ  is  a  yield  dependent  multiplier. 


t 
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kj  Is  ao  empirical  Inverse  time  constant,  0.011  sec*\ 

Vy  Is  the  viad  velocity,  usually  assigned  to  be  10  mph, 

P  Is  the  constant  In  the  empirical,  formula  of  Y^/V^,  0.95  whea 
Z  •  5,C00  to  50,000  ft.,  d  ■  200  to  1,200  microns, 

q  Is  the  constant  In  the  empirical  formula  of  V^/v^.,  1-02  *  10*^ 
when  2  •  5,CC0  to  50,000  ft.,  d  -  200  to  1,200  microns,  and 
®0  is  the  uajor  cesl-oxis  of  the  fireball  at  ground  zero. 

The  iterative  procedure  le 


0L  ■  Initial  estimate 

O 

(67) 

“W  *  S  ' 

(68) 

If  the  sequence  of  a’ a  approach  a  limit,  then  this  at  is  such  that  f(ot)  i  0} 

or  a  la  approximately  equal  to  o  or  a  .  depending  on  the  sign  of  the 

iwu  nu 

number  enclosed  in  the  absolute  value  symbol  of  Equation  (56). 

When  the  downwind  distance  )»  greater  than  or  cquil  to  the  X 
intersection  coordinate,  it  is  assumed  that  the  fallout  originates  from 
the  cloud.  Por  particle  groups  falling  within  the  cloud  area  at  any  (X,  Y) 
location,  the  maximum  and  minimum  a  values  can  be  calculated  byi 

a  .  — -  ♦  <<?  -  va)  h2  .  b2(i  ■  v3/»3) 

*  h2  -  b2  (1  .  y2/.2) 

where  Y  is  translated  to  y  by  Equation  (57). 

vv 

For  any  Ot  value,  the  fail.rate  is  given  by  •  JJ"  •  The  height  of  fall 
for  the  stem  fallout  ia 

<y  h  •  v> 


Z  m 


(70) 


For  the  parr.lcle  groups  that  originate  from  the  cloud  the  maximum  and 
minimus  height  of  fall  is  given  by 


Using  this  information  calculated  in  the  particle  size  parameter  program. 
Figure  7,  and  a  set  of  curves  prepared  by  D.  E.  Clark,  see  Figure  8,  the 
approximate  particle  size  distribution  at  (X,  jf)  can  be  estimated. 

For  a  10  MT  detonation,  and  for  various  downwind  distances  the 

maximum  and  nini=r_n  a;  the  falling  velocity,  Vf:  the  initial  altitude  Z 

and  particle  diaojter  u  are  given  in  lab  Lb  VII. 

(*»y 

Recent  vork  by  Clirk  and  Cobbin'  '  is  closely  related  to  the 
subject  of  this  report.  These  authors  used  the  same  fallout  model  aa 
described  herein  for  the  calculation  of  particle  size  and  radiation 
intensity  at  varying  distances  from  ground  zero.  The  results  obtained 
are  in  good  agreement  with  those  herewith  reported  for  weapons  of  similar 
yield. 
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Figure  ?•  Computer  Program  fnr  Parsunetera 
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Table  VII 

particle  Site  Distribution  at  Various  Downwind  Distances  for  5  WTT  Veapon  Yield 
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IV.  IKPXED  RADIOACTIVITY  EC  SOI  IB 

Elements  can  be  made  radioactive  by  the  action  of  thermal  neutrons 

which  excite  the  stable  atoa  making  it  unstable  and  radioactive.  Ac  thermal 

neutrons  are  liberated  in  a  nuclear  explosion,  there  will  be  some  induced 

radioactivity  In  the  soil  around  the  detonation  site.  A  small  portion  of 

this  coil  will  be  taken  up  into  the  fireball  of  the  bomb  (for  a  ground  burst 

or  near  ground  burst)  and  become  part  of  the  total  fallout. 

A.  V.  KlsmonV^  in  his  report  concerning  the  potential  radionuclides 

produced  In  weapons  detonations,  lists  the  activities  of  the  various  nuclides 

26 

that  are  induced  In  the  soil.  A  total  of  10  neutrons  rcr  KT  was  assumed  to 

27 

be  liberated,  which  agrees  within  a  factor  of  15  with  that  of  1.5  x  10  1 

(l4) 

neutrons  per  KT  calculated  from  Classtone  .  The  composition  of  the  coll 
V0.1t  van  used  for  those  calculations  Is  given  in  Table  VIII. 

Figure  9  gives  the  decay  of  the  Induced  radioactivity  In  coil.  The 
data  van  obtained  by  taklog  the  contributions  of  the  various  nuclides  given 
in  Table  IX  and  sunning  then  for  the  various  times  after  detonation  up  to  one 
year. 

Figure  10  ,  "Integrated  Neutron  Flux  as  a  Function  of  the  Slant 

Range  In  Air  of  0.9  Sca.Devel  Density  for  a  1-Klloton  Explosion"  chows  tho 

neutron  flux  verouo  distance  for  a  1  KT  bomb.  These  curves,  as  given  by 

Classtone were  Integrated  over  ar  area  vlth  a  2,500  yard  radius  ond 

were  used  to  correlate  Klement'o  values. 

Senftle  and  Charpioo^^  give  a  detailed  discussion  concerning 

Induced  radioactivity.  A  method  for  calculating  the  Induced  radioactivity 

Involving  tables  and  formulas  is  alto  given.  This  method  was  used  to  pre-paro 

'  13) 

the  data  In  Table  X  In  order  to  correlate  the  data  given  by  Klcnent' 


and  ehovn  In  Tabic  XX. 
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Figure  9.  Taui  Activity  of  Soil 


TABLE  IX 


Radioactivity  In  3olla  (Mc/MT) 


Hue lid* 


0 

1  hr 

1  day 

1  month 

1  year 

IW* 

1620 

1430 

542 

a 

a 

m/7 

504 

7.05 

a 

Al28 

83,670 

15.9 

a 

Si31 

991 

773 

1.78 

a 

P32 

1.21 

1.21 

1.15 

0.276 

a 

J4 O 

K 

a 

204.5 

194 

50.0 

a 

Ca** 

0.300 

0.300 

0.300 

O.263 

0.057 

cU9 

1100 

135 

a 

Ti51 

156 

0.0015 

a 

Mn56 

21 22 

1655 

3.62 

a 

P*55 

0.110 

0.110 

0.110 

0.104 

0.087 

Pe59 

0.014 

0.014 

0.014 

0.0095 

0.00006 

a  •  less 

than  one  curie  per  megaton 

The  formula  for  the  radioactivity  induced  iu  any  element  as  given 
by  Senf tie  and  Champlon^"^  isi 

A.  -  (o  flfh)(l  -  e'0,693  t/T)(e-0.693  »/*)  (72) 


where 


activity,  In  disintegrations  per  cccond,  after  the  nuclide 
has  been  removed  from  the  flux  for  a  period  0 
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SLANT  RANGE  FROM  EXPLOSION  (YARDS) 

Figure  10.  Integrated  Neutron  flux  as  a  Function  of  the  Slant  Kan go  la  Air 
of  0.9  Sea-Level  Density  for  a  1-Kllcton  Explosion 
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As  a  special  case,  -where  t  is  toova  to  he  less  than  15S&  of  T, 
a  further  simplified  formula  was  used. 

\  •  a:  '2=|2)  {.-0-693  0/,r)  (75) 

Values  used  for  Ar  obtained  by  Senftle  and  Champion^^  are  listed 

below; 


rtl. 

at 

Ha 

-  1-39  x  13i0 

_  M» 

Ca 

-  1.73  x  10® 

Ml27 

-  1.3  x  lY 

ca49 

.  2.54  x  107 

Al2® 

-  4.79  x  'Jp 

ti5L 

-  9-01  x  107 

SL'1* 

-  6. 74  x  10 7 

-  I.39  x  10U 

F32 

-  3-69  x  iP 

Fe55 

-  4 . 56  x  10® 

K* 

.  1.32  x  Up 

Fe^^ 

-  2.26  x  107 

In  calculating  the  data  of  Table  X  ,  a  neutron  flux  of  0.91  x  101® 
r.uutrons/sq  ee  per  megaton  via  assumed.  The  total  weight  of  the  coil  involved 
was  taken  us  5.92  x  10^  grssu.  IMs  value  was  calculated  using  a  2500  yd 
rudlu3,  a  l*ft  depth,  an  a verige  density  of  1.18  gni/cu  cm,  and  with  a  soil 
composition  as  that  given  !c  Table  VIH.  The  irradiation  time  (t)  who  taken 
03  or.e  minute. 

The  data  presentee  it  Table  X  agree  rcaeonubty  well  with  the 
data  presented  by  Klement  la  Table  IX.  This  agreement  is  within  an  expected 
factor  of  15  ss  the  neutron  flux  used  vac  a  factor  of  15  higher  than  that 
assumed  by  Klecvnt. 

Knnleviile^^  sta'.*s  that  the  radioactivity  of  the  fission  products 
from  u  '-DcL'.f.r  explosion  overrun  lows  the  InduceJ  radioactivity  by  a  fuctor  of 
approx imatvlo-  10®.  Ibren  tho^a  this  tends  to  make  the  contribution  of  the 
induced  radioactivity  almost  r.j.-liglble  there  are  conditions  under  which  the 


induced  radioactivity  night  play  a  detectable  part.  The  induced  nuclides 
that  form  soluble  salts  will  be  readily  dissolved  In  both  surface  and 
ground  water  and  therefore  will  add  to  the  contamination  of  the  water 

supply. 

It  has  been  generally  recognized  that  the  activating  neutrons 
will  penetrate  the  roil  before  the  blast  wave  arrives.  Some  of  the 
activated  coil  will  then  be  token  up  into  the  cloud  due  to  the  cratering 
effect  and  thereby  becomes  part  of  the  total  fallout.  According  to 
Dr.  C.  F.  Miller  ^  the  induced  activity  in  the  coil  will  be  about 
0.019  (l«9  pei  cent)  of  the  ionization  rate  due  to  normal  fission  (lfp(t)] 
1ht3  fraction  (fill  be  a  small  percentage  of  the  total  ionization  rate. 
TMo  estimate  of  the  induced  activity  is  applicable  for  n  10 Of,  risslon 
yield,  assuming  0.8  neutron  capture  per  ficslon  from  U-238,  rcgnrdleoa 
of  weapon  size. 

It  has  also  become  established  that  the  induced  activity  will 
predominate  for  a  total  fusion  bomb.  For  a  JO £  fusion  and  jjOjt  fission 
bonto,  the  induced  activity  will  be  relatively  small,  except  for  a  short 
period  Immediately  after  detonation.  However,  activated  6odlum  and 
mn^tcalua  may  contribute  materially  at  a  later  time^). 

It  la  therefore  concluded  that  the  contribution  of  Induced 
radioactivity  to  the  contamination  of  voter  supplies  will  not  be  signi¬ 


ficant. 
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V.  TRANSPORT  OF  FALLOUT  PARTICLES  FT  SURFACE  WATER 

Transport  of  the  insoluble  portion  of  fallout  may  be  divided 
into  four  general  phases:  transport  in  air*  transport  in  overland  flow, 
transport  in  stream  flow  and  transport  In  reservoirs.  Each  of  the  four 
phases  may  be  treated  as  a  separate  topic.  Of  the  latter  three,  only 
transport  by  stream  flow  nas  been  considered  in  any  detail)  overland 
transport  and  reservoir  transport  are  mentioned  only  briefly.  Transport 
in  air  is  discussed  in  Section  II  of  this  report. 

There  are  three  major  limitations  on  any  stream-flow  transport 
analysis.  These  are:  (l)  the  complex  nature  of  sediment  behavior  in 
streams;  (f?)  the  paucity  of  data  concerning  fallout,  especially  data 
concerning  particle  size  distribution  or  the  iclatlon  of  radioactivity 
to  particle  size;  and  (3)  the  obvious  impossibility  of  obtaining  utrcoa- 
f low  data  for  every  reach  of  stream.  These  limitations  precluded  an 
analysis  of  the  quantity  of  sediment  resulting  from  the  fallout  of  a 
given  burst,  although  work  on  the  quantity  of  fallout  la  now  being  done. 
It  could  not  be  determined  that  a  certain  weight  of  sediment  would  fall 
on  a  water  surface  or  what  the  particle  size  distribution  would  be. 

Computations  arc  available,  however,  from  Glass  tone  ^  and 
the  OCD  fallout  model  which  gives  certain  relationships  between  particle 
sizes,  downwind  ulctancc  carried,  and  total  activity. 

Glasstone  chows  that  99  per  cent  of  the  total  activity  of 
fallout  is  carried  by  particles  of  ^00  nieivno  or  less  in  diameter. 

This  relation  of  particle  size  to  total  activity  was  derived  from 
Figure  9-18?  of  CLaestone's  boo)/^^  and  is  chcrvn  as  Figure  11  of  this 
report.  The  method  contained  in  this  section  is  based  on  this  relation¬ 
ship. 
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Particles  of  bC 0  nieron  diameter  are  on  the  order  of  magnitude 
of  normal  stream  sediment*  This  means  that  the  particles  carrying 
essentially  all  the  activity  will  remain  in  suspension,  or  at  least 
will  net  settle  out  immediately.  It  follows  that  most  of  the  fallout 
particles  which  land  on  the  water  surface  or  find  their  vay  into  the 
stream  from  overland  transport  will  be  carried  by  the  stream  and  will 
present  a  hazard  at  6ome  distance  downstream. 

A  relationship  of  particle  size  versus  downwind  distance 
curried  was  computed  In  Section  II  of  this  report.  (Table  VIl).  Data 
from  that  table  has  been  plotted  in  this  section  as  Figure  12. 

Hue  following  general  assumptions  were  made: 

(1)  Uptake  by  biological  organisms  la  negligible. 

(2)  The  stream  chemistry  is  not  such  as  te  dissolve  ordinarily 
insoluble  particles. 

(3)  Effect  of  flocculation  on  eettllng  rates  is  negligible. 

(1)  Weapon  size  (5MT),  wind  speed  (l5mph),  arl  fallout  pattern 

are  used  to  correspond  with  other  cectlons  of  this  report. 

(5)'  Radioactivity  or  particles  over  500  microns  may  be  ignored. 

Further  assumptions  are  made  and  discussed  in  the  text.  Some 
of  these  asc  unptlons  depend  on  the  nature  of  n  specific  stream,  some 
require  further  investigation,  ond  rome  which  ore  considered  negligible 
in  an  approximate  analysis  may  assume  a  proportionally  greater  importance 
In  a  more  detailed  computation. 

An  estimate  must  be  made  of  stream  width,  volume,  and  surface 
area.  Ibis  cuiy  be  done  by  equating  the  cross-sectional  area  of  the 
.it ream  (at  &  gaging  station)  to  a  rectangle  of  equivalent  area  whose  length 
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Figure  11,  Total  Activity  Percentage  in  Relation  to  Particle  Diameter 
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figure  12.  Particle  Diameter  versus  Dovuvind  DlGtaocc  (MaxIbub  and  Minimum  Curves) 
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is  equal  to  the  strean  width.  nils  ia  shown  in  Figure  13*  length 

of  the  stream  between  any  two  sections  ic  considered  to  be  a  straight 
line  equal  to  the  distance  along  the  6treaa  channel  between  the  two 
sections.  Each  section  of  the  stream  nay  then  be  considered  to  be  a 
truncated  pyramid  lying  on  its  side  with  its  bases  the  equivalent  rec¬ 
tangles.  (Or  for  the  section  containing  the  source,  a  complete  pyramid, 

since  the  source  is  considered  to  te  a  point.)  Concentrations  nay  be 

(18) 

found  using  the  volume  of  the  pyramids. 

Intensity  (roentgen/hour)  contours  for  a  given  burst  are  super¬ 
imposed  on  the  watershed.  The  Intensities  may  then  be  converted  to  activities 
(atorac/sq.  ft.)  for  insoluble  (and  soluble)  portions  of  the  various  lcotopes 
under  consideration  by  use  of  Table  VI.  Details  of  this  method  are 
discussed  on  page  95  through  page  98  of  this  report. 

If  fallout  landing  on  the  water  surface  is  acoumed  to  mix  uniformly 
within  the  reach  of  channel  enclosed  by  a  given  contour,  an  activity  con¬ 
centration  <aay  be  computed  using  the  idealized  channel  shape.  This  concen¬ 
tration  will  move  downs t re nn  as  a  unit,  and  the  time  of  arrival  and  time  of 
lapse  of  this  concentration  may  be  found  at  a  point  downstream  by  concideiirg 
the  particles  to  move  at  mean  stream  velocity. 

By  plotting  the  various  concentrations  from  unlto  under  given 
intensity  contours  aa  ordinate  against  time  as  abcioca,  a  cerlea  of  horizontal 
bars  lo  obtained.  These  may  be  averaged  graphically  to  obtain  a  time- 
conccntrotion  curve  for  the  otrervn.  !tote  that  concentration  is  exprcc.-.c-d 
as  activity  rather  than  particle  concentration. 

By  use  of  Figure  12  a  maximum  particle  size  at  a  given  point 
downwind  from  the  blast  may  be  found.  The  maximum  size  particle  landing 
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oa  the  stream  reach  under  consideration  is  compared  to  the  maximum  particle 
size  which  the  stream  will  carry  In  suspension.  If  the  maximum  size  particle 
1»  carried  then  it"  la  to  he  assumed  that  smaller  particles  will  also  he 
carried,  and  the  method  remains  valid. 

If  the  stream  will  not  carry  the  maximum  size  particle,  then 
some  percentage  of  the  activity  must  be  regarded  as  going  into  bed  load, 
which  travels  at  a  much  6lcwer  rate.  The  only  way  now  available  to  find 
the  proportion  of  activity  which  would  go  into  bed  load  is  by  assuming 
that  if  particles  of  a  given  size  are  found  to  settle,  that  all  particles 
of  this  size  or  larger  will  settle  out,  and  all  smaller  particles  remain 
la  suspension.  The  relationship  of  particle  size  to  activity  given  by 
Glasstone^1*^  may  then  be  used  to  find  what  percentage  of  the  total  activity 
would  settle  out  with  the  larger  particles. 

lfce  ability  of  a  stream  to  carry  particles  of  a  given  size  or 
smaller  is  best  defined  by  what  may  be  termed  the  effective  carrying 
velocity  of  the  stream,  i.e.  that  velocity  below  which  settling  of  particles 
of  a  given  olze  is  most  likely  to  occur.  The  effective  carrying  velocity 
is  obviously  not  related  to  the  idealized  stream  channel  used  to  compute 
concentrations,  but  is  rather  a  function  of  both  channel  characteristics 
(actual  channel  width  and  depth)  and  flow,  as  well  as  particle  size  and 
sediment  distribution.  Stream  flow  data,  no  matter  how  accurate  or  com¬ 
prehensive,  cannot  cover  every  mile  of  channel.  Conditions  which  will 
produce  the  greatest  possibility  of  settling  must  therefore  be  generaliza¬ 
tions  of  what  is  known  of  the  nature  of  an  individual  stream. 

Particle  Gize  distribution  in  air  is,  unfortunately,  not  known, 
except  for  maximum  and  minimum  distance  which  a  particle  of  given  size  will 
travel.  There  is  no  way  at  present  to  compute  the  number  of  particles 


per  unit  area  landing  on  the  watershed.  Sediment  concentrations  cannot, 
therefore,  be  computed.  ThlB  is  the  reason  for  computing  concentration 
in  terns'  of  activity.  Additional  fallout  data,  particularly  on  particle 
distribution,  and  additional  study  or.  sedirr.cnt  concentrations  are  needed 
to  support  the  assumptions  that  the  fallout  particles  mix  .ir..u'or«ly, 
that  particles  will  go  out  of  cur pens ion  according  to  size,  and  that 
particleo  irove  downs t ream  at  mean  stream  velocity. 

In  order  to  approximate  actual  stream  conditions,  close  scrutiny 
of  m'-.ps  ar.d  existing  stream-flow  records  vill  be  required.  Large-scale 
topographical  mans  should,  yield  much  information  on  channel  conditions 
and  stro'-m-bed  slopes.  The  U.  S.  ecological  Surveys  and  Corp3  of  Engineers 
publish  or  have  on  file  not  only  gaging-etatlon  records  but  much  special 
measurement  data. 

So  analysis  has  been  made  for  transport  resulting  from  overland 
flow.  Fallout  landing  on  the  land  portion  of  the  watershed  is  conoidcrad 
less  of  an  Immediate  problem  than  that  lending  on  the  water  surface, 
and  will  probably  present  no  Immediate  problem  unless  It  Is  raining 
at  time  of  fallout  arrival,  (in  which  case,  incidentally,  the  fallout 
nvvy  be  intensified. )  Subsequent  rain  may  produce  a  second  activity 
peak  st  a  water  intake  by  washing  particles  Into  the  stream.  It  Is  of 
interest  to  note  that  at  Covcetn  Hydrologic  Laboratory  It  war.  found  that 
much  of  the  runoff  which  was  formerly  considered  to  flow  overland  actually 
gcc-u  under  the  surface,  or  alternates  between  sub-surface  and  above -surface 
fj  -•  In  work  with  radioactive  tracers  at  Ooweeta,  great  difficulty  was 
experienced  lr.  finding  a  tracer  which  would  follow  the  water  at  the  came 
ijved  r.a  the  water.  Most  of  the  Isotopes  used  at  Coveeta  tended  to  otay 
near  U.e  point  where  they  were  applied.  These  studies  were  carried  out 


75 


on  steep  slopes  on  which  there  was  a  higher  percentage  of  surface  runoff 
than  on  the  same  soil  at  a  lesser  slope;  however,  the  Isotopes  used  were 
not  ah  sorbed  within  partlclea.^^^ 

Partielts  which  go  into  bed  load  nay  present  a  long  term  hazard 
because  of  the  slower  and  more  Irregular  movement  of  bed  load.  A  bed 
load  analysis  will  be  required  If  It  Is  found  that  large  particles,  In 
spite  of  their  possession  of  only  a  small  percentage  of  the  total  activity, 
are  found  to  actually  carry  a  dangerous  amount  of  activity.  The  smaller 
the  particle  size  that  Is  found  to  settle,  the  more  significant  becomes 
the  bed  load  analysis.  Note  that  such  a  bed-load  analysis  complements 
assumption  5  on  page  86. 

Transport  In  reservoirs  Is  more  likely  to  produce  settling 
because  reservoirs  are  more  likely  to  approach  quiescent  conditions. 
Calculations  will  probably  be  simplified  because  the  conditions  of  reser¬ 
voir  transport  more  closely  approaches  the  condition  of  settling  in  a 
quiescent  basin  than  do  the  other  phases  of  the  transport  problem. 

This  method  was  presented  not  only  because  It  is  the  one  most 
likely  at  present  to  produce  usable  results,  but  because  It  provides  an 

outline  of  the  problem  of  stream  transport.  Assumptions,  questions  as 
to  their  validity,  and  areas  requiring  further  study  have  been  pointed 
out.  Computations  have  not  been  Included  because  it  Is  wished  to  ascer¬ 
tain  the  validity  of  some  of  the  assumptions,  and  because  sufficient 
study  has  not  been  made  on  specific  streams.  This  treatment  therefore 
constitutes  only  a  good  start;  however,  further  study  In  the  same  manner 
should  produce  rcoults  of  known  validity. 
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71.  EVAIUATICS  C?  IffSTTCIPAI.  WATER  StTCTIX  CCCuAMMAITIGB 
A.  Introduction 

An  approach  to  the  problem  of  water  contamination  by  radioactive 
fallout  was  taken  which  would  give  a  no  re  reliable  set  of  concentration 
values  than  those  calculated  previously.  The  values  presented  In  Quarterly 
Technical  Report  Bo.  -  -were  maximized,  in  that  they  represented  a  case  in 
which  the  highest  c  once  at  rat  ion.  that  could  reasonably  be  expected  at  H  +  1 
hour  froa  any  kind  of  attack  or  envlronaent,  except  for  possibly  a  very  high 
megatonnage  attack.  PoLlovizg  is  presented  a  more  detailed  and  realistic 
approach  to  the  problen.  The  values  herewith  reported  represent  specific 
caces  in  which  envirosneattl  and  other  factors  have  been  considered,  and 
It  is  believed  they  are  at  least  of  the  came  order  of  magnitude  as  vould 
be  expected  following  a  nuclear  attack. 

Only  tha  cities  of  Houston,  Texas  and  Beu  York  were  chosen  for 
a  full  scale  evaluation.  Since  the  values  calculated  for  these  two  water¬ 
sheds  were  in  good  correlation,  it  was  felt  that  the  results  could  probably 
be  applied  without  serious  error  to  the  other  cities  previously  considered. 
Excellent  watershed  data  was  available  on  both  cities,  thus  facilitating 
the  calculations  cons  lie  rally. 

A  physical  Integration  waa  performed  over  the  considered  watersheds 
and  reservoirs  to  give  reaJLirtic  specific  Isotope  concentrations  at  H  +  1 
hour  for  variously  directed  lj  aph  winds.  As  before,  tha  case  of  fallout 
contamination  from  that  falling  in  the  reservoir  alone  was  considered  as 
well  as  the  case  of  rmoff  ctatani nation  free  the  entire  watershed. 

B.  Calculations 

Since  large  scale  naps  of  both  the  Houston  and  New  fork  watersheds 
were  available,  transparent  or* -lays  were  superimposed  over  the  watershed 
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ureas  tr.d  appropriate  grids  scaled  ir..  Thete  grids  consisted  of  five  mile 
i-tj'-i  ire  a  for  water  :hed  areac,  ar.d  two  xlle  c  qua  raj  for  the  reservoir  areas 
when  large  es.o>igh  to  apply.  Figure  l.h  bfcovc  the  grid  on  the  watershed 
serving  hcu-ito a-  The  gro.ind  -^ro  of  a  f-  M?  weapon  detonation  via  located 
so  as  to  give  a  fnllo'-t  pa.hV.--n  which  fell  on  the  watershed  for  the  15  mph 
wind  .-colei.  Vr-.riow-  wind  directions  wore  also  chot.-n  so  as  to  give  some 
bvis  for  comparison  and  to  account  for  passible  seasonal  variation.  It 
is  obvious  fron  the  geography  of  How  York  and  Kou-uton  that  the  wind 
directions  frer  the  east  and  north  would  give  rise  to  a  negligible  aiourt 
of  fi.ilort  over  shell-  respective  watersheds,  so  that  these  wind  directions 
were  twitted. 

The  attacks  that  were  evaluated  were  U3  follows; 

Case  la  Houston  -  Crovr.d  zero  13  downtown  Houston  ur.d  the  wind  direction 
is  fro*  the  south. 

Case  Th  Houston  -  Ground  zero  la  downtown  Houston  tr.d  the  vlr.d  direction 
Is  from  the  vest. 

Case  I  la  hew  'fork  -  Ground  zero  is  Central  Park  and  the  wind  direction 
is  from  the  south. 

Case  lib  Hew  York  -  Ground  zero  is  Binghamton,  Kev  York,  ur.d  the  wind 
direction  is  from  the  vest. 

tv  pirate  evaliiat  Ions  were  made  for  ench  of  these  cases  with  respect 
to  cor.  feci  nan  Ion  from  tne  reservoir  and  that  from  runoff  of  the  entire  water¬ 
shed.  T:i  the  case  of  Houston,  the  concentration  of  activity  was  calculated 
for  one  ol  the  f-.ed  str:ons  i'rsr.  direct  fallout  co:,t. urination  for  comparison 
purposes. 

For  the  reservoir  study,  the  ^re-a  of  the  respective  reservoirs 
within  each  two  mile  grid  v»>  found  by  piacireter.  In  the  Houston  watershed 


Figure  14.  Method  of  Graphical  Integration  Used  for  the  San  Jacinto 

Watershed  Serving  Houston,  Texan 


there  Is  only  one  reservoir  and  in  the  New  York  watershed  all  six  rcoervoir* 
empty  into  one  so  that  the  problem  is  somewhat  simplified  if  complete  mixing 
is  assumed.  The  centroid  of  each  selected  grid  area  wsuj  estimated,  and  lta 
coordinates  measured  with  respect  to  ground  zero.  An  Intensity  was  then 
assigned  to  each  set  of  coordinate^  and  it  was  assumed  to  be  constant  over 
that  specific  area.  The  intensity  values  were  chosen  according  to  whether 
the  specified  area  was  either  covered  by  the  upwind  or  downwind  oif'ud. 

Those  values  related  to  the  upwind  cloud  were  arrived  at  by  interpolation 
between  previously  estimated  intensity  contours,  and  those  of  the  downwind 
cloud  were  arrived  at  by  interpolation  of  computer  calculated  contours. 

To  convert  from  the  intensity  over  each  area  to  activity  in 
atoos/sq  ft  in  each  area,  the  value  of  N*(A)/j^j  from  the  fallout  model, 
for  each  isotope  considered,  was  taken  as  its  value  corresponding  to  the 
downwind  distance  X.  By  multiplying  the  appropriate  n*(a)/j^j  value  by 
the  intensity  aver  each  square,  the  concentration  in  atomc/oq  ft  is  found 
for  each  area.  3y  then  multiplying  etch  atom  concentration  by  its  cor¬ 
responding  area,  the  total  number  of  atoms  in  each  square  of  the  grid  is 
found.  The  total  number  of  atoms  of  each  isotope  at  K  ♦  1  hour  in  the 
reservoir  v&3  then  found  by  adding  up  the  contributions  from  each  square 
over  the  reservoir.  By  assuming  that  complete  mixing  has  occurred,  the 
concentration  of  activity  in  the  reservoir  in  atoros/litcr  may  be  obtained 
by  dividing  by  the  total  volume  of  water  in  the  reservoir.  The  concen¬ 
tration  of  each  isotope  considered  at  H  ♦  1  hour,  for  the  Houston  and 
How  York  reservoirs  is  shown  in  Tables  >1  and  XT  along  with  the  equivalent 
activity  in  nc/al. 
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Table  XI 

Activity  Concentrations  for  Direct  CootaMnotioa 
of  the  Houston  Reservoir 


Couth  Wind 

West  Wind 

Isotope 

Atoits/iiter  uc/ir.l 

Atoms/ liter  uc/ml 

Sr-89 

^•6  x  1C5  1.9  x  10*5 

b.6  x  109  1.9  x  IO'5 

Sr-90 

10.9  x  IO9  2.3  x  10'7 

10.9  x  109  2.3  x  IO'7 

Ru-106 

5-1  x  1C9  3.0  x  10*6 

5.1  x  109  3.0  x  IO'6 

1-131 

9.8  x  IO9  2.7  x  10*U 

9-8  x  109  2.7  x  10'* 

Cs-137 

5-3  *  109  1.1*  x  10*7 

5.3  x  109  1.1*  x  IQ'7 

Ba-li*0 

13.7  x  109  2.3  x  10'“ 

13.7  x  IO9  2.3  X  10'* 

Table  XII 

Activity  Concentrations  for  Direct  Contamination 
of  the  New  York  City  Reservoirs 


■I 

South  Wind  I 

West  Wind 

Atoms/ liter  uc/ml  1 

Atoms/liter  uc/al 

Sr-89 

2.9  x  1010  1.2  X  IQ-* 

3.9  X  IO10  1.6  x  10"* 

Sr-90 

5.6  x  1010  1.2  x  10’6 

6.9  x  IO10  1,5  x  IO*6 

Ru-106 

2.5  x  10L0  1.5  x  io*5 

3.0  x  IO10  1.8  x  IO'5 

I-13I 

5.6  x  1C10  1.5  x  10'3 

6.8  x  IO10  1-9  x  IO’3 

Cr-137 

3.7  x  1010  7.3  x  10’7 

5.1  x  IO10  1.0  x  IO'6 

Ba-ll*0 

7.0  x  1C10  1.2  x  iO*3 

8.1*  x  IO10  l.U  x  IO'3 

.100. 


Tke  activity  concentrations  for  the  Houston  reservoir  (Table  XT)  are 
the  earne  for  south  wind  and  vert  wind  beenuoe  the  come  intensity  contour  falls 
across  the  center  of  the  reservoir  in  Doth  caoes.  When  the  wind  is  from  the 
south,  the  right  hand  portion  of  the  contour  crosses  the  reservoir.  With  wind 
from  the  west,  the  left  hand  pt rtion  of  the  contour  is  involved, 

A  separate  evaluation  of  stream  eontamlnution  vac  mode  using  a 
randomly  selected  creek  feeding  the  Houston  reservoir  for  the  purpose  of  com¬ 
parison  only.  This  selected  fred  stream  (Cyprce3  Creek)  ran  perpendicular  to 
the  downwind  fallout  axis  and  emptied  into  the  Houston  reservoir.  The  evalua¬ 
tion  was  carried  out  in  the  rare  manner  ao  that  done  on  the  reservoir.  The 
stream  va3  assumed  to  approximate  on  isocelco  triangle  with  the  base  repre¬ 
senting  the  width  of  the  streoa  south.  Two  mile  lengths  were  laid  off  and  the 
trapezoidal  areas  cnlculuted.  Intensity  values  and  activity  conversion  values 
were  chosen  as  before  to  obtaim  the  total  activity  falling  directly  on  this 
strewn,  To  present  an  estimation  of  the  concentration  of  activity  In  the 
stream,  complete  mixing  was  as  rosed.  However,  it  is  evident  that  the  activity 
will  really  move  In  a  elug  type  of  formation,  completely  mixing,  if  at  all, 
only  when  the  radioactivity  retches  the  reecrvoir.  Although  some  mixing  during 
stream  flow  will  occur  due  to  turbulence  and  diffusion,  the  presented  case 
will  never  be  actually  attained,  since  countercurrent  mixing  is  unlikely. 

From  the  presented  values  in  Title  XIII  it  may  be  seen  that  some  build-up  of 
activity  may  occur  in  areas  iirte-1  lately  surrounding  the  stream  fecd-ln  points. 
However,  this  build-up  will  r.ot  represent  a  curious  problem  in  the  studied 
cace  because  of  its  distance  from  the  municipal  intake,  and  the  relative  time 
elapsed  from  the  time  oC  dctor.v.  Ion  to  the  intake  time. 

The  second  case  studied  was  that  of  the  contamination  to  be  expected 
assuming  runoff  from  the  entire  watershed.  Tills  stuuy  was  made  for  both  the 
New  York  and  Houston  watersheds,  the  method  being  similar  to  that  of  the  first 
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Table  XIH 

Activity  Concent  rations  for  Cypress  Creek  (South  Wind) 


Isotope 

Atoms/ liter 

jir/  rr). 

Sr-89 

2.2 

X 

1012  * 

9.4  x  10*3 

Sr-90 

5-4 

X 

1012 

1.1  x  10 

Ru-106 

2.6 

X 

10* 

1.6  x  10“3 

1-131 

4.4 

X 

10* 

1.2  x  10*L 

Cs-137 

2.5 

X 

1012 

4.9  x  10"5 

Ba-luO 

6.8 

X 

10 12 

1.2  x  10"1 

case.  Five  mile  grids  were  used,  however.  Instead  of  two  mile  grids.  Wind 
conditions  were  selected  Identical  to  those  of  Case  I.  To  calculate  the  actual 
amount  or  radioactivity  (soluble  fraction  only)  that  reaches  the  reservoir, 
the  runoff  coefficients  supplied  by  the  respective  municipal  voter  works  were 
assumed  valid  In  that  they  represent  a  maximum. value  to  be  expected.  The 
calculated  activities  are  therefore  considered  to  be  maximum  values.  The 
actual  radioactive  runoff  coefficient  will  actually  be  lesB  than  the  aqueous 
runoff  coefficient  depending  on  (l)  the  Instantcous  moisture  content  of  the 
soil,  (2)  the  duration  of  time  froa  detonation  to  rainfall,  (3)  lon-exchange 
and  absorption  in  and  on  the  coll,  and  (4)  plant  uptake,  Thv  calculated 
values  of  activity  concentrations  of  the  selected  lsotopco  at  H  ♦  1  hour, 
assuming  complete  dilution  by  the  composite  likes  and  streams  In  the  watershed 
for  both  Houston  and  Hew  York  Is  presented  In  Tables  XIV  and  XV.  The  factors 
Involved  obviously  vary  widely  from  one  environment  to  another,  so  that  general 
assumptions  made  from  ony  particular  reservoir  should  ho  applied  with  caution 


to  other  situations 
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Table  XIV 

Activity  Conccr.t  rat  loss  for  Contamination  from  Runoff 
for  the  Houston  Watershed 


South 

Wind 

Vest 

Wind 

Isotope 

Atom/liter 

nc/ -1 

Atons/liter 

pc/ml 

Sr-89 

4.5  x  1012 

1.9  x  10‘2 

3.6  x  1010 

1.6  x  10*4 

Sr-90 

8.0  x  10 12 

1.9  x  10*4 

3.5  x  1010 

1.8  x  10*6 

Ru-106 

4.1  x  1012 

2.5  x  10*3 

4.0  x  1010 

2.4  x  10*5 

1-131 

9.2  x  1012 

2.5  x  10*1 

8.4  x  1010 

2.3  x  10*3 

Cs-137 

12 

f  5*3  x  10^ 

1.0  x  10‘4 

4.3  x  1010 

8.5  x  lu*7 

Ba-llO 

12 

11.3  « 

1.9  x  10‘1 

.  11,0  x  id10 

1.9  x  10*3 

Table  XV 


Activity  Concent  rations  for  Contamination  from  Runoff 


for  the  Key  York  City  Watershed 


South  Wind 

Wcct  Wind 

Ioot.ope 

Atoma/llter  Mc/ril 

Atomo/llter  uc/ral 

Sr-89 

0,36  x  1012  1.5  x  10*3 

1.6  x  1012  6.8  x  10"3 

Sr-90 

0.65  x  1012  1.4  x  10  5 

2.8  x  1012  5.9  x  10*5 

Ru-106 

0.29  x  1012  1.7  x  10’4 

12  -4 

1.2  x  10  7.1  x  10 

1-131 

O.65  x  1012  1.8  x  10*2 

2.8  x  1012  7.6  x  10*2 

Co-137 

0.48  x  1012  9-5  x  10*6 

2.2  X  1012  4.4  x  10'5 
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C.  Cone  las  Iona : 

It  la  not  the  purpose  here  to  draw  any  concLualona  In  regard  to 
the  biological  hazard  resulting  from  the  concentrations  of  activity  presented 
In  Tables  XI  -  XV.  This  aspect  is  being  considered  elsewhere  In  this  report 
and  also  by  other  investigators. 

It  may  be  prudent,  however,  to  consider  the  relative  value  of 
the  previous  calculations.  It  must  he  remembered  that  only  a  5  KT  weapon 
was  considered  and  that  all  concentrations  arc  for  H  ♦  1  hour.  Thus 
for  any  other  size  weapon,  appropriate  revisions  in  the  fallout  model  will 
...ive  to  he  made  before  it  may  be  applied.  As  for  a  multi -bomb  attack, 
one  cannot  but  accept  some  method  based  on  additivity.  Complications 
arrive  here,  aa  far  aa  forming  a  model  la  concerned,  becauao  of  ground 
zero  locations.  Since  H  ♦  1  hour  values  of  concentration  nro  given,  one 
may  apply  suitable  decay  curves  to  adapt  them  to  any  specifio  otroam  flow, 
reservoir  flow-throu£>,  or  intake  time,  etc.  for  a  particular  watershed 
aystem. 

It  Is  not  to  be  Implied  here  that  the  Houston  and  New  York  values 
nre  representative  of  all  municipal  reservoirs.  However,  they  do  present 
a  reasonable  value  to  be  expected  In  a  rcaliotlc  case,  and  also  demonstrate 
the  ease  with  which  an  evaluation  may  be  made.  By  making  ouch  an  Integrated 
evaluation  for  any  specific  watershed  and  applying  local  conditions,  the 
relative  radiological  hazard  may  be  estimated. 

The  Houston  and  flow  York  values  will  be  seen  to  vary  from  the 
potential  valueo  If  unusual  conditions  prevailed  such  as  would  occur 
during  a  dr/  spell.  However,  the  effects  vlll  be  sonewhat  counterbalancing 
in  that  high  reservoir  concentration  due  to  lav  volume  will  bo  somewhat 
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offset  by  the  lack  of  runoff  contamination.  It  la  doubted  that  unnatural 
conditions,  unless  very  severe,  will  change  the  concentration  of  activity 
by  no  re  than  an  order  of  nagnltwie. 
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vn.  ATIALVSIS  OF  RADIOACTIVITY  lit  WATER 

A.  Radiochemical  Methods 
1.  Introduction 

In  most  radiochemical  analyses  the  steps  are  concentration, 
separation,  purification,  and  counting.  Evaporation  and  Ion  exchange 
are  the  evo  principal  methods  of  concentration,  but  may  not  be  necessary 
If  the  level  of  activity  Is  aufflclcntly  high.  The  individual  radio¬ 
nuclides  sre  usually  separated  by  means  of  precipitation,  solvent 
extraction,  or  Ion  exchange  elution.  Purification  of  the  Individual 
radionuclides  usually  follows  stan'lard  chemical  procedures  vhlch  depend 
mainly  on  the  solubility  characteristics  of  compounds  of  the  radionuclides* 
Coryell  and  »<garman  have  described  (urlf lcatlon  of  many  of  the  more 
Important  radionuclides.  Radloicotopta  are  counted  using  appropriate 
instruments,  depending  on  the  specific  Isotope  and  the  desired  accuracy. 

"Carriers"  are  of  prime  importance  In  vet  radiochemical  procedures. 
The  mechanisms  end  usage  of  radioisotopic  carriers  have  been  veil  defined 
by  Overman  and  Clark  A  carrier,  vhlch  is  normally  the  stable  isotope 

of  the  clement  being  determined  or  a  stable  element  vlth  very  slmllaur 
properties,  lo  added  to  lncrcace  the  total  concentration  of  the  element, 
because  the  radioactivity  vlll  normally  be  present  In  micro  (juantltles  only. 
Hence,  if  the  carrier  is  not  added,  part  of  the  radioactive  material  will 
usually  oe  lost,  even  though  employing  the  moat  scrupulous  analytical 
techniques.  By  weighing  the  final  solid  to  be  counted,  and  knowing  the 
original  nnount  of  carrier  added,  the  percentage  chemical  and  physical  loaa 
r.ay  be  calculated.  By  assuming  complete  interchange  of  carrier  and  radio* 
nuclide,  the  came  percentage  of  each  vlll  be  lost,  and  the  original  concen¬ 
tration  of  activity  may  readily  be  calculated.  The  carrier  itself  Is 


normally  added  as  a  solution  of  a  soluble  salt  of  carrier  element 
prior  to  any  chemical  procedure.  Carrier  addition  la  used  extensively 
throughout  the  described  analytical  procedures,  although  scoe  carrier-free 
determinations  are  now  cooing  Into  use. 

There  are  several  compilations  available  vr.tch  describe  detailed 

(20  22  23  24} 

procedures  for  the  quantitative  determination  of  many  radionuclides  v  '  *  ' 

tost  of  these  methods  have  been  proven  and  are  In  general  use. 

2.  Procedures  for  the  Determination  of  Specific  Eallcnuclldea 
a.  Radtostrenttua 

Strontium -90  has  been  generally  recognized  as  being  a  very  bio¬ 
logically  hazardous  radionuclide  as  It  accumulates  In  bone.  IXte  to  its 
relatively  dangerous  nature,  vide  lute  rent  has  been  aroused  to  find  a  simple 
and  rapid  method  of  analysis  for  otrontlua-90.  Sb.'irlques  that  have  been 
In  use  for  several  years  are  Still  finding  widest  use  because  of  their 
accuracy,  the  time  consideration  being  less  important  la  peacetime. 

Strontium -69 >  perhaps  somewhat  less  of  a  hazard  than  strontlun-90, 
because  of  Its  shorter  half-life,  will  also  be  present  In  fallout  contaminated 
water.  Separate  analyses  of  strontlum-89  and  street ium-90  are  thus  desirable. 
Strontlum-89  activity  Is  normally  found  by  subtraction  of  strontlum-90  from 
total  radlostrontlum  concentration. 

(25) 

(1)  Radlostrontlum  and  Radlobarlua  by  Mtrate  Separation* 

Strontium  and  barium  carriers  are  added  v.d  the  Croup  II  cations 
are  precipitated  as  the  carbonates.  Partial  separation  from  calcium  Is 
accomplished  by  nitrate  precipitation  In  fuming  nitric  acid.  The  remaining 
calcium  nitrate  Is  then  extracted  with  acetone.  Rare  earths  ar.d  other 
trlvaler.t  cations  are  removed  by  t  vo  hydroxide  ocawr-glngs.  Barium  Is 


finally  separated  fron  strontium  by  precipitation  barium  chromate.  Stron¬ 
tium  Is  then  collected  as  the  oxalate  and  counted.  The  barium  Is  converted 
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f/'oin  the  chromate  to  the  chloride  for  counting.  The  lowest  level  of  detection 

*8 

f.,r  radiostrontium  by  this  method  is  U  x  10  uc/ml  and  for  radlobarlum, 
l0'7  nc/ml.  A  precision  cf  about  1C#  cay  be  obtained. 

The  calculation  of  the  total  amount  of  radioactivity  due  to  either 
jailostrontiua  or  radlobariuz,  using  an  internal  proportional  counter.  Is 
t„fc.ie  by  use  of  the  f  cl loving  equation: 


uuc  _  nut  cpta  f 

nr  * 


Vijare 


A  ■■  efficiency  factor, 

B  *  chemical  yield, 

C  «  aelf-acsorptton  factor,  and 

D  *  volume  In  liters. 


<iiila  formula  also  holds  true  for  determination  of  the  other  radionuclide* 

(26) 

A  variation  of  the  preceding  technique  vae  introduced  by  Koolv  1 
A/i.er  the  separation  of  strontium  by  carbonate  and  nitrate  precipitations, 
t, 6 plum  reaioval  is  effected  by  barium  chloride  precipitation  in  a  hydrochloric 
i Id-ether  system.  Any  lantha.num-140  vhich  may  be  present  fj-om  barlura-l^O 
(3r;i:ay  io  removed  by  a  ferric  hydroxide  scavenge.  Sensitivity  in  the  10*10uc/ral 
,.»ngo  is  obtainable  by  this  procedure. 

(f>)  Stront lur.-r*  by  Solvent  Extraction  of  Yttrium-  9 i27) 

Strontium  carrier  is  added  and  carbonate  precipitation  is  performed. 
j-,r.  carbonate  precipitate  Is  dissolved  in  hydrochloric  acid  and  the.  solution 
(avenged  with  hydroxide,  strontium  io  again  precipitated  as  the  carbonate. 
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vhich  18  then  allowed  to  stand,  thus  permitting  yttrium-90  ingrowth,  n* 

yttrlua  la  extracted  with  tributyl  phosphate  equilibrated  vlth  lb»  nitric 

«td.  Storing  extraction,  the  yttrlua  Is  washed  from  the  organic  phase 

vlth  0.^  N  dried  la  a  plnnchet,  and  counted  In  an  anti -coincidence 

beta  counter. 

A  deeontiu'lr.afcion  factor  of  approxlmately  ^  la  obtalne<i>  ^ 

chemical  yield  is  greater  than  70 f  for  strontium;  the  recovery  of  yttrlua 
by  the  solvent  extraction  is  85^, 

Strontium-90  determinations,  cuch  as  that  shown  above,  are  normally 
nade  by  measuring  the  activity  of  its  beta  decay  daughter,  yttrium-90.  After 
a  suitable  period  ^  elapsed,  the  lnerowth  Is  separated  from  the 

bulk  of  the  strontium  by  some  suitable  means.  Ohe  yttrluo-90  is  then  purified 

and  counted.  Zhe  strontium  may  then  be  determined  by  calculation  using  the 
following  fomula:^2®^ 

Strontlum-90,  uuc/i  •  -  net  cpm _  , 

(aT(b)(c)(D)Tb>[F;(2.22) 

•  efficiency  factor, 

•  per  cent  extraction  factor, 

•  per  cent  Ingrowth  factor, 

•  chemical  yield, 

•  sample  volume  (liters),  and 

•  decoy  factor;  calculated  from 


where: 

A 

fl 

C 

0 

E 

F 


vhero: 


activity  remaining  after  a  tine  interval  t, 
activity  of  cample  at  core  original  time, 

_ CJ21 _ - _ -  ,  and 

half -life  of  yttrlun-SW  (c2.“*  hrs) 


t  •  time  (hours)  fron  separation  to  counting. 

3y  knowing  the  total  ot rontium-90  activity  ar.d  the  total  activity  dee 
to  radiostrontium,  the  strontlum-89  concentration  may  be  obtained  from  the 
difference. 

(OQt 

(3)  Strontium-90  by  Direct  Precipitation  of  Vttrlum-90v 

Strontium-90  activity  may  atso  be  easily  measured  by  a  variation  of 
the  preceding  method.  After  a  suitable  time  has  elapsed  Tor  yttrlum-90  in¬ 
growth,  yttrium  carrier  is  added  and  subsequently  precipitated  as  the  hydroxide 
by  the  addition  of  codlura  hydroxide.  The  radloyttrium  la  then  purified  by 
standard  radiochemical  techniques,  which  will  depend  on  the  concentration  and 
species  of  the  other  rare  earth  radionuclides  present.  Resolution  of  the 
decay  curve  will  be  necessary  if  yttrium-91  is  present. 

(U)  Radlor.trontlum  by  Ion-Exchange  Methods 

Kahn  and  Reynolds^0  ^  utilized  ion-exchange  resins  for  the  concentra¬ 
tion  of  radiostrontlum.  Strontium  carrier  is  added  to  a  10-liter  water  sample 
and  the  oolution  is  pasced  upflov  through  a  cation  exchange  resin  in  hydrogen 
form.  The  strontium  is  then  eluted  with  l^K  nitric  acid.  The  radiostrontlum 
Is  then  separated  from  other  alkaline  earth3  and  purified  by  previously 
mentioned  methods, (2S)(27)(29)  Jepondlng  on  vhether  total  radiostrontium  con- 

centratlon  i3  desired  or  strontium-90  alone. 

By  means  of  concentration  of  activity  with  lor.-cxchonge  resins,  great 

-10  /  . 

sensitivity  can  be  obtained,  possibly  down  to  the  10  uc/ml  range. 
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A  similar  iou -ex  change  method  vas  described  by  Kahn,  Eastwood,  and 
Using  appropriate  radiochemical  techniques,  other  long-lived  radio¬ 
nuclides  could  also  be  determined-  In  every  case,  the  recovery  of  activity 
was  greater  than  $9jf. 

Bryant,  Sattiznhn  rnd  Warren^  ^  utilized  an  ion-exchtngc  procedure 
to  separate  yttrium  from  strontium  after  a  ten  day  ingrowth  period.  The 
yttri'#1  was  then  selectively  el^ced  from  a  cation  exchange  resin,  followed  by 
a  rca'isorptlon  on  another  resin  and  finally  counted.  The  radiochemical  re¬ 
covery  was  greater  than  97$*  Gravimetric  measurement  is  not  required  when 
using  this  procedure. 

(5)  Radlostrontlum  by  Evaporation  Methods 
(  13) 

According  to  Libby'  strontium-90  may  readily  be  determined  in  low 
concentrations  by  concentrating  the  activity  by  evaporation. 

The  sample  (l  liter)  is  evaporated  to  dryness,  and  dry  strontium 
nitrnl?  carrier  is  added.  Following  addition  of  calcium  chloride,  to  the  dry 
residual  the  sample  is  dissolved  in  phosphoric  acid.  Yttrium  ia  then  removed 
by  the  addition  of  lanthanum  carrier,  followed  by  hydroxide  scavenging.  The 
precipitate  is  redlseolved  and  the  milking  repeated  twice  more.  The  third 
precipitate  is  saved,  Ignited  and  counted.  The  analysis  itself  takes  two 
hours  td  perform.  It  is  necessary  to  utilize  decoy  curves  to  determine  the 
yttrlufii»9C  activity,  due  to  the  presence  of  yttrium-91  and  other  rare  earths. 

h.  Radloceslua 

Padiocesium  i3  considered  to  be  one  of  the  more  important  of  the  bio¬ 
logically  hazardous  radionuclides.  Cesium  is  chemically  similar  to  the  other 
alkali  petals,  such  as  sodium  and  potassium,  which  are  commonly  found  in  the 
internal  organs.  Thus  ionic  interchange  within  the  body  is  easily  accomplished. 
Co-137  “!>d  Cs-ljU  with  half-lives  of  33  year3  and  2.2  year3  respectively 
constitute  the  greatest  radiocesium  danger. 


-tit. 

Selective  alkali  precipitation  and  ion -exchange  elution  have  proven 
to  he  the  test  net  hod  a  of  separation  for  the  analysis  of  radiocesium.  Deter¬ 
mination  of  the  various  isotopes  appears  to  be  most  feasible  by  gamma  spectro- 
ftcanr- 

(1)  Radloceslura  by  Fhosphomolybdate  Precipitation^  ^ 

After  the  addition  of  cesium  carrier  to  the  sample,  phosphoric  acid 
and  ammonium  molybdate  are  added  and  the  cesium  is  precipitated  os  cesium 
ammonium  phosphonolybdate .  After  dissolving  the  precipitate  in  sodium 
hydroxide,  the  cesium  is  reprecipitated  as  cesium  eobaltinitrite.  This 
precipitate  is  dried  at  100°C  on  an  aluminum  planchet  and  counted  .in  an  in¬ 
ternal  proportional  counter.  A  93$  removal  of  cesium  activity  was  attained 
by  the  phosphomolytdate  precipitation.  Concentrations  as  high  as  10 pc/ml 
and  as  low  as  10”1®  uc/ml  verc  measured  with  good  accuracy.  The  too  principal 
cesium  radioisotopes  measured  by  this  method  vere  cesium-137  and  cesium-13b. 

(2)  Radiocesium  by  Cobaltlnitrlte  Precipitation 

After  the  addition  of  cesium  carrier,  the  cesium  and  other  Group  I 
cations  vere  precipitated  as  the  eobaltinitrite  after  the  addition  df  sodium 
nitrite  and  cobaltous  chloride.  After  vashlng,  the  eobaltinitrite  is  dis¬ 
solved  in  hydrochloric  acid.  Cesium  siiicotungstatc  is  then  precipitated  by 
the  addition  of  silicotungstlc  acid.  The  silicotungstate  is  then  dissolved 
in  dilute  base  and  the  solution  scavenged  vith  ferric  hydroxide.  The  excess 
tungsten  is  removed  as  the  insoluble  trioxide.  Cesium  and  sodium  perchlorates 
are  then  precipitated  by  the  addition  of  perchloric  acid  and  absolute  alcohol 
to  the  solution.  The  sodium  perchlorate  is  then  removed  by  washing  the  pre¬ 
cipitate  vith  absolute  alcohol.  The  final  cesium  precipitate  is  then  vashed, 
dried  and  counted.  This  method  has  been  proven  to  be  quite  accurate  and  is 
nov  in  vide  use. 
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Osmond,  et  ad,™  co-precipltated  cesium  cobaltinltrtte  on  potassium 
cobaltlnitrite  and  followed  vlth  a  similar  purification  procedure.  Thle 
method  Las  found  little  popular  -J«. 

(3)  Padloceslua  V.~  rillcotungst nlc  Precipitation^^ 

This  procedure  is  booed  re  a  method  originally  described  by  Yienagata 
and  Yamigata^"^ .  The  cesium  Is  first  precipitnted  as  the  sllicotur.gi  tnte. 
Dissolution  of  the  precipitate  is  then  followed  by  n  ferric  hydroxide 
scavenging.  The  cesium  is  then  ^precipitated  as  the  dlplcrylunlnnte.  This 
salt  Is  dissolved  In  4-!Bethyl-2-ps=.tanone  nnd  the  cesium  Is  extracted  by  moons 
of  2M  hydrochloric  acid.  The  cerium  13  finally  precipitated  ns  the  perchlorate, 
In  which  form  It  is  dried,  velghei,  and  counted.  The  chemical  yield  lo  0 Of>,  and 
eight  analyses  may  be  performed  —  eight  hours. 

(23) 

(4)  padloceslum  by  s  Co-crystnlllrntlon  Procedure 

A  method  has  been  describe!  by  which  cesium  13  separated  from  the  bulk 
of  the  alkali  elements  and  the  miasd  fission  products  by  eo-erystallUntion 
with  amrsonlum  aluminum  sulfate.  The  ammonium  salts  are  decomposed  by  heating, 
end  the  cesium  is  precipitated  for  counting  from  a  dilute  hydrochloric  ncld 
solution  as  the  chloroplntlnate.  This  method  has  been  in  use  as  standard 
procedure  In  son*  laboratories. 

(5)  Rndloceslum  by  I.-c.-Hxchangc  Methods 

*  IT  1 

Kahn,  Eastwood,  and  Lacy'  have  developed  n  separation  scheme  lor 
the  analysis  of  the  more  hazardous  radionuclides  by  ion  exchange.  The  npeclflc 
mdloelencnts  studied  were  cesium,  cerium,  cobalt  and  strontium.  The  lower 
limit  of  detection  waa  decrease!  a  hundredfold  by  concentrating  the  radio¬ 
activity  of  n  large  sample  with  s  cmtlon  exchange  resin. 

Cesium  was  selectively  el-ted  by  means  of  CH  hydrochloric  acid  1  rom 
the  cation  exchnnge  resins  stuuiei-  Standard  purification  procedures  were 


then  U6ed  to  prepare  the  sample  for  counting.  The  cesium  vas  separated  from 
the  other  alkali  cation  contaminants  by  precipitating  It  as  the  olllcotu-.g- 
state,  dissolving  the  precipitate  In  sodium  hydroxide,  and  then  repreclpltattng 
It  as  the  perchlorate  for  counting.  The  activity  recovered  for  all  radio¬ 
nuclides  vas  greater  than  9 9*2^. 

( 37) 

Tsubota  and  Kitano  found  that  an  ammonium  formate  -  formic  acid 
buffer  of  pH  3-2  selectively  eluted  the  alkali  metals  from  a  cation  exchange 
resin.  This  buffer  vas  employed  particularly  for  the  determination  of  radio- 
cesium,  being  somevhat  superior  to  hydrochloric  acid  or  citrate  buffer  as 
an  eluant. 


C .  Radio iodine 

Radiolodlne  may  be  either  beta  or  gamma  counted,  although  gamma 
scintillation  Is  preferred  for  determination  of  the  Individual  rndlolsotopea 
of  Iodine. 

(l)  Radlolodlr.c  by  Chemical  Methods 

(llendenln  and  Metcalf  determined  radiolodlne  activity  by  on 

extraction  jXiriflcatlon  procedure.  Carrier  sodium  iodide  Is  added  to  the  vater 

sample,  and  Interchange  io  accomplished  by  oxidation  to  the  lodato  vlth  sodium 

hypochlorite  In  basic  colutlon,  folloved  by  reduction  to  the  iodide  by  codlura 

bloulflto  In  acid  solution.  Sodium  nitrite  Is  then  added  to  oxidize  the  Iodide 

to  elemental  Iodine  vhich  is  then  extracted  into  carbon  tetrachloride.  The 

Iodine  Is  further  purified  and  concentrated  by  back-extraction  into  sodium 

bisulfite  solution  vhich  la  finally  gamma  counted  at  the  lodlnc-131  photo peak. 

•8 

The  sensitivity  of  the  method  Is  approximately  10*  c/ml  for  a  liter  sample. 

The  problem  of  Incomplete  carrier  Interchange  Is  overcome  by  this 
procedure  vhich  employs  oxidation-reduction. 
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The  same  authors  also  describe  &  procedure  in  which  the  iodate  Is 
reduced  directly  to  the  iodide  by  hydroxylaalne  hydrochloride.  The  iodide 
obtained  after  the  sodium  blsunflte  bach-extraction  Is  precipitated  vitta 
sliver  nitrate  as  silver  iodide  vhlch  is  then  dried,  veiled,  and  beta 
counted . 

A  similar  but  faster  nethod  for  radiolodlne  determination  va3  de- 
veloped  by  Levis  .  He  developed  a  continuous  extractor  employing  tvo 
centrifugal  pumps,  to  promote  mixing  ar.d  extraction.  The  iodine  was  extracted 
into  carbon  tetrachloride  and  backcxtraeted  into  bisulfite  solution  In  one 
operation.  TRie  gaaaa  emission  of  a  pipetted  sample  vas  then  measured.  The 
time  of  analysis  vae  cut  from  2  hours  to  30  minutes.  A r.  accuracy  of 
93.0  -  d.of  va a  achieved. 

d.  Total  Radio  Rare  Earths  Determination 

(1)  Total  Rare  Earth  Activity  by  Fluoride  Precipitation 

Home  and  Martens  developed  a  rapid  method  for  the  determina¬ 

tion  of  rare  earth  activity  by  which  the  beta  and  gamma  activities  arc  deter¬ 
mined  separately.  The  method  la  quite  rapid,  the  gamma  determination  taking 
only  a  half-hour. 

Tvo  rare  earth  fluoride  precipitations  serve  to  remove  zirconium  and 
niobium,  the  principle  gamma  emitters  of  fission  material.  Barium  and  stron¬ 
tium,  vhich  arc  heavily  eopreclpitatcd,  are  removed  by  hydroxide  precipitation 
after  the  addition  of  holdback  carrier.  Tne  cample  may  then  be  gamma  counted. 
To  estimate  tho  beta  activity,  the  rare  earths  must  be  precipitated  as  the 
oxalates  before  counting. 

(2)  Total  Hare  Earth  Activity  by  an  Alternate  Mjthod 

(1*1) 

ocldrldge  and  Hume  developed  a  method  similar  to  that  above 
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fcr  the  det*  —  ir.atlon  of  total  rare  earth  activity.  It  Is  somewhat  lengthier, 

Vat  zcre  arr.rate,  particularly  far  beta  counting. 

-'able  c-rium  Is  added  as  a  carrier  for  the  entire  group  of  rare 

earths.  !2at  rare  earths  are  then  precipitated  aa  fluorides,  rediasolved  in 

boric  and  -_i*.ric  acids,  and  repreelpltatcd  as  the  hydroxides.  Zirconium  la 

then  recover  by  precipitation  aa  the  lodate  after  reduction  of  all  the  rare 

eartha  to  '-it  trlvalent  6tate.  Following  hydroxide  scavenging  for  alkaline 

earth  remrrsl,  the  rare  earths  arc  precipitated  aa  the  oxalates  and  beta- 

counted  It  an  ir.tem&l  proportional  counter. 

<U2) 

Alctad  and  Pappas*  '  employed  a  similar  technique,  using  lanthanum 
as  a  carrier  for  all  the  rare  eartha.  Following  separat'on  and  purification 
of  the  gro-t  ;f  rare  earths,  individual  separation  vaa  attained  by  elution 
from  a  catiot  exchange  column  with  omonlura  lactate.  The  ion-exchange  separa¬ 
tion  step  vti  performed  without  additional  carrier,  and  the  recovery  of 
activity  ftc  each  radionuclide  was  assumed  to  be  the  same  as  that  for  the 
lanthanum-  Thr:  advantage  of  the  non-carrier  separation  Is  that  a  "weightless** 
sample  is  rd'air.ed,  thus  eliminating  cclf-aboorptlcn  and  self -scattering 
factors. 

e.  Kir—llancouq  Radionuclides 
(ij  Fad  lobar  jura 

>  termination  of  -adlobariun  r.ay  easily  be  made  during  the  radlo- 
8tro,ntlum  ar-vlyrls.  toth  t;ron»lum  and  barium  form  Inooluble  carbonates,  and 
subsequent  ::-*ration  may  bo  made  by  precipitation  of  barium  as  the  chromate. 

This  cethoi.  by  :->hn  and  Straub,  was  doccribcd  earlier  under  radio  strontium. 

(1*3) 

Hunter  and  r-.-oina  followed  the  tame  procedure  except  that  the  precipitate 
Is  counted  *j  tie  carbcnatc  Instead  of  the  chloride. 
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Mlnkklnen^  has  developed  a  procedure  for  the  analysis  of  radio* 
barium  In  vhlch  barium  chloride  carrier  la  specifically  precipitated  ns  the 
nonohydrate  from  a  concentrated  hydrochloric  acid-ethyl  ether  mixture.  This 
precipitate  Is  dissolved  In  water  to  vhlch  ferrto  hydroxide  io  subsequently 
added  at  a  scavenger.  The  barium  16  then  Isolated  ag  the  chromate.  The 
chenlcal  yield  of  barium  chromate  Is  about  70$,  Tne  barium  chroouite  precipi¬ 
tate  la  then  set  aside  for  13**  hours  to  nllov  the  barium-140  and  Its  daughter 
lanthar.um-lUo  to  come  to  equilibrium. 

(2)  RqJiocerlua 

( 143 

Borgus  and  Engelkemelr  '  have  developed  method  of  analysis  for 
radlocoi  1  uni  (cerlura-144),  employing  a  series  of  precipitations.  After  the 
addition  of  cerium  carrier  and  alkaline  earth  holdback  carriers,  cerium  was 
separated  along  vlth  the  rare  earths  by  fluorid«  p. ocipitatlon.  Alternate 
hydroxide  und  fluoride  precipitations  are  then  employed  to  remove  zirconium 
and  alkaline  earths  completely.  The  cerium  Is  then  oxidized  to  the  tetrnvnlent 
state  vlth  perchloric  acid  and  precipitated  as  the  lodate.  Thorium  activity  la 
then  removed  as  the  lodate  following  the  reduction  of  cerium  vlth  sulfur 
dioxide.  After  scavenging  vlth  thorium  end  further  rcpreclpltutlon,  the  cerium 
lo  finally  precipitated  as  cerouo  oxalate,  vhlch  Is  subsequently  Ignited  to 
the  oxide.  After  velghlng,  the  cerium  activity  io  courted  on  on  Internal 
proportional  counter. 

Tho  percentage  of  chemical  recovery  Is  only  about  7/0$,  but  the  radio* 
active  purity  of  the  final  proclpltnto  is  very  high. 

It  lo  likely  that  this  process  could  be  shortened  bccauce  of  the 

probable  absence  of  thorium  from  fallout  contaminated  vatcr. 

(40 

Ames  employed  a  similar  procedure  for  the  analysts  of  total 
radiocerium.  The  concentration  of  cerium-144  Is  obtained  by  counting  la- 
mediately,  using  a  217  mg/cm  aluminum  absorber.  This  shields  all  the 
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cerlua  betas,  counting  only  those  t rom  prascodymium-lUt,  the  daughter  product 
of  cerlua-lbU.  The  difference  between  the  total  radloceriua  concentration  and 
cerium-lW*  concentration  la  assum'd  due  to  cerlum-ll*l,  a  neutron  induction 
product. 

i  1  *  x 

Barnes'  7  describes  a  mel'iod  by  which  radiocerium  Is  icolated  by  an 
extraction  method  and  counted  for  0*riua-l4W.  Carrier  interchange  is  accom¬ 
plished  by  means  of  an  oxidation-reduction  cycle.  Two  dibutyl  phosphate  ex¬ 
tractions  are  then  carried  out  to  remove  the  heavy  elements.  The  cerlua  la 
then  precipitated  as  the  fluoride^  oxidized  to  the  tetravalent  state,  and  ex¬ 
tracted  into  hexonc.  Finally  it  id  converted  to  the  dioxide  and  counted  for 
cerium- lMr  by  means  of  u  beta  counter  using  a  217  mg/c»  aluminum  absorber, 
ttie  chemical  yield  varies  from  “pO* 

Kahn  and  Reynolds^  ^  have  developed  another  method  for  the  determina¬ 
tion  of  radiocerium  activity  employing  the  use  of  ion-exchange  reslno.  Con¬ 
centration  of  rudlocerlut  and  added  carrier  was  attained  by  passing  a  ono-llter 
water  sample  through  either  Dowex^t-O  on  the  hydrogen  cycle  or  1R-220  on  the 
sodium  cycle.  Cerium  was  oclcctively  eluted  from  the  reoln  with  25  ml  of 
3M  nitric  acid.  Purification  of  tit"  rndloccrium  is  obtained  by  successive 
precipitations  as  the  fluoride,  th<>  hydroxide  and  finally  as  the  oxalate.  The 

cerium  oxalate  is  weighed  to  dete»tnlr>*  carrier  loss  and  then  counted. 
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Sensitivity  of  detection  was  incrcdotd  to  approximately  10  pc/ml. 

(3)  Radlqruthcnlua 

Glendenli/1*^  has  developed  f>n  analysis  scheme  for  the  determination 
of  radloruthenlum  which  hoo  been  utilised  us  a  standard  method  of  procedure. 
Separation  of  the  radloruthenlum  vl<b  carrier  ruthenium  is  accomplished  by 
oxidation  with  perchloric  acid  to  volatile  tetrooxide  which  is  subsequent¬ 
ly  distilled  off.  Sodium  bismuthnl0  *6  added  prior  to  distillation  to  prevent 


volatilization  of  the  halides  by  oxidizing  then  to  the  oxyaclds.  The  ruthenium 
is  absorbed  into  a  sodiun  hydroxide  solution  and  precipitated  in  the  forn  of  its 
lever  oxides  toy  reduction  with  ethanol.  The  ruthenlun  oxides  are  then  dis¬ 
solved  in  hydrochloric  add,  and  rathe n iua  la  precipitated  in  the  metallic 
ctate  by  reduction  with  magnesium  metal.  The  ruthenium  metal  la  weighed  for 
carrier  loss  and  counted.  The  chemical  yield  is  about  65$. 

Melnlck^®^  ,  using  an  almost  identical  chemical  procedure,  describes 
a  process  for  the  determination  of  ruthenium-106 .  By  using  a  series  of  heavy 
aluminum  absorbers  with  a  beta  proportional  counter,  a  correction  factor  for 
the  presence  of  ruthenium-103  (half -life  -  U2d)  may  be  obtained  by  extrapola¬ 
tion.  Extrapolation  for  the  correction  la  necessary,  in  that  while  counting 

a 

the  betas  from  ruthenium- 106  through  a  210  mg/cmf  aluminum  absorber,  the  gansna 
rays  from  ruthenium- 10 3  are  also  counted.  If  the  activity  due  to  ruthenium- 103 
alone  la  desired,  gcuma  scintillation  may  be  used  with  a  2,000  mg/cm  aluminum 
absorber. 

Merritt^)  describes  a  method  or  radioruthenium  analysis  in  which  the 
ruthenium  Is  determined  In  the  presence  of  strontium,  cesium  and  cerium  by 
fuoion  and  extraction.  Following  the  addition  of  carrier  ruthenium  to  the 
water  sample,  the  sample  is  fused  with  potassium  hydroxide,  sodium  nitrite, 
and  sodium  carbonate  at  550°C  for  two  hours.  The  melt  is  leached  twice 
with  water  to  dissolve  the  ruthenato  and  cesium,  leaving  strontium  and  curium 
In  the  residue.  Strontium,  cerlim  and  cesium  are  then  analyzed  by  previously 
mentioned  procedures.  The  ruthenium  is  extracted  from  the  leach,  after  the 
addition  of  periodate,  with  carbon  tetrachloride.  The  ruthenium  Is  then 
stripped  from  tho  carbon  tetrachloride  with  6t!  hydrochloric  odd.  Magnesium 
Is  then  used  to  reduce  the  ruthenium  to  Its  elemental  metallic  state  in  which 
form  it  Is  subsequently  counted.  This  method  has  proven  satisfactory  for 
soli  leaching  analysis,  etc.,  but  could  veil  be  used  on  water  samples. 
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( *»)  Rad  lo  z  i  rcon  lua -Rad  ton  icb  lea 

Steinberg^®)  has  developed  &  method  for  the  determination  of  zirconium- 
niobium  activity  which  has  been  tested  ar.d  proven.  Following  the  addition  of 
both  zirconium  and  n  lob  lira  carriers,  the  two  elements  are  completed  by  the  ad¬ 
dition  of  oxalic  acid,  thorium  cxalate  scavenging  is  carried  out  by  the  ad¬ 
dition  of  potassium  chlorate  to  the  acidified  solution.  The  niobium  oxide 
is  redlssolved  and  reprecipitated  as  the  cxlde  hydrate.  The  water  is  then 
driven  off  by  ignition  and  the  nloblun  oxide  is  weighed  and  mounted  for  count¬ 
ing.  the  zirconium  Is  precipitated  from  the  original  solution  as  the  phen- 
ylarsonatc  by  the  addition  of  phenylarsonle  acid.  Thla  precipitate  is  dlosolved 
In  oxalic  acid,  forming  the  soluble  oxalate,  the  zirconium  la  reprecipitated 
as  the  phenylarsonate,  which  Is  subsequently  ignited  to  zirconium  dioxide, 
tills  precipitate  Is  weighed  and  counted.  3oth  precipitates  (niobium  and 
zirconium)  are  beta  counted  on  an  Internal  proportional  counter. 

Brady  and  EngelKcmcir^^  have  developed  a  phosphate  method  for  the 
determination  of  activity  due  to  zirconium -niobium  which  Is  somewhat  lengthier 
than  the  preceding  method.  Separation  of  zirconium  1b  based  upon  the  precipi¬ 
tation  of  zirconium  phosphate.  The  phosphate  Is  then  redlssolved  as  the 
fluoride  and  precipitated  as  the  phenylarsonate  which  Is  subsequently  Ignited 
for  counting.  Ilioblua  is  precipitated  la  acid  solution,  dissolved  a6  the 
fluoride  complex  and  finally  weighed  at.l  counted  as  the  oxide.  About  10  houra 


lc  required  for  the  entire  analysis.  The  chemical  recovery  of  zirconium  la 
about  75&  and  that  of  niobium  io  about  S0$. 


A  method  Is  alro  described  by  Stanley'"'  '  for  determination  of 


zlrconlum-95.  Carrier  Interchange  io  effected  by  the  formation  of  the  fluoro- 

.2 

zlrconate  complex,  ZrF^  .  Bare  earth  activities  are  removed  by  lanthanum 
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fluoride  scavenging.  The  zirconium  lo  then  separated  by  means  of  three  barium 
fluorozirconate  precipitations.  The  zirconium  is  finally  precipitated  vtth 
mantle  lie  add  from  hydrochloric  acid  solution  and  ignited  to  the  dioxide.  In 
which  fort  It  Is  weighed  and  counted.  A  chemical  yield  of  about  755&  io 
obtained.  It  is  essential  to  not  begin  the  analysis  until  zirconium-97 
{tjy2  ■  t?  hours)  has  had  time  to  decay  to  a  negligible  amount.  After  the 
final  precipitate  has  been  obtained,  counting  should  begin  Immediately  so 
that  no  appreciable  nloblum-95  has  time  to  grow  in.  A  beta  proportional 
counter  should  be  employed  for  the  counting. 

3.  3wiuj«i/  of  Radionuclide  Analyses 

For  the  convenience  of  the  reader,  additional  references,  as  well 
as  those  previously  cited,  are  summarized  In  Table  XVI, 


1.  Radiostrontium 


25,  26,  27,  29, 
30,  31,  32,  33 


1*9,  *>3,  54,  55, 
56,  57,  6l>  62 


2.  Radiocesium 


23,  31,  34,  35, 
36,  37 


29,  *i9,  54,  55, 

50,  r4,  65 


3.  Radioiodiae 


28,  38,  39 


55*  W 


1*.  Total  Radio  Rare 
Earth  a 


4o,  4l,  42 


55*  ^ 


5.  Radlobarlum 


24,  25,  43 


6.  Radioccrium 


30,  44,  45,  46 


29*  4»  ^9,  63 


Radioruthenium 


47,  48,  49 


8.  Radiozirconium 
Radioniobium 


50,  51,  52 


•  l?2- 


3.  Instrumentation 

Quantitative  analysis  of  radioactive  elements  in  trace  concen* 

♦  ration*  la  done  most  conveniently  by  meens  of  r&diat  ion -detecting  instru¬ 
ments.  .  Due  to  the  very  small  amounts  °f  radionuclides  present,  sueh 
Instruments  are  usually  superior  in  npei  and  sensitivity  to  conventional 
vet  chemical  analyses.  Hovcver,  wheil  complete  chemical  analyses  are 
required,  extensive  use  is  made  of  vet  chemical  procedures  for  the  separa¬ 
tion  of  various  isotopes,  as  descrited  °n  pages  105  -  121,  Such  separations 
are  not  required  vhen  only  gross  fission  product  activity  is  to  be  mea¬ 
sured  or  when  only  one  radioisotope  it*  present. 

If  the  activity  of  a  sample  18  f-ficiently  high,  this  activity 
ray  be  determined  directly  using  a  simple  detector  and  ratemeter.  At 
very  low  concentrations  in  aqueous  efeHipl*6  or  for  beta  emitters,  for  vhich 
self-absorption  in  the  cample  is  a  tiihjor  factor,  preconcentration  of  the 
cample  and  specially  designed  instruments  tuch  as  a  flow  counter  may  be 
required. 

A  survey  of  commercially  aV§lls!l*  portable  radiation  Instrument* 
hue  been  carried  out  and  the  information  obtained  is  summarized  in 
Tables  XI  ^  and  XjJC.  These  are  all  illPtruser.ts  of  the  ratemeter  type 
capable  of  recording  count  rates  dowll  t®  a  few  counts  per  minute.  To 
establish  their  usefulness  for  the  determination  of  fallout  product 
concentrations  in  water  it  is  necessary  to  correlate  the  sensitivity  of 
the  equipment  with  different  types  of  detectors  and  reasonable  sample 
volumes  with  the  level  of  contamination  concentration  expected  in  emergencies. 

For  rav,  untreated  uqueous  polutiona,  Figure  15  shows  the  relation 


of  nuclide  concentration  to  measured  d®&«  rate  for  various  isotopes  of 
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Interest  for  a  representative  eonsrorciHl  instrument.  It  ha?  toen  shown  ir. 
Tables  and  X'  that  fallout  product  cencshf rations  of  the  order  of 

10  5  •  |C  tic/r.I  r.i/  ce  expected  In  the  water  supplies  of  typical  watersheds 


durirg  the  early  path -attack 
'£*y 

Xshr»’  rhovs  *.i*  relative 
survey  meters  vita  a  lo.v-st- 


period.  Table  M’..  token  from  l*»cy  and 
detectability  of  selected  radioisotopes  vith 
scale  sensitivity  down  to  0.01  mr/hr  assuming 


ideal  counter  gi-tmetry. 

If  it  'J  merely  desired  to  determine  the  presence  or  absence  of 
these  feign- level  fallout  coacent  rations.,  It  is  evident  that  available 
survey  inatroett.*  sooul'd  be  capable  of  supplying  this  information,  pro¬ 
vided  that  they  are  Iftted  with  &  detector  head  suitable  either  for  immersion 
into  the  liquid  rumple  or  for  flow-through  measurement.,? . 


‘fr/vever,  in  sort  cases  the  concentrations  rw wared  ore  veil 
above  the  ICKF  it-stlmx:  permisaible  levels  of  concentration  in  voter-  ' 1  or 
close  tc  the  limit;  of  detectability  or  the  instruments  listed  in  Table  XV  a  1 1 . 
For  this  reason  it  seems  necessary  to  develop  end  make  available  procedures 
end  facilities  for  the  sore  accurate  meas'.urement  of  lav -level  fission  product 
concentrations  to  ensure  the  safety  of  available  sources  of  drinking  water* 
Thic  eevelopeeat  car.  follow  twe  paths: 

1.  f-eeccccctration  of  the  cample  to  increase  the  number  of 

radioactive  shea?  to  a  icvel  that  is  readily  detected  by  conventional 

•• 

survey  instruments;  or 

2.  Increas ic ;  the  eensitivity  cf  the  detectors  by  an  increase 
in  r.ample  volume,  cackgrcmd  reduction  by  edded  shielding  and  coincidence 
circuit  Arrangers*  tie,  or  more  intimate  contsct  between  cample  and  detector 
volumes,  t?  in  projor: lonal  flow  counters  or  .’n  liquid  scintillation 
counting. 


Best  Available  Copy 


Table  XVII 

Approximate  Minimum  t-f-tec table  Concent ;~atlor.s  In  Water  (uc/nl) 


Tor  Survey  Meters  vlth  Pull  Scale  Sensitivity 


idlolrotope 

0.01  mr/hr 

0.1  nr/hr 

1.0  rcr/hr 

10  mr/hr 

_  106 

Ru 

5  x  10*6 

5  X  10*5 

.1 

5  x  10 

5  x  10*- 

Pa1**0 

1  x  10*7 

1  x  10*" 

5  x  10  4 

5  x  10*J 

P32 

5  x  10*6 

1  X  10*U 

1  X  10' 3 

1  x  10*'' 

,20i» 

-i, 

-3 

,  -2 

-2 

T1 

5  X  10 

1  x  10  J 

1  x  10 

5  x  10 

Zr95 

1  X  10 

1  x  10*3 

5  x  10* 3 

5  x  JO*2 

I131 

5  x  10*4 

1  x  IQ*3 

.2 

1  x  10 

5  x  10*2 

Ka24 

1  x  10*5 

.U 

1  x  10 

1  x  10*3 

1  x  lo*2 

K1*2 

5  x  10*6 

5  x  10”5 

5  x  lo*4 

5  X  10*3 

Sc**6 

1  x  10*3 

5  x  10*3 

5  X  10*2 

1  x  10*1 

Cr'1 

5  x  10*4 

1  x  10*3 

1  x  10*2 

5  x  10*2 

CO60 

1  x  1 O*3 

5  x  10*3 

5  X  10*2 

1  x  10*1 

Sr90 

1  x  10’5 

_k 

l  x  10 

5  X  10*u 

5  x  10*^ 

Au193 

1  x  10*5 

l  x  10 

5  x  10'4 

5  x  10* 3 

c*137 

1  x  10*3 

5  x  10’3 

5  x  10*2 

1  x  10*1 

Values  extrapolated  from  Mcy,  W.  J.  anl  Kahn,  P.,  "Survey  Meters 
end  Electroscopes  for  Monitoring  Radioactivity  In  Wntcr,"  J. A.W -W.A. , 
*•6,  59  (1951*) 
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Preconcert  ration  Methods 

Many  method a  iucooucrui.tai.joa  have  disadvantages  which  limit 

their  use  in  on  er-ergeocy.  Concentration  by  evaporation,  for  example, 
requires  a  great  deal  of  care  and  skill  to  o'oialu  reasonably  accurate 
results.  The  use  of  Ion  exchange  resins  to  rerove  the  activity  from  the 
sample  can  give  very  good  rc-pt'oduceability  and  In  some  applications  can 

/ /fQ  \ 

be  alnost  completely  automated.  Basons  and  Lauderdale'  '  have  developed 
a  continuous  water  monitor  which  employs  a  long  C.M.  tube  at  the  center  of 
a  column  of  cation  exchange  resin.  By  employing  two  columns,  one  may  be 
used  while  the  other  Is  being  regenerated  with  strong  acid.  With  thlG 
system  sensitivities  of  around  10*7  uc/ml  have  been  obtained. 

The  mixed  resin  bed  13  potentially  capable  of  sensitivities  In 
the  10*7  -  10*®  uc/ml  range  when  used  with  law  background  counters. 

Although  the  measurement  of  radioactivity  on  resin  presents  problems  of 
self  absorption  in  the  cample,  reasonable  efficiencies  and  rcproduccabllitles 
can  be  obtained  if  sufficient  care  io  taken  In  sample  preparation. 

One  method  which  may  be  employed  to  overcome  self -absorption 
in  the  resin  lo  to  detect  the  beta  particle#  with  a  scintillation  resin. 
Little  work  hao  been  done  on  scintillation  resins,  but  when  technique# 
are  developed  this  should  become  a  very  valuable  tool  in  contamination 
evaluation  work. 

A  highly  accurate  way  in  which  the  radioisotope  content  of  water 
can  bo  determined  Is  by  radiochemical  analysis.  By  chemically  separating 
the  various  elements  of  Interest  and  then  measuring  the  uctivity  of  each 
constituent,  the  activity  concentration  as  well  as  the  Identity  of  the 
radio  contaminants  may  be  obtained.  This  type  of  analysis  requires  a 
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skilled  radiochealst  and  such  an  individual  nay  he  diff'lcalt  to  flnd  ln 
emergencies.  Solvent  extraction  methods  have  also  beof»  l,fled>  but  at 
present  a  method  involving  &  mixed-bed  ion-exchange  re?  in  6cee »  to  offer 
the  simplest  practical  approach^^. 

Detector  Development 

Improvement  ln  detector  sensitivity  can  be  aelti,iVcd  by  a  suitable 
choice  of  detector  characteristics,  detector  geometry*  ?l,d  electronic 
circuit  parameters.  In  order  to  measure  the  contamlnal *evel  vater 
supplies  the  radiation  emitted  from  the  water  oample  muef-  fce  determined 
as  efficiently  and  as  accurately  as  possible.  Because  uf  the  statistical 
nature  of  the  radiation  emission,  higher  accuracies  cat)  bc  obtained  only 
by  increasing  the  total  number  of  counts  regiotered  by  t-l*c  detector, 
while  at  the  same  time  keeping  background  counts  to  a  F’or  thi« 

reason,  at  low  levels  of  activity  individual  counting  dlrPuito  or  scalers 
aro  superior  to  integrating  counting  circuits  ouch  as  fal-croetera. 

If  the  radioactivity  is  composed  largely  of  K*l.«  articles,  a 
mnjor  problem  In  measuring  beta  activity  in  voter  is  the  short  distance  which 
a  beta  particle  can  travel  in  this  dense  medium.  It  is;  therefore*  desirable 
to  bring  no  large  a  volume  of  voter  os  possible  into  thb  Immediate  proximity 
of  the  detector. 

One  of  the  most  convenient  monitors  for  this  Jj'IT P°oe  1°  ^he  flow 
through  type  which  allows  the  water  to  come  ln  contact  the  detector  for 

a  brief  period  and  then  be  flushed  out.  The  moot  common  d'-tector  u=ed  ln 
flow-through  monitors  is  the  G.M.  iu'oe^^K  Various  coiH  lg>ra„lor.3  cf  the 
G.M.  monitor  type  have  been  used  to  check  cooling  water  *.lllusnt  at  rhny 
of  the  water-cooled  reactors  in  the  country. 


Lately,  the  geometrical  flexibility  of  plastic  selntllle.tor3 

ha3  increased  the  volume  of  vater  which  can  be  brought  close  to  the  active 

volume  of  the  detector.  Thor?  are  presently  on  the  market,  water  monitors 

which  use  large,  spherical  scintillators  capable  of  measuring  activities  in 

the  range  of  10*^  pc/nl^^.  Another  method  which  displays  much  of  the 

geometrical  advantage  of  liquid  scintillation  but  has  none  of  the  chemical 

(  7c) 

problems  is  the  use  of  scintillation  fibers  .  Although  this  method 

lh 

has  been  used  for  the  most  part  to  measure  C  and  tritium  in  water.  It 

should  yield  good  results  with  fission  products  and  extend  direct  monitoring 

-8 

capabilities  below  the  10"  yc/ml  range. 

Liquid  scintillation  techniques  have  yielded  very  accurate  results, 
high  efficiencies  and  good  sensitivity  in  the  determination  of  radioactivity 
in  water.  This  method  may  prove  quite  valuable  in  the  hands  of  an  experienced 
radiochemist. 

The  limiting  factor  In  the  sensitivity  of  a  detection  device 
Is  the  background  or  "zero  activity"  response  of  the  device.  Both  electrical 
and  physical  approaches  have  boon  pursued  in  the  reduction  of  this  background. 
Physical  methods  include  the  shielding  of  detectors  and  the  reduction  of 
detector  sensitivity  to  the  gorana  radiation  which  can  penetrate  the  shielding. 
Since  the  cross  cectlon  of  plastic  scintillators  is  quite  high  for  beta 
radiation  while  being  fairly  low  for  gamma  radiation,  very  thin  scintillation 
discs  have  been  used  to  produce  good  beta  detection  effioienceo  with  a 
background  of  less  than  one  half  count  per  minute  for  a  one  Inch  diameter 
sample 

The  electrical  method  of  reducing  background  Involves  a  shielding 
detector  ar.d  anli-colncider.ce  circuits.  The  detection  device  is  thus  fade 
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insensitive  to  most  of  the  background  radiation  which  enters  the  detector. 
Most  of  the  anti -coincidence  detectors  presently  on  the  market  reduce  the 
apparent  background  to  a  range  of  0.5  to  2  counts  per  minute. 

The  only  reliable  method  presently  explicable  to  civil  defense 

use  for  censuring  the  amount  of  specific  isotopes  in  water  is  gamma -energy 

analysis.  Most  gemma  analyzers  are  large  ar.d  quite  complicated  to  operute, 

although  son/C  that  are  portable  ar.d  easy  to  use  arc  r.ow  appearing  cn  the 

(74) 

market.  '  Vith  a  small  amount  of  training  civil  defense  personnel  should 
be  able  to  make  good  use  of  these  instruments. 

Summary 

Preparation  for  the  post  attack  monitoring  of  drinking  water 
supplies  for  radioactive  contaminants ,  may  now  be  divided  into  four  areas 
of  required  action.  These  include  the  installation  of  continuous  monitoring 
equipment  at  treatment  plants,  dissemination  of  information  on  the  use  of 
available  instruments  for  water  activity,  insuring  the  availability  of 
battery-per.-ered  equipment  designed  for  water  or  general  liquid  radioactivity 
assay,  and  continuing  the  research  on  equipment  capable  of  meeting  emergency 
requirements. 

A  large  number  of  water  treatment  plants  would  probably  still  be 
operational  after  a  nuclear  attack.  The  danger  of  contamination  from 
fallout  would  require  that  the  water  be  continuously  monitored  for  contamina¬ 
tion  even  at  levels  well  below  the  maximum  permissible  levels.  In  this  way 
alternative  water  sources  could  be  located  or  decontamination  procedures 
could  be  brought  into  operation  as  soon  as  the  possibility  arises  that 
existing  wa*er  supplies  become  unfit  for  consumption.  It  would  be  desirable 
for  this  permanent  equipment  to  be  able  to  determine  the  relative  concentra¬ 
tions  of  the  various  radlocontaalnants. 


wv.v.v.v.--.v.v.v.-v^ 
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Many  commercial  instruments  which  detect  nuclear  radiations  can 
he  set  up  to  taeasure  the  amount  of  radioactivity  In  voter  to  a  degree  of 
accuracy  suitable  for  emergency  ur£ «  This  measurement  would  require  a 
given  procedure  and  a  cct  of  conversion  tables  for  each  typo  of  Instrument 
or  class  of  Instruments. 

Of  the  widely  available  Instruments,  the  most  common  type,  the 
C.M,  survey  meter,  can  bo  adapted  for  emergency  uce  by  either  dipping  the 
probe  Into  the  water  or  holding  t-ho  probe  close  to  the  water  surface,  to 
measure  radlo-contamlnatlon  levels  down  to  10  ^  uc/ral  of  some  of  the 
Isotopes.  Some  scintillation  metcl'P  can  reduce  this  lower  limit  below 
10*6  yc/ml  with  the  absolute  limit  depending  on  background  radiation. 

Most  lor.  chamber  3ur/ey  meters  are  not  as  sensitive  as  the  G.M.  or  sclntll* 
latlon  types,  although  some  of  the  more  sensitive  ones  might  be  used  when 
no  other  equipment  io  available* 

Companies  marketing  nuclear  detection  equipment  were  asked  to 
supply  detailed  specifications  on  their  survey-type  instruments.  A 
summary  of  Information  received  by  January  30,  19^3  la  given  in  Table  XV.  'I. 
The  Instruments  are  listed  In  order  of  decreasing  sensitivity  (m4nlmua 
full  scale  readings). 

Pertinent  specifications  In  the  tablo  Include  the  maximum  full 
scnle  range,  battery  life,  type  of  ludlatlon  detected,  type  of  detector 
used ,  weight,  the  manufacturer's  model  number,  and  unit  cost  of  the 
instrument.  The  ability  of  on  Instrument  to  discriminate  between  various 
types  of  radiation  Is  also  ir.cluJ"!,  ao  veil  as  the  type  of  radiation  to 
which  a  given  Instrument  Is  sens!'.  Ive,  is  indicated  by  an  "X"  In  the 
appropriate  column.  If  an  "0"  appears  under  the  Indicated  type  of  radiation 
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it  simply  means  that  It  can  be  shielded  out  by  a  movable  absorber.  In  the 
event  an  instrument  may  be  powered  by  standard  'D"  cells,  the  number  of 
batteries  required  is  indicated.  Instnasento  which  arc  specified  to  be 
moisture  proof  or  immersion  proof  are  indicated  to  be  "weather  proof"  by 
an  "X"  in  the  appropriate  column.  It  is  believed  that  the  nbove  cited 
characteristics  are  the  tost  Important  to  be  considered  in  selecting  an 
instrument  for  surveying  and  measuring  water  supplies  and  for  other  general 
uses.  The  rrianuf acturer  of  each  instrument  may  be  determined  by  referring 
to  the  "Commercial  Source"  column  In  the  table  and  then  to  the  company 
number  listed  in  Table  XIX. 


Table  XIX 

Equipment  Manufacturers 


1.  Atomic  Accessories,  Inc. 

611  V.  Herrick  Road 
Valley  Stream,  'lew  York 

2.  Baird  Atomic,  Inc. 

33  university  Rood 
Cambridge  39,  Haas. 

3.  Eborline  Instrument  Corp. 

P.  0.  2ox  279 

Santa  Fc,  New  Mexico 

U.  Vlctorecn  Instrument  Co. 

5606  He ugh  Avenue 
Cleveland,  Ohio 

5.  Lionel/Anton  Electronic  Laboratories 
1226  Flushing  Avenue 

Brooklyn  37,  ilev  York 

6.  Nucleonic  CoiTpomtion  of  Am. erica 
19c  Dcgraw  Street 

Brooklyn  jl,  New  York 

7.  Radiation  Co>intcr  Laboratories,  Tnc. 
5121  W.  Grove  Street 

Skokie,  Illinois 


8,  Radiation  Equipment  and 
Accessories  Corp, 

£65  Merrick  Road 
Lynb  rook ,  Hew  York 

5.  Technical  Associates 

l'*0  W.  Trovldcncla  Avenue 
Burbank,  California 

10.  Trace  rlab 

lC'Ol  Trapclo  Rond 
Walthr.m  51* ,  Massachusetts 

11.  Iluclcar-Chlcugo  Corporation 
333  East  Howard  Avenue 

Deo  Plaines,  Illinois 

12.  Franklin  Systems 
273b  Hillsboro  PoaJ 
West  Pain  Death,  Florida 

13.  Nuclear  Corporation  of  America 
T.'.strvr.’  U*  and  Control  Division 

Rlchwood  Place 
Dcnv tile.  How  Jersey 
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Using  the  Information  in  Figure  If.  to  convert  the  instrument 

sensitivities  In  Table  XVIII  into  equivalent  concentration  value*,  Tatle  XV'.I 

has  been  compiled  to  chow  the  concent rations  of  typical  Isotopes  which  can 

•  f  C  (\ 

Just  be  detected  with  lnatrJM.-r.tc  of  the  sensitivities  stated'  .  It  mur-t 
be  bomt  in  mind  that  the  dote  calibration  in  mr/hr  for  mod  Instruments 
is  done  with  a  cobalt  or  radl-ja  ground  source,  ar.d  that  the  dorc  extrapoia* 
tlon  to  other  Isotopes  and  tr.crc'  rangoc  will  not  be  the  same  for  different 
detectors. 

It  is  evident  from  the  foregoing  diccur slor.  that  theix;  appear  to 
be  adequate  commercial  ln.  tru.ner. *.a  or.  the-  msrkc-t  which  can  be  adapted  to 
give  a  reliable  Indication  of  high-level  contamination  of  water  supplies. 

The  at  thy  of  icotopic  contamination  ar.J  the  accurate  detc ralnatlon  of  cafe 
concentrations  of  rad loi cotope- s  in  uater  by  field  instruments  under  emergency 
conditions,  however,  doeu  require  .-.orec  further  research  and  development 
work.  Such  equipment  vil).  have  to  be  made  available  in  sufficient  quantity 
for  extent- lee  water  unaly.-es  at.  soon  after  on  attack  as  possible-. 
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VIII.  DECONTAMINATION  OF  HATER  SUPPLIES 

A*  dlceuosed  la  Quarterly  Technical  Report  #3  (see  Appendix  A) 
rodloIogTcnUy  contaminated  waters  have  been  subjected  to  decontamination 
procedures  of  both  conventional  and  non-conventional  nature.  Of  the  six 
radioisotopes  of  biological  importance,  all  have  been  found  to  be  amenable 
to  removal  to  some  degree.  The  degree  of  treatment  required  has  been 
predicted  in  previous  reports,  predicated  by  the  Maximum  Permissible 
Concentration  for  peacetime  consumption  of  water  or  exposure  over  a  con¬ 
siderable  period  of  time,  These  requirements  are  very  unrealistic  when 
considered  In  the  light  of  the  occurence  of  a  nuclear  Incident.  The 
removal  efficiencies  for  the  most  part  have  been  derived  from  laboratory 
data  on  small  volumes  of  water  treated  under  rigidly  controlled  conditions. 
Of  practical  Importance,  only  those  data  obtained  by  actual  treatment  plant 
operation  are  to  be  considered  reliable  on  a  major  scale.  Even  in  this 
respect,  it  must  be  pointed  out  that  these  data  were  derived  from  analysis 
of  water  supplies  which  had  been  contaminated  by  long-range  weapon  debrie. 

It  has  been  discussed  elsewhere  in  this  report  that  the  long-range  fallout 
la  more  soluble  than  short-range  or  early  fallout.  It  is  for  this  reason 
that  we  may  expect  full-rcaie  operations  to  be  more  efficient  in  radio¬ 
isotope  removal  of  early  fallout,  than  the  presently  reported  efficiencies. 

It  appears  that  water  decontamination  in  the  early  phases  of  a 
nuclear  incident  aftermath  are  not  beyond  .the  realm  of  practical  application 
of  the  basic  sanitary  engineering  principles  of  water  treatment.  The  major 
problem  involved  in  voter  decontamination  will  be  the  long-term  exposure  of 
the  surviving  population  to  concentrations  of  the  long-half lived  isotopes 
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of  biological  Importance.  The  reduction  of  these  contaminants  to  (l>-f'treablc 
levels  vlll  constitute  the  primary  long-term  problem  to  the  water  tiVntia'.nt 
field,  but  Is  by  no  means  an  insurmountable  task* 
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IX.  BCC?  3UKDHI  FROM  C0KTAMINA7ED  WATER 

In  the  analysis  of  tody  burden  from  continued  ingestion  of  water 
contaminated  with  radioactivity  resulting  from  nuclear  attach.  It  has  been 
found  that  the  developr.er.t  of  a  rational  formula  is  desirable.  A  formula 
has  been  developed  in  this  report  which  gives  values  of  body  burden  com¬ 
paring  closely  with  those  obtained  by  empirical  means. 

Many  of  the  empirical  values  used  in  this  analysis  were  derived 
from  animal  studies,  and  cay  or  may  not  be  applicable  to  man  per  se; 
however,  the  values  obtained  by  formulas  developed  in  this  report  are  in 
good  general  agreement  with  data  collected  from  experiments  with  man. 

Data  from  man  has  been  used  wherever  possible;  unfortunately  experiments 
with  man  are  not  sufficiently  comprehensive  to  servo  as  the  entire  basis 
of  on  analysis. 

The  relation  of  body  burden  to  many  of  the  radioisotopes  found 
in  fallout  is  not  easily  understood  for  short  periods  of  time,  that  is, 
when  the  study  period  is  a  very  snail  fraction  of  the  physical  half-life 
of  the  radioisotope.  In  this  report,  in  addition  to  study  of  the  immediate 
effects,  several  cases  are  analyzed  for  periods  of  time  In  excess  of  fifty 
years.  (Figures  18  and  20) 

Hunwrous  studies  liave  been  made  on  the  affinity  of  certain 
types  of  tissue  for  particular  radioisotopes.  These  selectivity  studies 
ore  of  special  Interest  to  anyone  who  wishes  to  set  up  theoretical  values 
which  may  not  bo  exceeded  without  permanent  damage,  that  is,  such  values 
as  "l-taxlnuB  Permissible  Concentration,"  "emergency  limits,"  etc.  Although 
it  is  not  within  the  scope  of  this  report  to  determine  such  values,  it  is 
of  Interest  to  note  that  allowable  body  burden  ia  the  total  concentration 
of  a  6lveu  isotope  allowed  in  the  most  critical  organ.  The  moot  critical 
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organ  is  that  organ  which  receives  its  damaging  dose  first  during  the  chain 
of  biological  processes. 

The  phenomenon  involved  in  the  rate  of  turnover  of  various  isotopes 

by  a  particular  organ  Is  of  Interest  since  this  rate  nay  be  equal  to, 

( 7 

greater  than,  or  less  that  the  rate  for  the  entire  body.1''  a  sample  plot 
of  the  critical  organ  versus  physical  half-life  is  shown  in  Figuro  16. 

Both  biological  and  physical  decay  rates  are  considered  in  this 
report.  In  order  to  demonstrate  the  difference  in  body  burden  levels  ob¬ 
tained  by  (l)  considering  only  physical  decay  and  (2)  considering  both 
physical  and  biological  decay,  l.e.  effective  decay,  values  for  both 
cases  are  plotted  in  Figures  lo  through  21  as  curves  B  and  A  respectively. 

values  plotted  ao  curve  A  are  parameters  of  the  rate  at  which 
the  radioactivity  level  In  the  body  la  reduced  by  a  combination  of  bodily 
processes  and  physical  decay.  These  values  are  therefore  parameters  of 
the  effective  decay  cates  of  the  various  radioisotopes.  The  amount  of  any 
radioisotope  which  passes  through  the  body  must  further  be  related  to  the 
amount  retained  by  the  critical  organ.  This  relation  is  represented  by  an 
f  value  which  gives  the  amount  of  a  radioisotope  taken  up  by  the  organ 
critical  for  that  Isotope. 

In  the  lntcreot  of  generality  ve  have  developed  a  method  of  compu* 
tntion  (using  AR  values)  rather  than  a  Hot  or  specific  figures  for  a  limited 
number  of  cases. 

ilo  attempt  has  been  made  to  evaluate  cell  damage  resulting  from 
various  levels  of  body  burdens.  It  is  obvious  that  the  greater  the  body 
burden  the  more  extensive  is  the  radiation  damage.  Examples  are  cited  in 
this  report  of  sooe  of  the  biological  blocking  processes  which  may  bo 
applied  to  effectively  reduce  the  body  burden.  These  An  vbIucs  will  give 
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the  highest  burden  revolting  aiilJ  fro®  drinking  contaminated  vater  while 
eating  a  normal  uncontcminated  d*et  and  breathing  uncontaninated  air. 

With  an.  easily  derived  scaling  factor,  however,  these  same  graphs  may  be 
applied  to  ingested  solids.  A/(  values  nay  be  obtained  either  from  formulas 
to  be  derived  later  In  this  see'.  Ion  or  from  the  graphs  in  Figures  16 
through  21.  The  graphs  are  easier  to  use  and  are  included  as  an  aid  to 
computation. 

No  limitations  of  physical  death  or  early  genetic  death  of  the 
individual  under  study  have  be?fi  placed  on  the  data  shown  here,  since 
both  these  Limitations  require  Sft  analysis  of  damage.  We  cannot  overlook 
the  possibility  that  At  eertaih  levels  of  contamination,  or  consumption 
physical  or  genetic  death  will  l*mit  the  continuity  of  our  data,  as  will 
any  departure  from  the  size  and  nKa  of  a  standard  man. 

Establishment  of  a  Criteria 

(a)  The  data  to  be  presented  should  be  in  the  most  useful  and 
general  form  applicable  to  the  g'lbjcct  matter. 

(b)  Analysis  of  damage  that  rctijr  occur  at  the  various  levels  of 

body  burden  is  not  to  be  attempted;  however,  comment  regarding  damage  in 
general  will  be  incorporated  to  dcPtb  to  the  conclusions. 

(c)  The  development  Qr  a  rational  method  for  evaluating  the  body 

burden  is  to  be  outlined  and  e/plalr*cd'  correlation  of  the  values 

Obtained  by  the  use  of  this  im-tifOd  and  those  obtained  by  empirical  means 
will  be  discussed. 

(d)  Assumptions  oho/M  fce  QS  f<fv  end  as  reasonable  as  possible 
in  order  to  maintain  the  lntegfMV  of  the  data  to  be  presented. 

(e)  All  values  used  !'or  the  various  calculations  shall  be 
those  derived  from  the  most  o-j*. norltatlve  information  available. 


Assumptions 

(a)  The  biological  decay  rate  used  Is  that  of  a  standard  man. 

31r.ce  the  age  arid  overall  metabolism  of  the  Individual  will  affect  the 
turnover  and  replacement  rate  of  radioactive  materials  in  the  body,  and 
therefore  the  biological  decay  rate,  it  is  necessary  to  standardize  these 
parameters  for  the  present  study. 

(b)  The  contamination  of  an  individual’s  water  supply  takes 
place  all  at  one  time,  and  no  further  contamination  occurs.  The  analysis 
of  a  fluctuating  contamination  level  would  require  that  we  assume  ratea 
of  fallout,  wind  direction  and  numerous  other  meteorological  conditions. 

In  order  to  retain  the  generality  of  our  Investigation  we  have  taken 

a  case  in  which  weather  and  attack  data  were  not  opplicatle.  cur/es 
may  be  added  to  each  other  if  new  contamination  occurs  (possibly  from 
fallout  from  e.  different  burst  or  fallout  return  from  wind  reversal), 

(c)  The  radioactive  contamination  of  the  water  supply  is  reduced 
only  by  physical  decay.  This  assumption  is  conservative  in  nature,  since 
it  disregards  entirely  the  amount  of  activity  which  might  settle  out  or 

be  consumed  by  life  forms  found  in  the  water  supply,  e.g.  plankton,  algae, 
fish. 

(d)  Effective  decay  rate6,  even  if  extrapolated  from  animal 
data,  are  applicable  to  man.  In  the  absence  of  experimental  data  on  man, 
the  animal  data  'used  was  the  best  available.  It  nwy  be  of  inferest  to  note 
that  many  contemporary  investigators  are  not  satisfied  with  many  of  the 
extrapolations  of  animal  body  burden  data  that  have  heretofore  been  accepted 
as  correct,  or  at  least  useable,  for  the  case  of  humans. 

(e)  Tne  daily  intake  of  water  is  constant  and  takes  place 
instantaneously.  Virtually  all  of  an  Individual's  water  intake  occurs 


in  sixteen  hours,  and  over  half  occurs  in  eight  houro;  however,  the  error* 
associated  with  this  amplification  tend  to  he  compensating  if  the  individual* 
drinking  hah  its  are  consistent*  For  reasons  Of  generality  vc  axe  forced  to 
consider  that  the  man  under  study  is  a  creatul0  of  habit.  By  drawing  a 
smooth  curve  through  the  dolly  points  we  have  essentially  integrated 
over  the  time  periods  between  the  daily  values  and  have  described  a  contln- 
uous  function  as  if  the  individual  had  spread  hi®  intake  over  the  eutire 
day.  Tho  combination  of  the  two  effects  result0  in  a  close  approximation 
of  the  actual  build-up  of  body  burden. 

Development  of  a  Rational  Method  of  Analysis  i1*'  body  Burden  for  Continuous 
Consumption  or  Water  Contaminated  with  RadioRi'Uve  Element* 

Based  on  the  Gelger-Nutall  Law  of  radioactive  tlecay  we  may  say  that  the 
amount  of  activity  present  at  any  time  1*  equal  to 

At 

A  •  A  e  v 
o 

when  only  physical  decay  is  considered.  iS  physical  decoy  constant 
for  a  particular  radioactive  element,  and  Aq  la  amount  of  activity 
initially  present. 

Furthermore,  if  only  biological  decay  1*  considered  we  may  say 

that 


where  is  the  biological  rate  of  elimination  constant. 

If  ve  define  "Effective  Half  Life”  p«  th<?  time  required  for  a 
radioactive  element  fixed  in  the  tissues  of  nil  animal  body  to  be  diminished 
by  fifty  per  cent  as  a  result  of  the  combined  action  of  radioactive  decay 
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and  biological  elimination,  then  we  nay  say  that  the  effective  half-life 

»<7« 


Biological  half-life  x  Radioactive  half-life 
Biological  haU-life  +  Radioactive  half-life  *  eff 

An  effective  decay  constant  (keff)  may  be  defined-/^ 


Let  us  assume  initially  that  for  any  activity  taken  into  the 
body  the  entire  amount  will  be  fixed  in  tissue.  It  is  then  apparent  that 
immediately  before  initiation  of  consumption  of  contaminated  water,  th-'  body 
burden  iu  t-quai  to  whatever  radioactivity  is  present  a3  a  result  of  pre-  . 
vlous  activity.  For  the  development  of  our  argument  we  ahall  take  thi3 
Initial  body  burden  to  be  zero  and  correct  it  in  later  developments  for 
the  caoe  of  Initial  body  burden  not  equal  zero. 

The  following  reasoning  may  then  by  applied: 

(1)  Just  before  initiation  of  consumption  activity  in  the  body 
equals  zero. 

(2)  Immediately  following  the  3tart  of  consumption  the  activity 
in  the  body  equals  A^,  that  is,  the  activity  per  unit  volume  of  the  water 
(nc/ml,  ntorou/gal,  etc.)  times  the  volume  of  water  consumed.  Since  we 
have  assumed  that  the  amount  of  liquid  is  constant  we  moy  drop  the  volume 
term  and  ^incorporate  it  later  in  the  development.  Furthermore,  by 
considering  Aq  equal  to  unity  we  will  develop  an  equation  which  when 
multiplied  by  the  actual  Initial  activity  will  give  activity  present  at 
any  given  day  n.  Wn  wish  to  develop  a  plot  which  is  independent  of  both 
amount  consumed  and  initial  activity,  and  is  dependent  only  on  the  radio- 
lrotopc  ill  question. 
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(3)  Immediately  prior  to  consur.ptloa  on  the  ee»,c,n<l  3ay  the  body 
burden  would  be 

Ai  “  V 


and  Just  after  consumption  the  body  burden  is 

Aeff(l)  -*Jl> 

AJ  .  A*  eff  ♦  Ae  p 

L  O  Q 


(U)  Just  before  consumption  or,  the  third  day  \Ha  body  burden 


would  be 


-  AJ.  'k'frl  -  (Ao,  A«"‘  « 


and  immediately  following  ingestion,  the  expression  for*  t',vV  burden  may 
be  written 

•1  1  -A  l  A  l  A  1  A 

Ag  »  (Aq6  P  )e  P  +  (AQe  eff  ♦  Aoe  P  )« 


(5)  Continuing  this  reasoning,  a  recursion  fO)W>la  1°  developed 
of  the  fora 

A  .  A  I."*'"1  ♦  «  *  V‘ 

n  o 


♦  e 


-ueff  +  (n*l5\.lj  -n\_l 


•nil 

P"  ♦  e  ]  *  V  P 


However,  rather  than  deal  with  the  (n  +  l)  terms  nccebsHf/  1®  evaluate 
this  equation  at  n  days  ue  simplify  the  e-^.tion  os  folld'"'0 


hn  ~a 

A  .  A  b  (H-  '  .r.) 
n  o  b  -  p  ' 


where  h  •  A  t 
0 


and 


p  -  Aoe 


*yi> 
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(6)  This  general  expression  for  Aq  is  the  value  plotted  as 
curve  A  in  Figures  16  through  21. 

(?)  We  may  now  incorporate  into  our  formula  the  portions  of  total 
tody  burden  which  have  thus  far  been  neglected. 

(a)  Since  we  have  considered  Aq  to  be  unity,  we  may  find 
the  total  body  burden  by  multiplying  Aq  by  the  true  original  activity, 

Aq  tnte»  the  volume  consumed  daily,  D. 

Therefore  Body  Burden,  (ggS)  -  Aq  ^  (ggl)  x  D(M^)  x 


A  Of  consumption* 

n  body  ' 


(b)  The  assumption  that  all  the  activity  taken  into  the 
body  is  immediately  incorporated  into  the  intestine  is  highly  improbable. 
We  may  however,  consider  the  amount  Incorporated  into  tissue  ns  n  fraction 
of  the  ingested  activity  reaching  an  organ  (fy)  and  multiply  the  body 
burden  of  (a)  above  by  f^  to  obtain  a  more  reasonable  value. 

Thus  Body  Burden..  *  f  x  A  xDxA 

*1  w  o  true  n 

Some  values  of  t  for  Isotopes  under  consideration  are  given  in  Table 

(c)  While  reducing  the  bedy  burden  by  the  factor  f  we  have 
neglected  the  activity,  A^,,  which  is  present  in  the  intestinal  tract 
on  day  n.  Thio  activity  is  an  addition  to  the  body  burden,  therefore: 

Body  Burden, XI  -  x  *0  true  x  D  x  Ar  ♦  AQ, 


where  by  the  reasoning  and  definitions  previously  expressed 


o  true 


A  (n-0.25)  -1.25V 

P  )(e  ^ 


Since  the  tine  required  for  ingested  material  to  pass  through  the  body 
is  31  hours  or  about  1.25  days,  the  A^j  tem  is  strictly  transitory. 

Tt  nay  or  may  not  be  negligible  depending  on  the  period  in  which  it  takes 
place.  On  the  first  day,  for  example,  it  would  be  the  largest  portion  of 
the  total,  but  by  the  time  the  fractional  half  life  (n/T^)  is  equal  to 
orie,  the  term  may  in  some  cases  be  considered  negligible,  (n  equals 
time  in  days  since  beginning  of  consumption  and  equals  radioactive 
physical  half-life  in  days.) 

(d)  The  only  tern  still  missing  from  the  mathematical  expression 
Is  that  for  the  activity  present  in  the  body  at  the  lime  drinking  of 
contaminated  water  started  (Aj).  That  which  was  already  present  in  the 
tissue  at  time  of  initial  ingestion  will,  after  n  days  of  effective  decay, 
equal 


V 


A 

e 


rf° 


That  present  in  the  gaatro-inteetinal  tract  at  the  time  of  initial  ingestion 
will  give  an  additional  activity  equal  to 

AJ(  •  fy  x  Age  offU 

where  A^  is  the  pre-ingestion  activity  in  the  intestinal  tract. 

(8)  After  combining  the  various  components,  the  expression 

for  body  burden  from  any  particular  radioisotope  at  day  n  may  be  written: 

-\  (n-1.25)  -1.25V 

Body  Burden  ■  (f^)(Ao  tn true®  ^ 

(*!  •  Vo>*  •  *  Aot  *  \  *  \ 
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The  full  explanation  of  this  development  has  been  presented  in 
order  to  allow  the  user  of  the  quatlons  to  decide  which,  if  any,  terms 
may  be  neglected  for  any  specific  analysis.  Should  the  individual  In 
question  be  an  average  IT.  S.  citizen  who  has  sought  shelter  during  an  attach, 
the  term 

*i  *  fA  '  "  *T  *  *x 


will  for  practical  purposes  equal  zero.  As  previously  noted,  the  value 
of  AQI  will  generally  be  negligible  after  a  relatively  short  period  of 
consumption.  If  the  (A^  ♦  A^)  and  terms  are  dropped,  the  practical 
expression  for  body  burden  is  (fyH*0  true^D^An^*  117,6  ACI  tcrn»  h<*reve*“> 
may  be  of  major  ir.jvrtance  when  the  physical  half-life  io  much  greater 
than  the  effective  half-life. 

It  muot  nlco  be  streoced  that  to  evaluate  the  effect  of  Aqj, 

At  and  on  any  individual  that  individual’s  hlBtory  of  exposure  to 
radiation  must  bo  known.  It  is  only  by  thorough  Knowledge  of  the  develop¬ 
ment  of  this  rational  method  that  one  may  make  a  reasonable  calculation  of 
body  burden . 

Reduction  of  Body  harden  by  Blocking 

It  has  Iron  found  that  in  certain  coses  uptake  of  radlolcotopos 
may  be  reduced  by  biological  blocking.  For  example,  by  administering  a 
soluble,  non-toxic,  inert  compound  of  iodine  one  can  prevent  uptake  of 
radioactive  iodine  Vy  the  tlyroll.  Ten  to  one  hundred  milligrams  of 
potassium  iodide  token  at  one  time  will  prevent  the  uptako  of  a ay  more 
iodine,  radioactive  or  otherwise,  by  the  thyroid  for  about  one  week.  To 
continue  this  effect,  it  is  believed  that  an  additional  '.00  milligrams  evory 


two  days  would  be  required 


(77) 


Blocking  lovers  the  voluo  of  f  in  body 
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burden  calculations. 

Id  the  case  of  Strontium  89  and  Strontium  3C,  which  behave  like 
calcium,  it  1*  believed  that  blocking  nay  be  completely  or  partially 
accomplished  by  consuming  a  calcium  compound  far  in  excess  of  normal 
requirements.  Ibis  in,  however,  a  conclusion  made  on  the  baaia  of  extrapo¬ 
lated  data;  the  only  case  in  which  it  has  definitely  ceen  proved  that  blocking 
Is  effective  is  that  of  iodine. 

Should  the  body  be  Initially  deficient  in  a  given  stable  isotope 
(possibly  from  deficiency  of  a  common  mineral  in  the  pre-ingestion  diet) 
the  value  of  f  for  the  radioactive  isotope  of  that  mineral  or  a  chemically 
.timllar  mineral  may  be  lncreaoed  because  of  the  increased  affinity  of  the 
body  for  that  mineral. 

Blocking  should  not  be  considered  a  panacea  to  the  problem  of 
body  burden,  since  it  may  result  in  the  cure  being  worse  than  the  disease. 

In  the  iodine  example,  for  instance,  the  prescribed  mounts  of  potassium 
iodide  constitute  a  daily  Input  of  iodine  into  the  blood  one  thousand 
times  normal.  Although  it  is  believed  that  for  the  majority  of  the  U.  8. 
population  this  treatment  would  not  be  harmful  for  to  to  three  weeks, 
the  secretion  rate  of  the  thyroid  is  decreased.  Extended  treatment  will 
produce  other  still  more  harmful  side  effects  vhich  bare  not  yet  been 
fully  evaluated.  Blocking,  however,  may  play  an  important  role  in  com¬ 
puting  body  burden  valuco  and  should  not  be  overlooked  when  selecting 
f  values. 

Discussion 

(l)  In  figures  16  through  21,  curve  A  gives  the  value  resulting 
fron  the  effective  decay  rate.  Curve  B,  included  for  purposes  of  comparison, 
represents  the  rate  at  which  physical  decay  alone  take*  place.  It  accns  to 


point  out  that  because  of  the  body's  tendency  to  eliminate  certain  radio¬ 
isotopes,  or  not  to  incorporate  then  into  tissue  at  all,  body  burdens  are 
significantly  lowered. 

(2)  T^ble  XXTY  gives  a  fev  selected  values  of  f  for  six 
radioisotopes  and  various  organs.  These  are  values  which  result  from  a 
normal  uncontaminated  diet.  Values  of  f  nay  vary  from  this  table  according 
to  the  excess  or  deficiency  of  certain  minerals  in  the  body.  Note  that  the 
concentration  of  specific  Isotopes  must  be  known  before  body  burdens  can 
be  computed. 


Table  XX 

Selected  Values  of  t  for  Various  Organs 


Isotope 

Organ 

fv 

Sr09 

Bone 

2.5  X  10* : 

Sr90 

Bone 

2.5  x  10*: 

Ru106 

Kidney 

2.0  x  10* ! 

x13l 

Thyroid 

2.0  x  10* 

Whole  Body 

1.0  x  10° 

C.137 

Liver 

U.8  x  10* 

B.1U0 

Bone 

7-0  x  10*: 

Motes: 

(1)  The  f  values  shown  are  for  critical  organs.  The  critical 

organ  will  determine  the  A  since  It  Is  the  one  used  in  computing  X  ... 

n  cx  1 

(2)  The  f  value  for  the  entire  tody  is  the  ruunnation  of  f 

v  w 

values  for  all  the  Individual  organs  concerned. 


I  J 

FRACTIONAL  HALF-LIFC 


Figure  l6.  vereuo  Tlisc  for  Iodine  1J1 
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1  2 

FRACTIONAL  HALF-LIFE 


Figure  18.  versus  Tine  for  Strontium  90 


Figure  19.  Aq  versus  Time  for  Ruthenlua  106 


IV  kvr  »>  VMUJir,* 
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FRACTIONAL  HALF-LIFE 


Figure  21.  An  versus  Time  for  Barium  l40 
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A 

\j 

H 

Ao 

a; 

\ 

\ 


a 


enpirical  constant  for  sublimation  or  vaporization  reaction 


radioactivity  after  come  period,  of  decay 

radioactivity  already  present  in  the  Intestinal  tract  at  time 
of  ingestion  of  contaminated  water 

radioactivity  present  In  the  intestinal  tract  at  some  time  after 
Ingestion  of  contaminated  water 

radioactivity  already  incorporated  Into  body  tissue  at  time  of 
Ingestion  of  contaminated  water 


a  parameter  of  the  effective  decay  rate  in  the  body  of  a  given 
isotope 


initial  radioactivity  per  unit  volume  of  contaminated  water 
taken  Into  the  body 


(o  f  Nk) 

t-%1  A 

MT* 


■ 


radioactivity  present  In  body  tissue  resulting  from  decay  of 
that  activity  initially  present  in  tissue  (A^;  after  some  time 
of  effective  decay 


activity  at  time  t 

activity,  In  disintegrations  per  second,  of  the  radionuclide 
In  a  target,  after  the  nuclide  has  been  removed  from  the  flux 
for  a  period  9j 

radioactivity  incorporated  into  tissue  from  that  activity  originally 
present  In  the  Intestinal  tract  <V 

cloud  horizontal  semi-axis 


a  fireball  horizontal  eemt-axis  at  ground  zero 

o 

*2  fireball  horizontal  semi-axis  at  2 

B  empirical  constant  for  sublimation  or  vaporization  reaction 

B  ratio  of  fission  to  total  yield 

Body  Burden  radioactivity  incorporated  into  body  tissues  as  a  result  of 
Ingestion  of  radioactivity 


b 


cloud  vertical  semi-axis 


c 

D 

D 

DX 

d(t,  1) 
FDr(t) 
FPr(t) 
f 

fw 

e 

b 

yt) 

ljft) 


k 


k 

k 

a 

\ 


G10GSAXC 

(Coatlaitd) 

fireball  vertical  send-*/!*  at  ground  zero 
fireball  vertical  semi~tala  at  2 

empirical  constant  for  zublination  or  vaporization  reaction 
instrument  response  faetor 
daily  water  intake,  volume 
exposure  dose 

decay  correct  ion  factor  for  H  ♦  1  hour 
fraction  of  device  contour  ratio 
f ia6ion-produet  contour  ratio 

thermal  flux  in  neutronr  per  square  centimeter  per  second; 

fUgacity  of  the  element  in  the  liquid  phase 

a  coefficient  relating  invested  radioactivity  to  Body  Burden 

acceleration  due  to  grertty 

cloud  center  height 

fallout  intensity  at  time  t 

air  ionization  rate  per  fission  at  3  ft  above  an.  infinite, 
ideal  plane  for  a  uniform  distribution  of  the  normal  fission 
product  mixture 

air  ionization  rate  per  fission  at  3  ft  above  an  infinite, 
ideal  plane  for  a  uniform  distribution  of  the  condensed 
fission  produdt  mixture 

air  ionization  r*ta  per  fission  at  3  ft  above  an  infinite, 
ideal  plane  for  a  uniform  dlstrioution  of  the  neutron  Induced 
activities 

relative  abundance  of  tie  isotope  from  which  the  radionuclide 
is  formed 

arbitrary  constant 

fireball  horizontal  set! -axis  expansion  constant 
fireball  vertical  semi -axis  expansion  constant 
Henry's  Law  constant 


maaatam 
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Mr(t) 

a 

m 

N 

V 


»J(A) 


n 

n(i) 

n. 


n 


i 


n,5A(t) 


n 


HT 

inverse  time  constant 
mass  contour  ratio 

mass  of  fallout  per  unit  area  at  any  downwind  distance,  X 
mass 

disintegration  multiplier 

the  total  number  of  atoms  of  the  element  In  the  target; 
number  of  atoms,  or  moles,  of  f lesion  products  per  unit  area 
mole  fraction,  of  element  J  la  the  liquid  phase 
mole  fraction  of  element  J  in  the  vapor  phase 

number  of  a  tome  (of  the  end  member  of  cass  chain  and  condenccd  era 
the  outside  Of  tho  particle)  which  land  per  square  foot  of  ground 

number  of  moleo  of  vapor 

total  moles  of  liquid  carrier 

abbreviation  of  nj ^(t) 

amount  of  element  J  condensed  on  the  surface  of  the  solid 
particles 

amount  of  element  J  In  the  vapor  phase 
abbreviation  of  "Ja«> 

number  of  moleo  of  element  J  with  mass  number  A  dissolved 
In  the  n(i)  moleo  of  liquid  carrier,  which  is  the  particle 
In  the  liquid  phase,  prior  to  solidification 

number  of  moleo  of  element  J  with  mass  number  A,  which  Is 
mixed  with  nolcn  from  other  mass  chains  to  form  n  moles  of 
vapor  and  not  condensed  in  the  liquid  carrier 

empirical  constant  In  particle  falling  velocities  relation 

sublimation  or  vaporization  pressure 

partial  pressure  of  the  liquid 


<1 

<*X 

B 


ro(A,  t) 
r;(A,  t) 
rx(0 
T 
T 


t 

t 


V 

V. 


w 

x 


sublimation  pressure 

<V9 »; 


empirical  constant  in  particle  falling  velocities  relatioa 

terrain  shielding  factor 

molar  gas  constant 

initial  fireball  spherical  radius 

fractionation  number  of  the  first  period  of  condensation 

fractionation  number  of  the  second  period  of  condensation 

gross  fission  product  fractionation  number 

absolute  temperature 

half-life  of  the  radionuclide  formed 

time  required  for  a  radioactive  element  fixed  in  the  ti66ue 
of  an  animal  body  to  be  diminished  50  per  cent  as  a  result 
of  combined  action  of  physical  decay  and  biological  elimination 

radioactive  (physical)  half-life 

time  of  irradiation 

time 


time  of  arrival  for  particle  falling  out  of  the  cloud 
time  of  cessation  for  particle  falling  out  of  the  cloud 
particle  falling  time  in  the  fireball 
particle  rising  time  in  the  fireball 
molar  volume 

average  particle  falling  velocity 
vind  velocity 

instantaneous  particle  falling  velocity 
weapon  yield  in  Xlloton  TNT  units 
longitude  axis  in  ground  coordinates  system 
longitude  axis  in  cloud  coordinates  system 


x 
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JrM(t) 


latitude  axla  In  ground  coordinate  system 

total  amount  of  radionuclides  for  the  entire  chain  yield  of 
=«33  meter  A 

total  tis,y. int  of  radionuclides  of  element  J  present  at  time  t 

total  number  of  atoms  of  element  J  of  all  mass  numbers 

maximum  half -width  of  the  crossvind  distance  on  a  1  r/hr 
contour 

latitude  axis  in  cloud  coordinate  system 

amount  of  radionuclides  of  element  J  and  mass  number  A 
present  at  time  t  after  fission 

number  of  atoms  of  element  J  of  mass  number  A  not  condensed 
in  the  liquid  carrier 

altitude  axis  of  ground  coordinate  system 

yield  dependent  multiplier  in  fireball  altitude  equation 

altitude  axis  in  cloud  coordinate  systea 

scaliest  Intercept  of  particle  falling  slope  with  cloud 

largest  Intercept  of  particle  falling  slope  with  cloud 

particle  size  parameter 

maximum  or  minimus 

X/h 

a  constant  expressing  the  biological  rate  of  elimination  of 
radioactivity 

a  constant  expressing  the  rate  of  attenuation  of  radioactivity 
in  the  tody  as  a  recult  of  both  physical  decoy  and  Moiogical 
elimination 

disintegration  constant}  the  fraction  of  the  number  of  atoms 
of  a  radioactive  nuclide  which  decay  in  unit  time 

activation  cross  section  in  square  centimeters  for  00  m/sec 
neutrons 


0 


tine  of  decay 
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APPan>.TX  a 

Ecrr-pts  f’-o-a  Qu-.rt  riy  Te.obn-.caI  Prpcrt  tl 

F •  rroVi.nltJ~.ry  Evaluation  of  Wotcr  Supply  Ccyitamlnatlrr.  for  Selected  Clttea 
1.  Ir.t  redder  lor. 

Tr.  order  to  have  some  taels  fc-r  evaluating  tho  relative  hazard  of 
nuclear  fallout  to  the  water  supplies  of  different  cities,  attacks  have  been 
assumed  on  c-ach  city  thv-t  would  result  In  the  ssrr.e  intensity  over  each 
watershed  invol/ed. 

From  this,  the  centred  ra  tion  of  'he  voter  due  to  the  suluble  portion 
oi  the  fallout  has  beer,  calculated  for  K  1  hour.  Simplifying  assumptions 
were  made  in  regard  to  an  average  intensity  over  the  specific  watersheds  and 
c.j  to  the  time  factors  involved.  There  ssctjnptior.;  were  made  to  provide  a 
hasie  fc-r  the  calculations  of  the  activity  concentrations  shown  in  Table  IV. 
There  values  are  shown  principally  to  d-.  mono  t  rote  tV.c  approach  being  taken 
5nd  to  give  s.-a  iiua  of  the  maximum  contamination  that  could  be  expected 
immediately  following  a  nuclear  emergency. 

The  approach  to  be  taken  in  future  study  will  be  more  specific. 
Factors  that  will  be  taken  into  account  to  give  comeuhat  more  reasonable 
values,  will  be  .  .'l)  a  realistic  sup.-rimposlt  ion  of  intensity  contours  upon 

tho  specific  vote rshe Jt  with  regard  to  eu.sor.ai  wind  directions  and  target 
sites,  (i)  consideration  of  decay  factors  to  gi.c  c.vprcted  activity  at  times 
hater  than  ’1  +  l  hour,  ;  3)  consider*  f  lor.  of  vat-,  r-hed  characteristics  such 
as  runoff  co^t  l  lc lent,  feed-stream  •ciocltiee,  ivser/oir  draft,  transport 
pV  :.c-.w-n» ,  nnd  reservoir  mixing  propvities,  'L)  integration  of  intensity 
contours  to  give  h  no-~e  precise  value  for  up.  clfic  isotopic  activities  or. 
the  watermen  e.r  d  (5)  a  conversion  of  ictetialty  to  activity  for  sites  off  the 
dc.ro  1  r. t  ar, lu  er.p toying  a  known  rcia*. ;or.uhlp  to  tho  predetermined  values 
directly  do-no ir.d . 
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T&fcle  IV 

Activity  C.ar.ccnt  rations  lr.  V»'j t<*r  Selected  Isotope? 
Folio'- .;i»t  3  MI  s.r  attacks 

Hailonnclide  lorctr.trat Ion  in  Vinter  (uc/ml) 


**-  i  v.’ 

Sr-09  I:-90  R-i-15c  I-l 

31 

0.-13-  *a-l:-0 

T  ir.vcr 

i 

2.7  x  10*“'  3.-  X  10*-’  3.0  x  10*-  3.0  X 

10*1 

1.9  >-  n**  3.2  x  10*1 

1’  ->lo . 

ii 

f.l  x  10* 1  t.5  >•  10* 5  5.6  v  10*'  5.7 

3.5  X  11*-  6.1 

" *  *  1 ii h 

•» 

7.3  x  10'h  7.3  x  10*6  e.2  x  10* 7  8.0  x 

X0*J 

4.9  >.  1C*-  5.7  X  10* 3 

Mi  t.n . 

lot  Applicable 

to:: 

- 

6.5  x  10*J  6.1  x  10*5  5-0  x  10’*  5.0  x 

io*2 

3.1  x  10*"  *.6  x  10*? 

le/s.- 

ii 

6 . 5  x  10  1  1 . 1  x  10  ■  7.0  x  10*  '  3-0 

3.1  x  10*'  3-6 

'  <?w  ':ork 

i 

10.1  x  10*’’  9-  3  x  I0’6  i.l  x  10***  1.2  x 

10*£ 

7.2  x  10**  8.3  x  10* 3 

:j.  y. 

ii 

3.2  x  iO*2  2-9  x  10*U  3.5  x  10* :  3.5  x 

10* 1 

2.e  x  i:*~  2.6  x  10*-' 

^  C  |i  ^  ^  # 

.  prir.gf  ield  I  1.2  x  .10* 3  l.j  x  10*  '  1.6  x  10*  1.6 -x  10*  8.5  >.  10*''  9 -9  x  10  “ 

X  «  p  i  -  * .  m  % 

vi  9.3  x  i.0  8.5  x  .0*  i.l  v  10*-  1.1  6.5  x  10**  7-6  x  10  * 
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By  employing  the  aforementioned  characteristics  and  variables.  It  la 
believed  that  activity  concentrations  will  be  reached  that  vould  approach  a 
realistic  case. 

2 .  Calculation  of  Ac*  Wtty  Concentration 

The  calculation  of  concentration  of  activity  used  In  the  preliminary 
study  Is  based  upon  n  cons ^ration  of  the  area  of  the  water  surface,  the  volume 
of  vater,  and  the  total  watershed  area.  In  the  case  of  Duluth,  Minnesota,  where 
the  entire  Municipal  water  e'ipply  is  obtained  from  Lake  Superior,  only  the  depth 
of  the  lake  a.  the  Intake  point  was  considered. 

Current  Information  on  the  characteristics  of  the  watersheds  was  obtained 
in  each  case  from  officials  of  the  municipal  vater  supply  system  Involved  or 
their  engineering  consultant .  Data  used  in  this  report  lo  tabulated  In  Table  VI. 

TABLK  VI 

Characteristics  >»f  Watersheds  Serving  Selected  Cities 

Total  Mu torched  Total  Water  Total  Water 


Denver,  Colo. 

Area 

jf 

hi 1 

Surface  ( 

ft.2) 

Volume 

(liters 

8.29 

X 

1010 

U.36  r 

V? 

2.34 

X 

10U 

Houston,  Tex. 

7.8 4 

X 

10l° 

7.85  x 

108 

2.50 

X 

101L 

New  York,  H.  Y. 

4.24 

X 

1010 

1.37  x 

109 

1.90 

X 

1012 

Springfield,  Mass. 

1-30 

X 

109 

1.68  x 

107 

1.97 

X 

1010 

Ac  an  example  of  »•!•«  procoduro  used,  the  derivations  of  the  concentra¬ 
tions  in  the  Denver,  Coloi**'  water  supply  prior  to  any  municipal  treatment  are 
presented  here  in  detail.  'iV  other  cities  listed  in  Table  TV  vere  analyzed  In 
the  same  manner  us  Denver,  »-n:ipt  for  Duluth  which  is  presented  separately. 
Danver  Case  I 

In  case  I,  only  radioactive  fallout  deposited  directly  on  water 
surfaces  was  considered  ai  ",u  «’urco  *-he  water  contamination.  Tho  following 
assumptions  were  made  prior  10  *•*««  calculation. 
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1.  5  HT  bomb  attack. 

2.  Vlnd  speed  of  15  mph. 

3.  ->n  intensity  of  5,500  r/hr  over  the  entire  watershed. 

4.  Complete  homogeneous  solution  of  the  soluble  portion  of  the  fallout  In  the 
voter  within  one  hour. 

5.  A  ratio  of  ooluble  radioactive  ato ns  per  square  foot  to  intensity  was 
chosen  as  that  value  for  each  isotope  at  51*5  miles  directly  downwind. 

C.  That  there  will  be  no  change  in  concentration  of  radioactivity  either  duo 
to  reservoir  draft  or  feed-in. 

7.  That  only  Sr-89,  Sr-90,  Ru-106,  1-131,  Cs-137,  and  Ba-llO  vould  be  conoldcred. 

ft.  That  the  volume  and  area  of  feed  st roans  will  not  affect  the  final  concentra¬ 


tion  appreciably. 

The  coluble  activity  density,  for  each  of  the  six  isotop co  being  cod- 
flldcred,  in  atoms/ft2  waa  calculated  as  follows: 

Given:  l(l)  >  5,500  r/hr, 


•  6.17  X  10  atomo/ftyr/hr. 


fy.(90)  iq  ,  2 /  a. 

a  11.1  x  10  atomc/ft  Jr/hr, 

JI.(106)  10  2 1  . 

"  ^.87  x  10  atoms/ft  Jr/hr, 

<(137)  10  .  2/  , 

-  -y^yy  -  a  10 . 7  X  10  QtOKs/ft  Jr/hr, 

<(l37)  10  2/, 

■  ■  a  9.19  X  10  atonc/ft  j r/hr,  and 

N.  (l**0)  i_o  ,  2  1  . 

—  a  12.-  x  Kj  atoos/f ty r/br, 

where  the  subscript  4  is  the  value  of  the  particle  site  parameter,  a,  correcpondirg 


HuU37)  10 

'  -y^yy  -  a  10.7  X  10  at  OK. 


to  the  downward  distance. 


By  multiplying  5,500  r/hr  by  <(A)/i(l)  for  each  isotope,  the  total 
number  of  soluble  atoms  of  each  isotope  per  square  loot,  K^(A),  will  be  obtained. 
This  value  for  each  isotope  was  found  to  be  ns  follows: 

<(A)  -  KD 


<(89) 

m 

3.1*0  X  10l4  Btoms/ft2 

<(*> 

m 

6.10  x  10l4  (4 toms/ft2 

<(106) 

• 

2.68  x  10l4  ntoms/ft2 

<(131) 

■ 

5.e9  x  10l4  ntoms/ft2 

<(137) 

• 

5.05  x  10l4  ntoms/ft2 

<(140) 

■ 

6.82  x  10l4  atoms/ft2 

Multiplying  the  <(A)  value  for  cocli  of  the  isolopt-s  by  the  total  area 
of  water  surface  will  give  the  total  number  of  atoms,  N(A),  of  each  isotope  in 
soluble  form,  deposited  in  the  water. 

<(A)  x  Area  .  HU)  atoms 


N(89) 

■ 

3.1*0  x  10l4  X  <1.36  x  108 

■ 

11*  .0  x  1022 

N(90) 

• 

6.10  x  10l4  x  (|.36  x  108 

■ 

26.6  x  1022 

K(  X06) 

■ 

2.68  x  10l4  x  »* . 36  x  10® 

« 

11.7  x  1022 

N(130 

9 

5.89  x  10l4  x  (t.36  x  10® 

m 

25.7  x  1022 

N(137) 

m 

5.05  x  10l4  x  ><.36  x  10® 

m 

22 

22.0  x  10 

N(lUO) 

m 

6.82  x  10l4  x  <*.36  x  10® 

m 

22 

29.7  x  10 

atone 
atoms 
atoms 
atoms 
atons 
a  tons 


Thus  Knowing  the  total  number  of  soluble  storm  of  each  isotope;  divlolon  by 


tba  total  vater  roliaae,  VV  In  liter*  will  giro  the  concentration  of  each  Isotope 
la  atooa  per  liter. 


C ,  m  ll (a)  atoms 
V  liter 


* 

Ca 

9 

14.8  x 

1022 

• 

6.33  X  1011 

9? 

2-3*-  x 

10U 

# 

9 

26.6  X 

1022 

m 

11.4  x  1011 

90 

2.34  X 

To53- 

* 

c  ^ 

11.7  X 

1022 

9 

4.99  x  10U 

106 

2.34  X 

To11 

* 

25.7  X 

1022 

. .  11 

C131 

■ 

2.34  X 

m 

11.0  x  10 

c* 

22.0  x 

io22 

m 

9.40  x  1011 

137 

2.34  x  10U 

C*| 

29-7  x 

,«22 

10 

m 

12.7  x  1011 

l40 

2.34  x  IQ3* 

ter 


<  J 


We  may  now  convert  the  activity  In  atoms/liter  to  mlcrocuries/mllliiiter 
(pc/nl)  tiding  an  appropriate  conversion  factor  for  each  Isotope. 

First,  to  change  from  atomb/llter  to  disintegrations  per  minute  per  liter, 
dP°?  >  Me  following  relationship  Is  used. 

°«  *  CAX  (22) 

Where:  Cn  -  Concentration  of  activity  In  d ym/t. 

CA  “  Concentration  of  activity  In  atoms/i,  and 

^  *  Character lstlc  time  constant  defined  by  the  e<juatioo 

\  m 

V2 

Where:  •  physical  balf-llfe  of  tb#  specific  Isotope  In  Mnutes. 
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1?  , 

Site*  there  ere  2.22  x  10  '  dpn/curle,  a  constant  for  all  Isotopes, 
division  of  CQ  'ey  this  valu*  vili  give  the  concentration  In  curleo/llter. 


CA* 


f  2.22  X  10U 


2.22  x  10 


where  Is  th=  concentration  of  radioactivity  In  curies  per  liter. 

Multiplication  of  hy  10^  will  convert  curlcs/liter  to  mlcrocurie.?/ 
milliliter.  To*  entire  conversion  factor  for  etch  Isotope  la  then 

X. 


K. 


2.22  x  1CT 


The  conversion  factor,  X^,  for  ench  of  the  six  isotopes  being  studied 


ate; 


Kq9  .  4.25  x  10' 15 


90 


^106 


K13l 


2.12  x  10 


-  5-95  x  10 


-17 


■16 


2.72  x  10 


-14 


KL37  •  1.93  x  10 


K140  "  1'70  x  10 


■17 


■14 


By  Eulsiplylng  each  inversion  factor  hy  the  respective  concentration 
of  each  isotope  ic  atoms  per  liter,  the  following  values  were  obtained. 

0(69)  -  2.7  x  10*3  uc/al 
C(50)  -  2.4  x  10*5  wc/ml 
C(l06)  •  j.C  x  10*1*  uc/ml 
C(l3i)  -  3-0  x  10*2  itc/tal 


Case  II 


* 

. 

C(137)  -  1.9  *  10*5  Hfi/ml 

C(lUo)  -  2.2  x  10*2  uc/nl 

In  this  case  the  soluble  fraction  of  the  fallout  particles  deposited  over 
the  entire  watershed  was  assuracd  to  have  been  transported  to  and  homogeneously 
dispersed  in  the  streams  and  reservoirs.  Other  assumptions  were  the  same  as 
for  Case  I.  To  continue  the  Denver  example.  In  Case  II,  the  watershed  area  of 
8.29  x  1010  ft2  was  multiplied  by  the  IT^  value  for  each  isotope  to  give  the 
total  number  of  atorto  of  each  Isotope.  The  remaining  calculations  were  carried 
out  in  the  sane  fashion  aa  for  Case  I. 

Duluth.  Minn. 

The  entire  water  supply  for  the  city  of  Duluth  Is  obtained  from  Lake  Superior, 
and  hence,  this  situation  la  different  from  that  at  the  other  cities  studied.  The 
Duluth  Intake  line  extends  out  llj6  feet  Into  the  lake  and  terminates  at  a  lake 
depth  of  72  feet.  The  fallout  deposited  over  this  Intake  was  considered  to  be 
evenly  dispersed  throughout  this  depth  of  water,  and  from  this  assumption  the 
number  of  atoms  per  unit  volume  and  subsequently  the  activity  In  yc/ml  was  calcu¬ 
lated. 

j.  Conclusions 

In  all  of  the  caocs  thus  far  studied,  the  time  factor  has  not  been  considered. 
Even  at  on  assumed  wind  velocity  of  15  miles  per  hour,  a  number  of  hours  would 
have  to  elapse  before  the  arrival  of  particulate  mntter  at  the  far  reaches  of  a 
diotant  watershed.  Even  if  the  material  landed  directly  in  the  water,  much  more 
time  would  be  involved  in  transportation  In  streams  and  through  reservoirs  to  the 
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intake  of  the  water  wOi’Ks.  There  would  be  further  delay  la  getting  hnto  streaas 
thooe  particles  deposited  Initially  on  diy  surfaces.  Another  delay  would  be  In 
the  occurance  of  rainfall  of  sufficient  intensity  to  dissolve  and/or  wash  this 
naterlnl  into  streams.  All  of  this  time  will  alio,  radioactive  decay  to  proceed 
thus  ameliorating  the  hazardous  effects  of  the  fallout. 
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I.  Decontamination  of  Water  Supplier* 

1.  Removal  by  Conventional  Wp.ter  Treatment  Processes 

Conventional  municipal  water  purification  processes  generally  Include 

aeration,  chemical  coagulation  with  sedimentation ,  rapid  sand  filtration  and 

(79} 

chlorination'  .  To  a  lesser  extent,  lime-cods,  osh  softening  and  :,lov  rand 
filtration  are  used.  Fxcept  for  c-eintlon,  the  processes  are  capable  of  removing 
radioactive  contaninatlon  to  some  decree,  either  singularly  or  in  combination. 

The  decontamination  capability  of  each  type  of  process  is  discussed  separately 
below,  with  particular  reference  to  the  six  elements,  barlun,  cesium.  Iodine, 
lanthanum,  ruthenium  and  strontium. 

a.  Chemical  Coagulation 

The  most  common  coagulants  in  water  treatment,  aluminum  and  Iron  salt*, 
form  aluminum  or  ferric  hydroxides  which  precipitate  as  floe.  The  chemical 
floe  acts  as  an  efficient  scavenger  by  ndcorblng,  entrapping  or  otherwise 
bringing  together  suspended  matter,  particularly  that  which  Is  colloidal  in 
nature.  The  artificial  increase  of  the  alkalinity  in  water  may  also  form  the 
hydroxides  of  heavy  metals,  which  co-precipitate  with  the  aluminum  or  ferric 
hydroxide. 


ORNL 
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Coagulation,  followed  by  sedimentation,  has  been  extensively  studied  by 
with  rcp-'i-tcd  removals  of  only  36$  and  51$  for  Cs  and  Sr,  respectively. 


It  appears  that  the  method  Is  effective  for  the  removal  of  suspended  or  colloidal 
material  ar.d  for  most  cations  of  valence  3,  •*  or  5,  including  the  rarc-earths 

group.  Removals  In  excess  of  98$  have  been  reported  for  r- 32  as  the 

.  .  ,  <8l)(82) 

orthophosphate'  '  . 

(83) 

Removal  of  fa,  I  and  La  by  coagulation  was  investigated  by  Lacy  , 
with  a  maximum  of  70.7$  found  for  both  Ba  and  La,  but  only  ^5$  for  I. 


References  begin  on  page  1J5 
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Matsuraura^®^ 

coagulation. 


states  that  Ru  can  be  removed  up  to  JA  by  chemical 


b.  Rapid  Sand  Filtration 

The  amount  of  radioactivity  removed  by  filtration  vlll  vf.ry  dep»  r.ili'.i!  ur. 
the  nature  of  the  material.  ^  cad  Zr,  probably  present  in  the  colloidal  ~ta*e, 
veie  removed  from  up  to  935“  by  sand  filtration  alone  as  reported  in  a  com- 
preher.s lve  C.RTIL  report^®0^,  whereae  materials  such  as  Sr  )  and  C~  (’r-Oi) 
present  in  true  solution,  were  cot  greatly  reduced  by  passage  through  siund 
filters.  They  further  reported  that  Pa  and  la  could  be  removed  up  to  sv.d 
iH,  respectively . 

No  data  were  reported  on  the  removal  of  !  or  Ri  by  rapid  suad 
flit  ration. 


c.  Chlorination 

Hannah,  et  studied  various  methods  for  the  removal  of  1*131 

from  water,  and  found  that  email  dosages  of  chlorine  in  the  presence  of  100  ppm 
activated  carton  produced  up  to  8056  removal.  The  authors^®'’ ^  concluded  that 
the  only  effective  method  found  for  removing  1-131  with  materials  normally 
available  in  water  treatment  plants  involves  chlorination  followed  by  adsorption 
of  liberated  iodine  on  activated  carbon.  The  optimum  chlorite  dosages  were 
quite  small,  1 ,c . ,  0.05  to  0.1  ppm.  The  removal  of  1-131  decreases  to  less 
than  ZOf>  when  the  dosage  of  chLorir.e  increases  to  ippra.  Therefore,  normal 
pre -chlorination  employed  by  the  water  treatment  plant  could  not  be  used  for 
lodlr.e  removal  because  the  chlorine  residuals  would  generally  exceed  i  no  re¬ 
quired  dosage  for  activity  removal.  Static  iodine,  in  dosages  greater  than 
0.01  ppm,  Inhibited  removal  of  1-131  with  chlorine  ar.d  Ajua  liuchiar  A 
(activated  charcoal).  Variation  of  pH,  achieved  by  adding  sulfuric  acid  and 
rodlun  nydroxlde,  was  found  to  have  little  effect  on  the  removal  of  iodine 
by  chlorine  and  activated  charcoal. 


d.  Lime-Soda  Ash  Softening 

Softening  plants  remove  excessive  amounts  of  scale-forming,  soap- 

consuming  compounds,  chiefly  comprised  of  th*  cations  of  calcium  and  magnesium, 

by  the  addition  of  lime  and  soda  ash  which  precipitates  calcium  as  a  carbonate 

and  magnesium  as  a  hydrate,  Because  of  the  similarity  In  the  chemical  picperties 

between  strontium  and  calcium,  the  former  ion  is  co-preclpitated.  For  most 

h  5 

satisfactory  removals  of  strontium,  Poyl  found  that  excess  dosages  of  both 

,  (g„  ) 

lime  and  soda  ash  were  required.  Removal  efficiencies  of  99*7?  are 
possible  under  favorable  conditions.  Other  Isotopes  which  can  be  effectively 
removed  by  the  lime-soda  ash  are:  3a,  Cd,  f ,  Sc,  Zr  and  Nb.  The  removal  of  Ba 

/p-  \  /r  a  \ 

(99$)  and  La  (90$)  was  reported  by  McCauley,  et  al  1  ,  0FC.L  reported 

that  Cs  could  be  removed  by  lime-soda  softening,  tut  not  efficiently  (25$). 

e.  Slew  Sand  Filtration 

Downing,  et  &1  reported  that  the  slow  ^and  filter  was  very 

effective  in  the  removal  of  radiostrontiutt  in  the  first  few  days  of  operation, 
but  that  the  rff  iciency  decreesed  rapidly  and  essentially  drops  to  zero  at 
days.  In  another  paper,  presented  by  Eden,  et  al  ,  on  experiments  carried 
out  in  the  sane  Water  Pollution  Research  Laboratory  In  Britain,  the  slow  sand 
filter  appeared  to  be  effective  for  only  a  few  hours.  The  activity  of  the 
effluent  rose  steadily  and  reached  30$  of  the  initial  (unflltered)  water  in 
one  day,  70$  in  2  days,  and  95$  after  7  days.  Eden,  et  nl  ^  also  reported 
that  iodine  was  not  efficiently  (50$)  removed  by  the  slow  sand  filter. 

2.  Removal  by  Non-Convent  lonal  Treatment  Methods 

Below  are  listed  several  methods  of  decontamination  of  water  not 
commonly  employed  in  municipal  treatment  plants.  Economic  considerations  nay 
limit  the  use  of  many  of  these  methods,  except  in  emergency  cases, 
a.  Ion  Exchange 

A  few  water  treatment  plants  use  loc  exchange  as  a  softening  measure. 
Resins  have  been  found  to  provide  one  of  tbs  cost  effective  methods  of 


decent tni nation.  Over  r/)-9f>  pc-bcvsI  of  3a,  C&,  Cc,  Cs,  I,  La  ead  over  99^ 

reaeval  of  V  and  Y,r  vuo  i4uotad  by  OKru  ^  |  ar«l  Airphlett  and  S&mmoa  ^  ^ 

U  s 

reported  vork  using  nix-d  bed**  A  derontamiration  factor  Of  10  to  10'  may  be 

fe-  ) 

achieved  if  Indulge  ie  «.  j  i  situated  .  The  natural  greens  an  da,  used  In 
tome  municipal  coitculrtg  Mietcllaticne,  nre  cot  effective  La  the  removal  of 
anions. 

Although  the  Ion  oxtfi&uge  process  offers  one  of  the  isost  promising 

methods  for  rc-jcvlng  tv  1  hi-contaminanta,  the  coat  may  preclude  widespread 

application  unless  sulLet lo  regeneration  techniques  ore  developed. 

The  po.  sii.-mty  of  using  hone  vatcr  softeners,  of  the  ioa  exchange  type, 

A>?  5 

has  teen  discurru-d  in  th-  literature  v’’ 
b»  Phosphate  Q.vij’-ilatloa 

Seca-use  the  conventional  method  employing  alum  or  Iron  salt  as  the 

(o'*  ) 

coagulant  does  not  affectively  remove  Sr  and  many  other  Isotopes,  Lauderdale 

Invest igated  phosphate  o?(lfi'1l«tton  as  a  mean*  of  increasing  removals.  A  more 

iy  ) 

comprehensive  study  yoe  ir*»'!v  later  by  h’esbitt,  et  ai  ^  .  Based  oa  the 

f  s' 

theoretical  conslderatlonf*,  ftfl  veil  05  experimental  observations,  Nesbitt  ?y 

concluded  that  the  me.d.or  1  t*mn  of  the  removal  of  rudloactivity  arc  ccprccipitatlon 

and  adsorption.  Ke.-.bitt  C'  p.'.’tcd  a  strontium  removal.  Udie  three  variables 

vhlch  wort  found  t.c  exert  Mi®  greatest  influence  on  cocgHstion  «re  pH, 

calcium  cor. cent r-xt ion,  w.d  j'ho'iphr.t.r  concentration.  Careful  control  of  there 

vuriablto  vac  necessary  t.c*  achieve  optimum  strontium  rcaovaln. 

(yi) 

Tan’ira  and  StruMnf b  reported  that  low  cost  natural  phosphates 

could  be  uroJ  in  plane  of  phosphate.  The  authors  reported  that  over 

yCji  removals  vo re  (•chic*.'-.’'1  Wth  these  minerals. 

Removal  of  (Jv.?'/,  fVi  :-n,  r/}f>  for  I/i  rr.d  for  Ru  have  been  reported 
by  Katcunrum  \  Cow-rel-  a-  ^  ^  report  only  yrft  removal  of  Cs  using 


phouptiate  coagulation. 


At  leaut  two  modified  coagulation  processes  have  teen  studied,  o«  of 
„H‘h  Is  the  continuous  addl.tion  of  city.  One  modification  of  the  piccese  la 
the  addition  of  small  amounts  of  silver  rltiate  to  increase  the  rcmovil  of 
lc-11dc.  The  removal  of  95-5$  1-131  has  teen  reported  by  Eden,  et  aly  . 

Lacy  and  de  lagu"'  rcpoi’t  rcocvnlc  of  9?-7$  for  Ba  cr.d  Ce,  and 

.97) 

yy,li$  for  Su  when  coagulation  is  supplemented  by  cloy  olditlon.  Cowsor,  .et  al 
u.ijr.g  a  llme-c'oda  ash  softening  process  with  the  addition  of  ciuy,  reported 
t  .t'/Kenlum  removals  of  76$-  The  removal  of  cesium  averaged  86$  with  clay  doses 
of  "'00  ppm. 

4.  Metallic  Dusts 

lacy'1'1  '  conducted  laboratory  Jar-teet  studies  for  the  removal  of 
radioactive  contaminants  by  employing  various  concentrations  of  pcnnlered 
nl' /Tilr.um,  copper,  iron  and  sine.  la  general,  highest  removals  were  obtained 
vlth  Iron  powder.  Removals  of  37-2$  of  Z  ve.re  reported  In  19&,  a 

pi.t  ciit  van  grented  to  Sllker  ^for  removing  ?.  As,  Kn,  the  rare  earth  metals, 
and  e.ctlr.ides  from  aqueous  eolations  by  sorption  on  particles  of  aluminum  metal. 

e.  Olay  Materials 

Straub,  et  al  and  others  vl0tv  .1C1*)  the  .^c  of  various 

t.ypeu  of  obeys  (kaolinite,  montmoi'illor.ltc  nrid  shales)  for  the  removal  of 
rpu'-lfi'.’  radioisotopes  from  water  and  vacte  solutions.  The  clays  vere  found 
to  he  especially  suitable  for  the  removal  of  Cs,  7. r  and  lib.  With  the  addition 

/  Qa\ 

of  , COO  ppm  of  clay,  removal  of  over  9&$  of  t’s  vao  report’d  '  v  .  Amphlctt  ar.d 
S.-jiiium  ^  reported  93$  removal  of  Sr.  This  method  r.ot  be  economically 
j  u  t ll'icd  because  large  volumes  of  clay  must  be  handled  both  initially  and  as 
c.  ■•.i;.huni:mtcd  sludge.  T':ii3  method  has  its  greatest  potential  use  where 
gr  >*'iglc  ar.d  hydrologic  conditions  arc  such  that  the  ground  itself  may  be 
v»i  for  disposal  of  the  eludge.  Studies  have  also  beet  made  on  the 

ujo  of  clay  for  high  level  vaetc  dispoc:d  by  fixing  specific  radioisotopes 
t„'  o  clay  by  means  of  .lgn  temperat-ir-  trevtscent . 
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f •  Oxidation-Deduction 

This  method  was  especially  designed  to  improve  !.u  removal.  Runhenlum 
Is  often  present  in  the  fora  of  a  nitrosyl-ruthenium  conplex  which  ia  difficult 
to  remove  ty  coagulation,  procedures.  An  approach  to  the  problem,  described 

/  qjA 

by  Amphlett  and  Samson  '  ,  involves  pre-treatment  v  ith  an  oxidising  agent  to 

destroy  the  ruthenium  complexes,  iolio.ci  uj  a  calcium  phoephate-ferrouu 


phosphate  floe  in  the  presence  of  to  maintain  she  ruthenium  in  a 

reduced  form.  As  high  as  99^  of  ruthenium  was  removed  when  100  ppa  KaHSO^  was 
used. 


g.  Foaming 

A  foam  separation  technique,  based  on  surface  activity  phenomena,  waa 

recently  developed  by  Radiation  Appli  eat  tons,  Inc.,  Long  Iolund  City,  Lew  Yorh^0^. 

A  reduction  of  the  concentration  of  radioactive  strontium-yO  and  cesium-137 
-5  -12 

from  10  to  10  molar  in  lav-level  nuclear  waste  was  claimed.  The  concen¬ 
tration  of  the  foaming  agent  is  one  of  the  most  critical  factors  in  determining 
the  separation  factor  (surface  adsorption  factor).  A  detailed  discussion  of 
the  theory  as  we  LI  ao  experimental  observations  of  the  foam  separation  technique 
was  presented  by  Schoen  and  Mazzella  '  . 

Laqy'i ^  Investigated  the  flotition  process  in  1957  and  reported 
removals  of  up  to  78$  for  oa-lbO-La-lbO,  8[#  for  Cs-137>  7*+$  for  Ru  and  68 In  for 
Sr-90.  These  studies  were  conducted  using  cetylpyridinlum  chloride  bg  the  gut- 
face-active  agent  and  Teaming  by  aeration. 

3.  Kncrgency  Methods  for  the  Dec  or.  Lamination  of  Radioactive  Water 

According  to  the  literature,  very  few  emergency  decontamination  unlt3 
are  now  available  for  municipal  use,  although  the  U.S.  Arsy  has  developed 
several  mobile  units  which  could  find  vide  application  for  public  use.  Several 
smaller  decontamination  unit6  are  corrercially  available,  most  of  which  employ 
ion  exchange  techniques. 


Supplementary  treatment  to  the  conventional  water  treatment  processes 


would  probably  he  the  cost  efficient  menns  o:'  supplying  potable  water  to  a 
larje  number  of  people,  following  radioactive  contamination  of  the  water. 

a.  Municipal  Gita  teconts'B.tr.fttton  Unite 

*  A:  o) 

Woodward  and  Pnbeck  v  '  report  on  ar.  ion  cjcehcrge  process  tnat  mr.y  be 
u£ti  to  supplement  the  normal  water  treatment  procedure.  An  Ion  exchange 
cclvnui  Le  the  form  of  a  cartridge  is  inserted  into  .10  system,  and  i dl.: paced 
of  after  'n-eriX'/hroo/h  occurs.  The  type  cf  rer In  v-lli  vary  according  to  the 
time  after  detonation  that  the  water  will  bo  p-t  into  use.  A  mlx>- 1-bod  :-.-cnin 
’•til  be  necessary  for  irmrediate  u.-o  of  tbo  voter,  whereas  a  notion  rvslu  will 
.'j'fl'.e  for  long-term  usage,  due  to  the  decay  cC  short-lived  iodine-131. 

Tee  life  of  t-a:h  lor.  cartridge  will  deptri  upon  ',‘c  total  sciido  present. 

Although  large,  cent  rally  lo-uhed,  da  font  xninat  ion  ur.it:  an:  more 
c:jr.omical,  smaller  decentralized  units  would  b?  preferable  os  the  problem 
cf  water  diatributioa  .xmld  be  lessened. 

b .  hie  Id  Decontamination  Tilts 

fti  1)  dlfl 

JThe  1 .3.  A r*v  has  developed  e:\v~rnl  mobile  decontamination 

utits  prlnm  lly  for  use  by  troops  in  Mv*  field,  but  these  could  also  be  -.mr.i  to 
s.pply  water  for  small  population  group:  in  case  of  nuclear  attack.  One  such 
unit  is  comprluod  of  a  flocculator,  filters,  dual-bed  Ion  exchange  column 
and  a  cblc-rlnnfcoi,  'ihe  entire  unit  is  mounted  or.  two  2-1  fi  ton  trucks.  A 
1,500  gpfc  output  la  obtainable,  regeneration  ol  the  ion  exchange  l-eclno  la 
effected  by  hydrochloric  acid  and  soda  ash. 

c .  Mlrc-llftntons  lei  or.  bend  nation  'Jr.lt 3 

laaotts*  Ohenlccl  Company  ha:  Ic.'eiOpvi  a  mixed  ion  exchange  jvcfn 


f 
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vater.  »e  -jolt  coeulrif  of  a  ted  of  Ariterlit-  realn  !3-3  ClH-120  and  J/'A-kiaJ 
and  la  capable  of  producing  10  gallons  of  vater  of  triple  distilled  quality. 


A  -f.lt  dt.-.e loped  by  Lauderdale  and  Etcori 


consists  of  tv! (i 


columns  of  steel  root,,  burnt  clay  cxd  activated  carbon  In  one  col — a  a:.<1  a 
ul.-i-l  bed  ion  exchange  resin  in  the  other.  Approximate!;.'  jO  Ute.-o  cl'  outer 
can  be  dei:c>>A«.,<no*«^  Vith  as  exceeding  99.9$. 


N..-CSO 


<5ii.  lA-y  ~  r\tzz  rite  e  rr/^o-J  cf  rrtHo^ciive  vntci'  >i<?con- 


taainatlon  whereby  the  water  is  passed  through  a  column  of  natural.  rv.Wtalo, 

.-,ueh  as  clay,  leaves,  huff'^,  gravel  and  sand.  Pesovals  of  over  9^  were 
reported.  O.e  problem  Wild  be  to  obtain  ur.cor.tuainate-d  materials  from  the 
environment  for  use  In  ouch  a  col-nan. 

** •  Diaeussian  and  Cone ImlonB 

— e  more  important  radioisotopes  that  appear  initially  it  ccritOfiMated 
vc.ter  are  Ba,  Cs,  I,  La,  Ru  ar.d  Sr.  T.ot>oratoiy  data  indicate  that  the 
combination  of  conventional  water  treatment  nvithods,  nt-t.’ly,  chemical  coagu¬ 
lation,  lime -coda  ash  soft  ning,  and  sued  filtration,  Is  capable  cf  rcw>/lng 
59$  ctrontium.  The  rctayval  of  iodine  con  be  accomplished  by  the  addition 
of  silver  ions.  86$  cesium  can  be  renu>vcJ  by  the  llme-.-oda  ash  softening  procedure 
when  proper  amour,  to  of  clay  are  added. 

'fhe  efficiency  of  gross  activity  removal  lice  generally  tetve'.'i  '/)  and 
75$  in  plant  scale  cporntlo:.***^.  The  lower  plant  cce.ie  efflcienolcc,  In 
comparison  to  the  Irtborntoiy  Jesuits,  are  bvll-.-vod  to  result  from  the  rant 
t.lvit;  (a)  only  one  or  two  decontamination  procc-ii;en ,  .-jsccg  those  Ber.t'ovid, 
tre  employed  by  the  ur.nnl  vater  t rent  sunt  plart,  aul  t.hc  operating  ci/vU*. Irmti 
m.-\y  not  be  optimised;  aj-.l  (b )  in  the  laboratory,  the  radioactive  rsitarial  >  arc 
generally  present  lr.  the  form  of  simple  salts.  In  contrast,  the  rad  lot."  Uve 


*15'' • 


materiad-i  reaching  the  vcter  treatment  pit et  are  generally  believed  to  be  ir. 

ehesieel  forms  which  are  mors  difficult  to  remove  beeducc  the  portion  which  ifl 

readily  remove!  has  already  been  eliminated  by  natural  decor. taetnat ion 

(llA 

mechanisms  before  reaching  the  water  plant.  Setter  and  .boa 3 ell''  stated  that 

the  longer  the  activity  has  been  subjected  to  natural  purification,  the  more 

difficult  it  is  to  remove  the  remaining  soluble  activity  by  coagulation. 

~b  rho'sld  b-c  pointed  out  here  that  the  decor. lamination  data  reported 

ir  the  literature,  both  on  laboratory  and  plant  scale,  ar«  based  on  or.o  of  the 

follcviug  three  scarcer  of  radioactivity: 

(l)  added  radioactive  salts 
(2J  lev-level  radioactivity  waste 

(3)  long-range  (vorl’.-vide)  fallout  from  atomic  bomb  teats 
Che  chemical  and  physical  characterlotlea  of  radian  chive  materials  fw.t 
any  of  the  above  sources  is  different  from  that  proceut  ir.  local  fallout.  The 
latter  is  characterized  by  Its  lor  water  solubili  ty  of  1$  to  23C  in  comparison 
to  the  or  higher  solubility  of  radioactive  materials  cf  the  above  three 
categories.  Low  water  solubility  generally  will  incresr.e  the  removal  of  local 
fallout  by  conventional  water  treatment  mtbods  which  ara  efficient  in  dealing 
with  particulate  or  colloid  materials. 

The  atxUnuin  decontamination  factor  for  etch  of  six  selected  elements 
above  ere  listed  in  Table  X.  Thcce  values  arc  baced  on  the  best  available  data 
in  the  literature.  It  is  emphanieed  that  the  percentage  of  removal  Is  greatly 
dependent  on  the  chealcnl  form  and  physical  state  of  the  radioactive  element  as 
well  2.3  the  concentration  of  the  treatment  additives  and  other  related  conditions 
such  ar-  p'i  ml  temperature  of  the  water.  As  an  example  of  the  decontamination 


factors  required  to  reduce  activity  to  the  Mc’J,  based  on  hypothetical 
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Tnble  X 


Mnxlmun  Decontamination  Kac‘ors 

Peported 

in  the  Site  rv 

-r?  for 

Selected  Iso'-rc* 

Decontaai*. 
nut  too. 

Method 

Decontamination  Factors  (D.F.)  and  Pe-f 
B'i-ltO  Cs-137  1-131  fa-lltO 

*  recces 

106 

Sr-89,90 

1.  Che-raical 

3-4 

1-5 

1.8 

3.4 

10  c 

2.0 

Coagulation 

(83) 

(80) 

(83) 

(8f) 

(84) 

(So' 

2.  Llme-3od» 
Softening 

ioa 

(87) 

(83) 

m 

10l 

(87) 

ifR 

3-3  x  i;L 
(So) 

3.  Rapid  Sand 

Pi It ration 

2  x  101 
(80) 

(30) 

m 

3-8 

(80) 

SP. 

1.04 

(So) 

4.  Ion 

10* 

10* 

103 

103 

ICR 

3.3  x  1:3 

Exchange 

(80) 

(80) 

(9i) 

(80) 

(80) 

3 •  Phosphate 

3.2 

1.5 

NF 

5  x  101 

■c-7 

4.5  t  i?L 

Coagulation 

(84) 

(9?) 

(84) 

(34) 

(94) 

6.  Cloy 

ifR  5 

X  10* 

(fio) 

tf) 

m 

ICR 

2  x  101 
(9U 

7-  Metallic 

Dusts 

NR 

hr 

1.6 

(80) 

NR 

ICR 

HR 

8-  Coagulation  v / 
Silver  Ions 

jrR 

HR 

2.2 

(96) 

HR 

SR 

NR 

9*  Oxidation- 
Reduction 

im 

HR 

NR 

HR 

102 

(91) 

HR 

10.  Coagulation 
with  Clay 

3.3  x  10l 
(99) 

3-3  x  iO1 
(99) 

'  HR 

nr 

3.3 

(99) 

HR 

It.  Lime  Soda 
with  CLay 

nr 

7.1 

(97) 

IfR 

NR 

4.2 

(97) 

HR 

l? .  Foaming 

4. 5 

6.7 

HR 

•*•5 

3.9 

3.1 

(Flotation) 

(109) 

(109) 

U09) 

(109) 

(109) 

13.  Glow  Sand 
Filtration 

NR 

HR 

(90) 

IfR 

SB 

HR 

Notes; 

(1)  Percent  Removal  »  100  .  (gi--)ioo 

(2)  RtJ  -  not  reported  iu  literature  revieved 


.  »  1 


■  'stji  %j> ' 


of  the  Springfield,  Haas.  vatershcd  (aa  reported  eloevhera  la 
ir.it  rrpiirO,  Table  XT  vna  prepared.  It  la  apparent  from  thin  tools  that  ion 
Grange  u.  the  only  £l_^le  decontamination  procesa  that  will  adequately  remove 
■<v*  alnlsna.  concent  rati  on  reported  for  the  six  celeettd  Isotopes*  and  that  no 
pro-'eaa  .111  adequately  remove  the  taxlcroa  radlocontaminant,  namely  I-I.31. 


-195 


AFFKNpIX  B 
Distribution  [lot 


rro.  of 

Adjvessce  Copies 

rfi'ict  of  Civil  lefease,  DCC  38 

Washington,  D.  0. 

A' 1  *-■'» :  Director  for  Research 

A*’"'  .Utrury,  Civil  Defense  Unit  3 
pentagon,  Washington,  D.  c. 

Assistant  Secretary  of  the  Array  1 

(rt  &  n) 

V/poi-.lr.gton  25,  U.  C. 

At1'-'*:  Assistant  for  Research. 

(ji.  i  <-f  of  ItHval  Research  1 

(t'tf'H*  IOU) 

j>.  pur  Ur.  vat  of  the  Navy 
Washington  25,  D.  C. 

Cl-i uf  of  Unvul  Operations  1 

(ifrorno) 

jv.jv.j-tment  of  the  f.avy 
Wayiiii.gton  25,  D.  C. 

Chief#  Bureau  of  l! aval  Weapons  1 

IUv<K-5) 

C«.fte.i't,rr.ent  of  the  -;avy 
We-r. tr.gton  25,  D.  0. 

Chiff,  Bureau  of  Medicine  and 

swgcry  ^ 

he^x-rtment  of  the  Davy 
WSniiir.gtori  25,  D.  C. 

Ch! <uf,  I.ur'-eu  of  Supplies  and 

/' ‘.count  0  ^ 

(•oidc  L12) 

1  ament  of  the  Navy 
!  I'.gf  Oil  f'5,  T:.  '.* . 

Chi'-f,  Bureau  of  Yards  &  Docks  1 

C trice  of  Research  (Code  7!) 
r>.  *  meat  of  the  Navy 

Woci  'ngt on  25,  D.  C. 

U,  »S.  [.'aval  Civil  Engineering  1 

;/-i  oratory 

pc,  fivnTe,  California 


Ko.  of 

addresses  Copies 

Advisory  Committee  on 

Civil  Defense  ^ 

riationtl  Academy  ;.f  Sciences 
2101  Constitution  A.o .  X,  w, 
Washington  25,  D.  c. 

Atten:  Mr.  Richard  Park 

Defense  Document-. tion  Center  20 
Arlington  Kf.ii  Station 
Arlington,  Virginia 

Chief,  Bureau  of  Ships  2 

Department  of  the  Ya/y 
(Code  33?) 

Chief,  Bureau  of  Ships  1 

(Code  203) 

Department  of  the  Navy 

Principal  investigator  1 

Office  of  Civil  Defense 
Contract  QCP-CS- 62-205 

General  Dynamics  Ccip. 

P.  0.  Box  5 

Old  San  DIcgo  Station 

San  Diego  10,  California 

Principal  lnveetigetor  1 

Office  of  Civil  Defense 
Contract  CCD-CS-62-166 
Technical  Operations,  Inc. 
Burlington,  Massachusetts 

Frir.cipr.1  Investigator  1 

Office  of  Civil  Defense 
Contract;  O:0-CS-62-135 
Stanford  Runv-.-.reh  institute 
Menlo  Park,  California 

Principal  Investigator  1 

Office  of  Civil  Defense 
Contract;  tCD-'JS-OV-ltS 
Naval  Civil  Engineering  Lab. 

Port  Iiuerx-.Tc ,  California 


Addressee 

Principal  Investigator 
Cfflce  of  Civil  Defense 
Contract;  'XD-CS*62-193 

Isosores,  Incorporated 
ID  3  Wood lawn  Avenue 
Westwood,  New  Jersey 

Principal  Investigator 
f : ice  of  Civil  Dcrense 
lent  root:  CDM-SR-59-51* 

;,val  Philological  Defense 
laboratory 

Can  Francirco,  California 


Irir.cip.i  Investigator 
y.fice  of  civil  Defense  Research 
Army  Chemical  Corps 
Deps.rtment  of  the  Army 
Ite-mlcal  Center,  Maryland 


No.  of 
Coplea 


Addressee 


No.  of 
•'oplcu 


Principal  Investigator 
Office  of  Civil  Defense 
Contract :  XD-C6-6C-279 

Curt  its -Wright ,  Tncorpo rated 
Caldwell,  ilc-v  Jersey 

Principal  Investigator 
Office  of  Civil  r c Dense 
Contract :  •  •O.'-cC -037 

iiorf  FIT  l.'.oru'.orle; ,  Ir.c. 
1922  East  107  Street. 
Cleveland,  Thio 

Prir.cip.l  Investigator 
Of  it co  of  Civil  reverse 
Contract;  CC  D  -  IS  -  6?  - 106 
Engineering  Science,  Inc. 
ItO  E.  foothill  Elvd. 
Arcadia,  Calif. 


Principal  Investigator  1 

Du  ice  of  Civil  Defense 

C  or.l  ract :  0CD*0S*62  *206 

lor. let ,  Incorporated 

132  -  bth  Street 

c>.  abridge,  Massachusetts 

Principal  Investigator  1 

If  lice  of  Civil  Defense 

Contract:  wCD»Ci% 62*73 

..  S.  Army  Material  Office 
Directorate  of  Research  and 
Development 
Weshlr.gtor.  25,  D.  C. 

r rlrcipal  Investigator  l 

II  l  ice  of  Civil  Defense 
Cornell  Aeronautical  Laboratory 
U<*55  Ifrr.esee  Street 

iulfalo  21,  Mew  York 

•  r  urcipal  investigator  1 

III  ice  of  Civil  Defense 
Cor.ric.et:  CCD-0S-62-ll<4 
:-.;arcti  Triangle  Institute 

c.  Pox  U90 

Durham,  North  Carolina 


Principal  Investigator  1 

Office  of  civil  Defense 
Contract;  IvD  VC'S  *62  *202 
Armour  Research  Foundation 
10  West  3;th  Street 
Chicago  16,  Illinois 

Homer  I.  McConnell,  Director  1 

Civil  Ceter.se  Training  Prcgrm 

fda/dhew 

Washington  2?»  D.  C. 

Lt.  Col.  Converse  P..  Lewis  1 

Office  of  the  Surgeon  Ccneral 
Main  Navy  building 
Washington  2?»  D«  C. 

Dr.  Lur.a  I-eopold,  Chief  1 

Water  Resources  Division 
v.  o.  Ceoioglcal  Surrey 
Washington  2j,  D.  u. 

Cordon  E.  Mela Hum,  Chief  1 

Water  Supply  *  lollutlon 
Control 

Dept .  oi  Health,  Education  and 

Welfare 

Public  Health  Service 
Washlrg-.cn  2$,  D.  0. 


-197- 


Addressee 

G« orgc  W.  Burke,  Jr. 

Dlv.  of  Water  Supply  and 
Control 

Publl^  Health  Service 
Wash  lepton  25,  D.  C. 

Fro.  THonas  deS.  Fun  urn 
Civil  Engineering  Dept. 
University  of  Florida 
Gainesville,  Florida 


Ko.  of 
Copies 

1 

Pollution 


1 


c'O.  of 

Addressee  Copies 

If.  Mai con  Childers  1 

70s  e'orth  West  Street 
Alexandria  la,  Virginia 

Dr.  Edward  Ray,  Director  1 

Kuclear  Science  &  Engineering 
Corp. 

P.  0.  Box  10901 
Pittsburgh  3^,  Pennsylvania 


